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ABSTRACT 
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•>  a 


subsonic  wind  tunnel  investigation  was  conducted  at  the  Naval, Postgraduate 
School  Wind  Tunnel  Test  Facility  to  examine  the  effects  of  grid  generated  turbulence 
on  the  high  angle  of  attack  aerodynamics  of  a  vertically  launched  missile.  Four  turbu¬ 
lence  generating  grids  produced  turbulence  with  length  settle  to  missile  model  diameter 
ratios  of  1.05,  0.89,  0.62,  and  0.15  with  respective  turbulence  intensities  of  3.3 1%. 
2.78%,  1.88%,  and  0.47%.  The  test  model  is  representative  of  current  cruciform  wing 
missiles  with  low  aspect  ratio  wings  and  long  root  chords.  The  tangent  ogive  nose  has 
a  nose  Fineness  ratio  of  2.29.  Three  separate  body  configurations  were  tested  with  and 
without  the  turbulence  generating  grids.  One  configuration  was  a  body  in  isolation  and 
the  other  two  were  wing  and  tail  configurations  at  0^  and  45c  roll  angles.  All  test  runs 
were  conducted  at  over  an  angle  of  attack  range  of  -56  to  95^.  Results  in¬ 

dicate  that  as  turbulence  length  scales  approach  body  diameter  size,  the  angle  of  attaek 
for  the  onset  of  asymmetric  vortices  was  delayed  and  the  side  force  magnitude  was  re¬ 
duced.  The  vortices  generated  by  the  nose  of  the  missile  continue  to  dominate  the 
afterbody  vortices  for  body  configurations  with  and  without  wings  regardless  of  the 
turbulence  conditions.  p  J  (~ 
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I.  INTRODUCTION 


A.  BACKGROUND 

The  development  of  missiles  capable  of  operating  at  high  angles  of  attack  has 
prompted  numerous  investigations  and  experiments  in  high  angle  of  attack  research. 
Much  interest  is  due  to  the  demands  for  highly  maneuverable  missiles  at  supersonic 
speeds.  In  cases  where  the  missile  is  launched  fiom  the  ground  or  at  sea  from  a  ship,  the 
missile's  initial  velocity  for  all  practical  purposes  is  zero.  From  this  initial  state,  the 
missile  must  fly  over  a  wide  range  of  Mach  numbers  from  subsonic  to  supersonic  speeds. 
The  aerodynamic  characteristics  of  the  missile  may  vary  greatly  as  it  transitions  through 
these  operating  speeds.  All  surface-to-air  missiles  follow  a  three  phase  operation: 
launch,  nudeourse.  and  terminal.  Each  phase  represents  a  different  missile  flight  regime 
owing  to  the  changing  flight  control  requirements  and  airflow  about  the  missile,  thus 
presenting  a  challenge  to  the  missile  designer  who  must  develop  a  missile  capable  of 
operating  throughout  each  regime. 

Vertical  launch  for  missiles  is  by  no  means  a  novel  idea.  It  has  been  the  primary 
means  of  launching  large  rockets  and  missiles  from  land.  Vertical  launch  of  surface-to- 
air  missiles  at  sea  is  a  recent  development.  Conventional  launch  systems  utilize  a 
trainable  launcher  consisting  of  a  rail  or  tube  adjustable  in  azimuth  and  elevation.  The 
vertical  launch  surface-to-air  missile  (VLSAM  f  uses  a  canister  container  which  stores 
and  launches  the  missile.  The  canister  system  eliminates  the  cumbersome  trainable 
launcher  and  associated  specialized  storage  facilities  and  loading  mechanisms.  Other 
advantages  of  the  VLSAM  are  the  360°  coverage  provided  without  interference  from  the 
ship's  superstructure,  a  higher  firing  rate  since  launcher  availability  is  no  longer  a  factor, 
and  space  and  weight  savings  [Ref.  1[. 

VLSAM s  have  their  own  unique  trajectories,  allowing  them  to  point  at  their  targets 
after  launch.  Two  flight  control  methods  are  currently  used:  thrust  vectored  control,  and 
aerodynamic  surface  control  [Ref.  2|.  Thrust  vectored  control  is  used  with  short  range 
missiles  which  must  pitch  over  quickly  to  engage  targets  at  close  range.  Long  range  mis¬ 
siles  maneuver  with  aerodynamic  control  surfaces.  In  each  case,  as  the  VLSAM  exits  the 
canister  it  enters  the  open  ocean  environment  at  low  velocity  and  is  subject  to  potentially 
significant  crosswinds.  The  resultant  of  the  missile  and  crosswind  velocities  is  a  potential 
high  angle  of  attack  flow  about  the  missile  [Ref.  3].  This  high  angle  of  attack  flow  may 


cause  asymmetric  vortices  to  form  on  the  missile  nose  and  afterbody  which  induces  side 
forces  as  discussed  later  in  this  chapter.  These  side  forces  could  give  lhght  control 
problems  to  the  VLSA.V1  at  launch. 

The  launch  environment  may  contain  some  degree  of  turbulence,  both  from  tire  at¬ 
mospheric  boundary  layer,  and  the  airflow  about  the  ship.  Understanding  the  elfect  of 
this  turbulence  on  the  aerodynamic  characteristics  of  a  YLSAM  during  launch  is  the 
goal  of  a  continuing  effort  of  research  and  experiments  conducted  at  the  Naval  Post¬ 
graduate  School  (XPS).  In  earlier  work.  Roane  {Ref.  3]  developed  a  system  of  modelling 
the  flowfield  turbulence  by  a  series  of  four  grids  to  generate  turbulence  in  the  NPS  low 
speed  wind  tunnel.  The  goal  of  this  thesis  is  to  experimentally  examine  the  effects  of  this 
turbulence  on  the  asymmetric  vortex  induced  forces  generated  on  a  YLSAM  model. 

B.  HIGH  ANGLE  OF  ATTACK  AERODYNAMICS 

Flow  separation  from  the  body.  wing,  and  tail  surfaces  is  an  important  phenomenon 
encountered  in  high  angle  of  attack  aerodynamics  since  asymmetric  vortices  may  form 
to  induce  side  forces.  Primary  factors  which  influence  flow  separation  from  the  missile 
body  are:  the  angle  of  attack,  crossflow  Reynolds  number  (Rcj  l.  nose  geometry,  and 
nose  length  to  diameter  ratio  (also  called  the  nose  fineness  ratio).  Secondary  factors 
include  roll  angle  and  roll  rate,  free  stream  turbulence,  surface  roughness,  acoustic  en¬ 
vironment  (pressure  waves),  and  vibration  {Refs.  4,  5.  6].  This  section  discusses  how 
the  factors  mentioned  above  influence  the  mechanisms  responsible  for  vortex  generation 
and  disposition. 

1.  Vortex  Regions. 

As  a  slender  unyawed  bodv  changes  angle  of  attack  a  from  O'  to  90\  the  airflow 
over  the  missile  body  transitions  through  four  regimes  as  described  by  Ericsson  and 
Reding  {Ref.  7:  pp.  246-247].  These  regimes  {Figure  1]  are  characterized  by  the  chang¬ 
ing  vortex  patterns  and  the  resulting  induced  in-plane  (normal),  and  out-of-plane  (side) 
forces  and  moments  which  act  on  the  missile  at  zero  side-slip  and  at  high  angles  of  at¬ 
tack  2  [Refs.  2,  4,  S.  9  ].  These  regimes  are  bounded  by: 
asv  -  the  angle  of  attack  at  which  steady  symmetric  vortices  are  formed. 
aAV  -  the  angle  of  attack  where  steady  asymmetric  vortices  are  formed. 

«LV  -  the  angle  of  attack  where  unsteady  vortices  are  formed. 


1  The  crossflow  Reynolds  number  is  set  by  the  base  diameter  of  the  missile. 

2  Without  wings  and  tails. 

1 


a.  Regime  I  -  [0°  <  x  <  asv] 

The  axial  (low  dominates  and  flow  around  the  body  remains  attached  for 
pointed  noses.  If  the  nose  is  blunted,  a  separation  bubble  may  occur  (discussed  below). 

b.  Regime  II  -  [xsv  <  x  <  xAV] 

The  crossflow  begins  to  push  the  boundary  layer  to  the  leeward  side  where 
it  separates  and  a  system  of  two  or  more  pairs  of  symmetric  vortices  are  formed.  These 
vortices  originate  at  the  nose  of  the  body.  For  long  bodies,  more  pairs  of  symmetric 
vortices  are  generated  along  the  body  length.  The  number  and  strength  of  the  symmetric 
vortices  increase  with  angle  of  attack.  Symmetric  vortices  typically  start  at  asvs:50 
(Ref.  4J. 

c.  Regime  III  -  [xAV  <  «  <  xLA] 

Crossflow  begins  to  dominate,  shedding  asymmetric  vortices  which  induce 
out-of-plane  or  side  forces  on  the  body.  These  asymmetric  vortices  are  relatively  steady, 
but  may  change  from  side  to  side  at  higher  angles  of  attack  near  the  maximum  side  force 
magnitude  [Ref.  11).  In  some  cases  the  side  forces  may  exceed  the  normal  force.  The 
vortices  which  are  the  most  asymmetric  yield  the  highest  side  force  magnitude 
[Ref.  12].  This  regime  typically  starts  at  «AV^. 20°  [Ref.  4]. 

d.  Regime  IV  -  [xt>  <  x  <  9Q°] 

The  crossflow  dominates  completely  and  flow  separation  becomes  unsteady 
and  dependent  on  Reynolds  number.  At  the  critical  Reynolds  number  the  boundary 
layer  separates  then  re-attaches  as  a  turbulent  vortex  sheet  (discussed  later).  The 
boundary  layer  is  shed  as  a  random  wake  at  supercritical  Reynolds  numbers,  or  a  von 
Karman  vortex  sheet  at  subcritical  Reynolds  numbers.  Normal  force  magnitude  will 
plateau  and  side  force  magnitude  decreases  to  zero.  This  regime  has  been  observed  to 
begin  at  [Ref.  4]. 

Asymmetric  vortex  induced  side  forces  are  dependent  on  the  strength  and 
disposition  of  the  vortex  patterns  which  are,  according  to  Ericsson  and  Reding  [Ref.  7j, 
a  result  of  boundary  layer  transition  on  the  surfaces  of  the  missile.  The  analogy  most 
often  used  to  describe  flow  separation  deals  with  boundary  layer  transition  on  a  two- 
dimensional  cylinder  in  a  normal  crossflow.  Three  dimensional  effects  due  to  nose  and 
afterbody  geometry  must  also  be  understood  as  they  dictate  the  formation  of  the  vortex 
structure. 
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2.  Two  Dimensional  Crossflow. 

Airflow  oser  the  missile  body  can  he  divided  into  normal  and  axial  components. 
The  axial  flow  component  follows  along  the  missile  body  length  and  the  crossflow  is 
essentially  a  two  dimensional  flow  normal  to  a  cylinder.  The  crossflow  analogs  prosidcs 
information  for  cylinder  lift  and  drag  which  act  in  the  direction  of  the  crossflow.  Side 
forces  may  exist  at  right  angles  to  the  crossflow  depending  on  the  type  of  flow  separation 
on  opposite  sides  of  the  c>  Under.  Mechanisms  behind  boundary  layer  transition  and 
separation  provide  an  explanation  for  flow  separation  and  subsequent  asymmetric 
vortex  generation  effects. 


REYNOLDS  NUMBER 


Figure  2.  Two  dimensional  crossflow  about  a  cylinder:  [Refs.  13,  14]. 

The  crossflow  Reynolds  number  is  the  primary  factor  which  influences  the  sep¬ 
aration  point  of  the  boundary  layer.  Viscous  flow,  surface  roughness,  and  turbulence 


are  other  factors  which  influence  boundary  layer  separation.  Achenbach  classified  the 
flow  around  a  cylinder  into  lour  distinct  regions  delimated  by  dillering  tlow  separation 
and  drag  behavior  [Ref.  13].  These  regions  are  shown  in  Figure  Hand  discussed  below. 

In  the  subcritical  range,  the  boundary  layer  is  laminar  and  How  separation  oc¬ 
curs  close  to  to  the  lateral  meridian  where  $.  defined  as  the  angle  from  the  direction  of 
the  crossflow,  varies  from  S0°  to  90°.  Asymmetric  subcritical  separation  occurs  at 
$5:90°.  where  the  side  force  just  begins  to  reach  a  noticeable  magnitude. 


REYNOLDS  NUMBER 

Figure  3.  Side  force  to  normal  force  ratio:  [Ref  7:  p.  260]. 

When  the  Reynolds  number  enters  the  critical  range,  the  laminar  boundary  layer 
separates  from  the  body  at  $=590°  and  re-attaches  as  a  turbulent  boundary  layer  which 
separates  again  from  the  surface  at  $5:140°  .  Between  the  laminar  separation  and  the 
turbulent  re-attachment,  a  laminar  bubble  is  formed.  In  this  region,  the  flow  separation 
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can  easily  change  back  and  forth  from  critical  to  subcritical  for  small  changes  in 
Reynolds  number.  The  changes  occur  on  opposite  sides  of  the  cylindrical  section  and 
when  critical  separation  exists  on  one  side  of  the  body  and  subcritical  separation  on  the  I 

opposite  side,  a  large  difference  in  <p  from  opposite  transition  points  is  possible.  When 
this  critical  and  subcritical  combination  occurs,  an  asymmetry  is  observed  in  the  gener¬ 
ated  vortices.  The  pressure  difference  on  opposite  sides  of  the  lateral  meridian  will 
produce  the  side  force.  For  the  largest  difference  in  </>.  the  vortices  will  be  at  maximum 
asymmetry,  and  the  side  force  at  the  highest  magnitude  as  seen  in  Figure  3.  The 
asymmetric  vortices  may  also  become  susceptible  to  altering  back  and  forth,  hence  the 
change  in  the  direction  of  the  induced  side  forces.  Maximum  side  force  occurs  at  the 
critical  Reynolds  number,  making  it  an  important  parameter.  [Ref.  15] 

Next,  as  the  Reynolds  number  increases  into  the  fully  critical  range,  laminar 
transition  moves  towards  the  direction  of  the  crossflow.  Turbulent  separation  will  mot  e 
forward  to  <f>  <  140°  and  the  laminar  bubble  is  no  longer  formed.  Asymmetric  critical 
separation  occurs  simultaneously  on  opposite  sides  of  the  body.  The  asymmetric  vortices 
produce  nearly  equal  surface  pressures  at  the  lateral  meridian  because  the  separation 
asymmetry  only  affects  the  surface  pressures  at  4>  >  140°.  In  this  range  the  vortices  are 
ineffective  at  producing  a  significant  side  force  and  a  sudden  decrease  in  magnitude  oc¬ 
curs  as  shown  by  the  curve  in  Figure  3.  The  boundary  layer  thickens  and  viscous 
properties  of  the  fluid  are  the  dominant  factor.  Drag  reaches  a  minimum. 

Finally,  as  the  Reynolds  number  increases  further  to  supercritical  and  into  the 
transcritical  range,  the  laminar  transition  point  works  towards  6^0°  and  turbulent  sep¬ 
aration  occurs  at  100\  The  asymmetric  transcritical  separation  point  moves  towards 
the  lateral  meridian  where  the  asymmetric  vortices  once  again  produce  a  significant  side 
force.  Drag  will  continue  to  rise  and  settle  out  at  a  level  lower  than  subcritical  condi¬ 
tions.  [Refs.  7:  pp.  247-261.  13] 

The  Reynolds  number  dictates  the  greatest  influence  on  the  normal  force  and 
drag  characteristics,  especially  within  the  critical  range.  Maximum  normal  forces  on  the 
missile  should  occur  when  the  boundary  layer  flow  is  completely  laminar  or  completely 
turbulent,  because  crossflow  drag  is  high.  Surface  roughness  and  viscous  effects  are  also 
present.  Experiments  by  Lamont  [Ref.  16J  and  Wardlaw  [Ref.  17]  show  that  maximum 
asymmetric  loads  follow  the  same  behavior  as  the  normal  force,  with  the  lowest  vortex 
induced  side  forces  occurring  at  transcritical  conditions.  The  largest  lift  to  drag  ratio 
occurs  at  the  critical  Reynolds  number  where  drag  is  a  minimum.  Maximum  vortex 
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asymmetry  has  been  observed  to  occur  at  the  critical  subcritical  Reynolds  number  ex¬ 
plained  above.  Thus  the  maximum  side  force  to  normal  force  ratio  may  be  expected  to 
be  found  at  the  critical  Reynolds  number.  (Ref.  14) 

3.  Three  Dimensional  Vortices. 

The  initial  separation  point  of  the  shed  vortices  occurs  at  the  nose  of  the  slender 
body.  Many  researchers  have  observed  that  these  vortices  shed  from  the  nose  will  dom¬ 
inate  other  vortices  shed  along  the  body  length  [Refs.  18,  19,  20).  Thus,  the  nose  ge¬ 
ometry  becomes  an  important  factor  in  vortex  generation  and  disposition.  Axisymmetric 
missile  nose  geometry. falls  into  two  shapes:  cones  and  ogives,  both  pointed  and  blunt. 
Noses  are  dimensionalizcd  for  comparison  by  their  length  to  base  diameter  ratio  or 
fineness  ratio,  (1N  d)no„  . 

For  pointed  conical  and  ogive  noses,  observations  indicate  that  aAV  occurs  as  a 
function  of  the  semi-vertex  angle  0A  .  At  all  Mach  numbers,  asymmetric  vortex  shedding 
begins  when  the  the  angle  of  incidence  or  is  greater  than  the  apex  angle  of  the  nose,  where 
0 A  =  dc  for  a  conical  nose. 


«a\*2*a 


For  a  tangent  ogive  nose?,  with  length  /N  and  base  diameter  d\ 


0,x  =  tan' 


*N 

d 


4- 


{i} 


{2} 


A  close  approximation  for  slender  bodies  is: 


Blunting  the  nose  may  cause  a  second  type  of  asymmetric  vortex  shedding  to 
develop  at  the  rear  of  the  nose.  The  blunt  nose  causes  a  separation  bubble  to  form  on 
the  leeward  side  of  the  nose.  Asymmetric  vortices  originate  at  the  rear  of  the  blunted 
nose,  and  tend  to  dominate  the  induced  effects  from  the  afterbody  vortices  (Ref  21). 
The  asymmetric  vortices  generated  by  a  blunted  nose  do  not  alternate  from  side  to  side 


3  A  tangent  ogive  nose  is  constructed  from  a  constant  radius  arc  with  the  center  of  the  arc  lying 
in  the  plane  of  the  base  of  the  nose  and  some  distance  away  from  the  nose  axis. 
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as  readily  as  those  generated  a  pointed  nose.  Side  force  magnitudes  are  also  reduced. 
Afterbody  vortices  from  a  blunt  nose  missile  remain  closer  to  the  body  and  are  easily 
shed  throughout  the  missile  length  [Ref.  4).  Thus,  blunting  the  nose  will  reduce  the  ef¬ 
fects  of  vortex  asymmetry  and  the  induced  side  force. 

Nose  fineness  ratio  has  been  shown  to  have  an  effect  on  the  asymmetric  vortex 
induced  side  force  [Refs.  18,  19,  22].  In  general,  side  force  magnitude  increases  with 
larger  nose  fineness  ratio  [Ref.  18],  As  the  nose  fineness  ratio  increases,  the  nose  apex 
angle  decreases  and  the  angle  of  attack  for  the  onset  of  asymmetric  vortices  will  de¬ 
crease.  Thus,  a  missile  with  high  nose  fineness  ratio  may  become  more  susceptible  to 
induced  side  forces  at  a  lower  angle  of  attack  than  a  missile  having  a  lower  nose  fineness 
ratio.  Jorgensen  demonstrated  that  decreasing  the  nose  fineness  ratio  is  more  beneficial 
in  reducing  side  force  than  blunting  the  nose  [Ref.  IS). 

The  presence  of  the  afterbody  at  a  constant  roll  angle  has  almost  no  effect  on 
vortex  generation  from  the  pointed  nose  of  slender  bodies  [Ref.  23 J.  The  first  pair  of 
vortices  generated  by  a  pointed  nose  will  set  the  pattern  for  vortices  shed  along  the  body 
[Ref.  24].  The  afterbody  will  also  shed  asymmetric  vortices,  but  these  will  not  effect  the 
nose  vortices  generation,  although  the  afterbody  forces  will  contribute  to  the  overall 
magnitude  of  induced  forces  from  the  nose-body  combination. 

Should  the  aferbody  roll  angle  change,  with  no  sc  roll  angle  remaining  constant, 
the  afterbody  vortices  will  change  the  side  force  magnitude  and  direction  as  shown  in 
research  by  Kruse  ]Rcf.  25],  Since  the  nose  and  body  of  a  missile  do  not  rotate  inde¬ 
pendently  of  each  other,  the  vortex  system  generated  by  the  nose  should  still  dominate 
afterbody  vortices  on  a  axisymmetric  missile  configuration. 

Experimental  results  for  pointed  ogive-cylinders  with  and  without  afterbodies 
show  that  small  variations  in  roll  angle  alter  the  asymmetric  vortex  structure  and  side 
forces.  Researchers  believe  that  surface  imperfections  and  deviations  in  the  nose 
axisymmetric  geometry  are  the  sources  of  this  phenomenon.  The  direction  and  magni¬ 
tude  of  the  side  force  becomes  unpredictable,  a  detriment  to  anyone  attempting  to  gen¬ 
erate  a  model.  [Ref.  26] 

The  preceding  section  on  the  two  dimensional  flow  separation  analogy  discussed 
four  flow  separation  mechanisms.  Keener  identified  that  three  of  these  flow  separation 
mechanisms  occurred  simultaneously  along  the  length  of  a  nose  and  forcbody  at  angles 
of  attack  greater  than  20°  and  transcritical  Reynolds  number  of  Red  =  O.SxlO6  .  These 
mechanisms  were  observed  from  surface  flow  patterns  on  tangent  ogive  noses. 


9 


1.  Primary  laminar  separation  near  the  nose. 

2.  Primary  transitional  separation  at  the  nose  midsection  consisting  of  a  laminar 
bubble. 

3.  Primary  turbulent  separation  near  the  base  of  the  nose. 

The  extent  of  the  transitional  separation  region  grew  with  increasing  angles  of  attack 
and  constant  Reynolds  number  while  the  turbulent  separation  region  decreased. 
[Ref.  11| 

Ongoing  investigations  have  not  fully  revealed  the  exact  cause  of  vortex 
asymmetry.  The  general  consensus  among  high  angle  of  attack  researchers  is  that  vortex 
asymmetry  is  induced  by  the  location  of  boundary  layer  separation  [Refs.  7,  9j.  Keener 
also  considers  the  possibility  that  vortex  asymmetry  may  be  caused  by  an  inviscid  in¬ 
stability  in  the  pair  of  separated  vortices.  [Ref.  11]  There  are  many  other  three  dimen¬ 
sional  effects  which  must  be  considered  to  complete  a  flowfield  analysis  which  include 
roll,  pitch,  and  yaw  rate;  and  yaw  angle. 

C.  LIFTING  SURFACE  EFFECTS 
1.  Wings. 

Compared  to  aircraft,  missiles  generally  use  low  aspect  ratio  wings.  Mutual  ef¬ 
fects  from  a  wing-body  combination  are  important  since  in  some  cases  the  missile  wing 
span  approaches  body  diameter. 

As  discussed  earlier,  nose  vortices  dictate  flow  behavior  over  the  missile  bodv 
at  high  angles  of  attack.  In  turn,  the  effects  may  also  be  felt  by  the  wings.  Without 
wings,  the  nose  and  body  vortices  tend  to  move  awuy  from  the  body  with  increasing 
angle  of  attack.  The  addition  of  wings  causes  the  vortices  to  take  a  different  trajectory, 
moving  them  closer  to  the  body.  As  the  vortices  move  closer  to  the  body.  the  result  is 
comparable  to  increasing  the  effective  angle  of  attack  causing  unsteady  asymmetric 
vortices.  [Ref.  27] 

A  significant  fraction  of  the  normal  force  from  a  missile  wing  is  proportional  to 
the  square  of  the  angle  of  attack.  This  is  attributed  to  the  vortex  lift  effects  from  flow 
around  a  sharp  leading  edge  delta  wing  4,  As  vortices  curl  around  the  wing  leading  edge, 
a  low  pressure  area  forms  on  top  of  the  wing  providing  lift  as  shown  in  Figure  4.  Since 
a  cruciform  missile  has  an  equilateral  wing  arrangement,  equal  lift  may  be  provided  in 
the  normal  and  side  directions  depending  on  the  wing  angle  of  attack  to  the  flow.  For 


4  A  common  wing  planform  shared  by  numerous  missiles. 
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wings  with  low  aspect  ratio,  a  large  propoition  of  the  lift  generated  h>  the  wing  will  be 
due  to  vortex  lift.  [Ref.  2S] 

The  net  e/lect  of  the  wing-body  combination  seems  to  be  a  reduction  in  the  ef¬ 
fective  angle  of  attack  for  the  onset  of  asymmetric  vortices  and  side  forces. 


Figure  4.  Flow  over  a  delta  wing:  The  left  drawing  shows  the  vortices  curling 
over  the  sharp  leading  edge  lor  low  Reynolds  numbers.  The  right  draw¬ 
ing  is  the  cross  section  B-B.  [Ref.  29]. 

2.  Strakes 

For  an  aircraft,  the  addition  of  strakes  is  to  improve  the  aerodynamic  charac¬ 
teristics  by  enhancing  lift  at  high  angles  of  attack.  Vortices  generated  by  Bow  separation 
from  the  aircraft's  forebode  interacts  with  the  flow  over  the  mam  wing.  The  addition  of 
strakes  will  generate  additional  vortices  improving  the  interaction  between  body  vortices 
and  airflow  over  the  wings.  The  effects  of  wing  strakes  arc  different  on  a  low  aspect  ratio 
wing  missile  since  missile  forebodies  are  comparatively  longer  and  wing  half  span  is 
smaller  when  compared  to  body  radius.  The  flow  becomes  separated  from  the  body  at 
a  moderate  angle  of  attack  before  passing  over  the  strakes  [Ref.  >].  The  strakes  on  a 


11 


missile  may  lessen  the  effect  of  asymmetric  body  vortices  since  the  strakes  would  cause 
interference  with  the  crossflow  component  around  the  body.  "1  he  real  benefit  of  incor¬ 
porating  strakes  with  low  aspect  ratio  wings  is  from  the  strong  vortices  produced  by  the 
strake  which  enhance  vortex  lift  effects. 

3.  Tails 

The  tails  of  a  missile  have  little  influence  on  flowfields  on  the  forebody,  espe¬ 
cially  at  low  angles  of  attack.  At  high  angles  of  attack  interference  on  the  flowfield 
caused  by  wings  and  the  afterbody  may  effect  tailflow  which  depends  a  great  deal  on  the 
placement  of  the  wings,  the  length  of  the  afterbody  between  the  wings  and  tails,  and  the 
missile's  angle  of  attack.  For  higher  angles  of  attack,  nose  and  wing  vortices  may  be 
sufficiently  displaced  away  from  the  body  to  have  only  a  slight  effect  on  the  tailflow. 
Should  a  large  distance  between  wings  and  tails  exist,  afterbody  vortices  may  dominate 
the  airflow  over  the  tail  surfaces.  Hence,  most  missile  configurations  with  long 
afterbodies  utilize  canard  or  wing  control  while  missiles  employing  tail  control  have 
short  afterbodies  and  low  aspect  ratio  wings  with  long  root  chords.  The  additional 
asymmetric  vortices  from  the  afterbody  will  also  contribute  to  the  overall  side  force  on 
the  missile.  [Ref.  28:  p.  172] 

D.  TURBULENCE 

1.  Intensity  and  Length  Scale 

Turbulence  denotes  an  irregular  fluctuation  or  disturbance  of  a  fluid  changing 
motion  when  superimposed  on  a  flowfield,  In  terms  of  scale  (size  and  energy),  the 
flowfield  may  be  magnitudes  higher  than  turbulence.  The  disturbance  exists  when  lluid 
particles  overcome  the  constraining  effects  of  viscosity  and  are  able  to  fluctuate  in  the 
flowfield.  Turbulence  exists  for  large  Reynolds  numbers  where  longer  llowficld  length 
scales  are  present.  A  horizontal  wind  or  crosswind  may  be  considered  as  a  large  scale 
flowfield  and  turbulence  becomes  a  smaller  component  of  the  crosswind.  Further,  when 
a  flowfield  contains  turbulence,  it  can  be  considered  a  turbulent  flowfield  and  so  a 
crosswind  can  be  though  of  as  a  turbulent  crosswind.  It  should  be  understood  that  any 
wind  for  practical  purposes  is  turbulent,  but  since  crosswind  size  and  energy  scales  are 
magnitudes  larger  than  the  turbulence  scales  contained  in  the  crosswind,  we  can  treat  the 
crosswind  as  homogeneous  and  two  dimensional. 

As  a  large  scale  homogeneous  flowfield,  the  crosswind  can  be  viewed  as  exerting 
a  uniform  flow  over  the  missile  body.  Flow  separation  mechanisms  discussed  earlier  for 
a  two  dimensional  cylinder  in  a  steady  flowfield  apply  in  the  crosswind.  Understanding 
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the  effects  of  small  scale  turbulence  on  the  flow  separation  mechanisms  is  important 
since  vortex  generation  and  asymmetry  determine  side  force  direction  and  magnitude. 
Small  scale  turbulence  effects  are  dependent  on  turbulence  length  scale  and  intensity 
which  will  be  discussed  below.  Turbulence  intensity  and  length  scale  components  exists 
in  the  streamwise  and  spanwise  directions  relative  to  the  flowfield.  For  purposes  of  this 
and  following  discussions  only  the  streamwise  component  in  the  flowfield  will  be  dis¬ 
cussed. 

Turbulence  intensity  is  the  measure  of  the  relative  magnitude  of  velocity  fluc¬ 
tuations  in  the  flowfield.  Mathematically,  for  a  horizontal  flowfield  or  crosswind,  the 
turbulence  intensity  Tu  is  the  ratio  of  the  root-mean-square  (;ms)  velocity  fluctuation 
u',  to  the  mean  velocity  component  in  the  flowfield.  UM. 


Since  the  turbulence  intensity  measures  the  magnitude  of  velocity  fluctuation,  a  higher 
intensity  indicates  higher  kinetic  energy  and  correlates  to  more  turbulent  flow. 

Turbulence  length  scale  Lu  describes  the  time  average  measure  of  the  fluid  dis¬ 
turbance  eddies.  The  measure  is  time  averaged  since  these  disturbances  constantly 
change  size.  The  size  of  the  disturbances  are  determined  by  external  conditions  associ¬ 
ated  with  the  How.  An  example  of  large  scale  turbulence  are  flow  disturbances  caused 
by  flovvfields  over  large  buildings  in  a  city,  while  small  scale  turbulence  is  generated  bv 
flow  through  mesh  screens.  [Ref.  31:  p.  557] 

Turbulence  length  scale  to  body  size  ratio  may  dictate  the  manner  in  which 
turbulence  affects  the  VLSAM  flowfield.  The  turbulence  length  scale  can  be  compared 
to  missile  length,  LyL,  .  or  missile  diameter,  Lu:L,j  .  When  the  turbulence  length  scale 
is  many  times  larger  than  the  missile  length,  LU$>L„  the  effect  on  the  missile  is  somewhat 
like  a  uniform  time  independent  steady-state  flowfield.  especially  at  low  turbulence  in¬ 
tensities.  The  effects  of  the  flowfield  on  vortex  development  are  largely  dictated  by  those 
conditions  and  factors  discussed  earlier  for  a  two  dimensional  cylinder.  When  the  tur¬ 
bulence  length  scale  is  comparable  to  the  body  length,  LU*L,  ,  the  flowfield  is  distinctly 
non-uniform  and  may  cause  unwanted  rolling,  pitching  and  yawing  motion  [Ref.  32J. 
In  both  cases,  the  missile  body  motion  results  from  the  actual  flowfield  conditions  about 
the  missile  body. 
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Perhaps  the  most  significant  turbulence  condition  is  when  the  turbulence  length 
scale  is  smaller  than  the  body.  LU<^L,  ,  and  most  important,  when  it  is  comparable  or 
smaller  than  missile  diameter,  Lu<Ld .  The  small  scale  turbulence  present  in  the  fiowfield 
may  change  the  boundary  layers  present  on  the  missile  surface,  particularly  when  the 
turbulence  length  scale  is  the  same  size  or  smaller  than  the  boundary  layer  thickness. 
An  important  parameter  for  small  scale  turbulence  is  the  turbulence  intensity,  since  the 
energy  in  the  turbulence  may  trigger  the  boundary  layer.  In  contrast  to  larger  turbu¬ 
lence  scales  where  the  primary  influence  is  on  the  fiowfield,  the  effects  of  small  scale 
turbulence  is  felt  on  the  boundary  layer  and  flow  separation  on  the  surface  of  the  missile. 
[Ref.  33[ 

Both  large  and  small  scale  turbulence  length  scales  actually  coexist  in  the 
fiowfield  since  the  turbulence  length  scale  is  the  mean  of  the  disturbances  in  the 
fiowfield.  The  combination  of  different  length  scales  may  add  complexity  to  the  qual¬ 
itative  analysis  of  the  fiowfield.  The  cascade  effect  explains  how  turbulence  in  varying 
length  scales  exists  in  the  fiowfield.  As  the  disturbances  in  the  fiowfield  experience 
strain,  the  disturbances  will  break  up  into  finer  disturbances  of  smaller  scale  and  de¬ 
creasing  energy.  The  process  repeats  itself  or  cascades  until  the  small  disturbances 
eventually  disappear  owing  to  viscosity.  When  the  disturbances  decrease  in  size  and  en¬ 
ergy,  their  individual  intensities  will  decrease  at  a  faster  rate.  Thus,  the  turbulence  length 
scale  average  will  be  biased  towards  the  larger  turbulence  scales  in  the  fiowfield  and  gives 
an  accurate  quantitative  representation  of  all  the  length  scales  within  the  fiowfield  since 
the  energy  and  intensity  of  the  larger  scale  turbulence  predominate.  [Ref.  34] 

2.  Boundary  Layer  Effects 

Turbulence  length  scales  on  the  order  of  boundary  layer  thickness  may  change 
the  boundary  layer  on  the  missile  surface.  This  small  scale  turbulence  will  be  referred  to 
as  boundary  layer  scale  turbulence,  and  aithough  three  dimensional  and  inhomogeneous, 
can  be  visualized  as  two  dimensional  crossflow  components  over  a  cylinder.  When  the 
boundary  layer  scale  turbulence  in  the  crossflow  reaches  the  missile  surface  boundary 
layer,  it  contributes  energy  as  it  dissipates  into  the  boundary  layer.5  This  turbulence 
interaction  results  in  an  earlier  laminar  transition  and  delays  separation  of  the  turbulent 
boundary  layer.  This  causes  flow  separation  to  occur  below  the  critical  Reynolds  num¬ 
ber  for  laminar  flow,  postponing  asymmetric  vortex  shedding.  Reynolds  number  effects 
dominate  at  low  turbulence  intensities  for  smooth  surfaces.  When  surface  roughness  in- 

5  In  accordance  with  the  laws  of  momentum  and  energy  conservation. 


creases,  the  boundary  layer  separation  becomes  less  dependent  on  Reynolds  number  and 
more  dependent  on  turbulence  intensity.  The  net  effect  postpones  the  side  force  induced 
by  the  asymmetric  vortices  and  allows  the  missile  to  operate  at  higher  angles  of  attack 
without  induced  side  forces.  [Refs.  33,  35] 

The  application  of  boundary  layer  scale  turbulence  to  add  energy  to  the 
boundary  layer  is  a  well  known  method  of  inducing  earlier  flow  transition  and  delayed 
flow  separation  to  eliminate  or  reduce  the  effects  of  asymmetric  vortex  induced  side 
forces.  Among  the  numerous  methods  implemented  to  add  turbulence  to  the  boundary 
layer  are  placing  roughness  strips  in  various  patterns  at  noses  of  the  missile  models  uti¬ 
lized  in  wind  tunnel  experiments.  [Refs.  11.26j 
3.  Vortex  effects 

When  the  vortices  fiom  the  missile  body  are  formed,  they  become  detached  from 
the  body  at  high  angles  of  attack.  The  interaction  of  turbulence  on  the  order  of  the 
vortex  s’zc  may  give  additional  insight  to  the  disposition  of  asymmetric  vortices.  This 
small  scale  turbulence  will  be  referred  to  as  vortex  scale  turbulence  and  is  magnitudes 
larger  than  boundary  scale  turbulence.  Vortex  scale  turbulence  may  be  compared  to 
body  diameter  size  Lu^Ld  .  Objects  in  the  flowfield.  much  larger  then  roughness  strips, 
may  be  used  to  generate  vortex  scale  turbulence.  Vortex  scale  turbulence  may  be  gen¬ 
erated  by  strakes.  support  hardpoints.  or  inlets.  The  effects  of  vortex  scale  turbulence 
may  influence  the  position  of  the  asymmetric  vortices  and  may  contribute  dissipated 
turbulence  energies  as  with  boundary  layer  turbulence.  Since  the  relative  positioning 
and  magnitude  of  the  vortex  structure  governs  the  pressure  distribution  on  the  surface 
of  the  missile,  the  vortex  scale  turbulence  could  dictate  the  position  and  magnitude  of 
the  asymmetric  induced  side  forces.  Research  into  this  phenomenon  has  been  con¬ 
ducted  as  another  means  of  alleviating  asymmetric  vortex  induced  side  forces.  Methods 
include  adding  nose  strakes,  nose  helical  vortex  generators,  and  nose  booms  to  missile 
models.  [Refs.  11,  18,  36] 

E.  VS  LAM  LAUNCH  ENVIRONMENT 
1.  The  Marine  Environment 

Conditions  in  the  marine  atmospheric  environment  may  have  a  significant  effect 
on  the  VLSAM.  In  the  atmosphere,  turbulence  conditions  exist  almost  everywhere,  in 
particular,  within  the  atmospheric  boundary  layer  (ABL)  which  is  the  lowest  part  of  the 
atmosphere  [Ref.  37:  p.  179].  The  conditions  in  the  ABL  are  the  result  of  atmospheric 
flow  over  the  sea  surface.  The  flowfield  cannot  be  strictly  classified  as  a  boundary  layer 
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j  flow  since  non-uniform  boundary  conditions  under  constant  fluctuation  and  changes 

i  almost  always  exist  [Ref.  32]. 

j  The  lowest  ten  percent  of  the  ABL  is  referred  to  as  the  surface  layer  and  is 

typically  50  meters  in  height  but  may  vary  from  5-200  meters.  The  surface  layer  is 
characterized  by  distinct  variations  in  wind  speed,  nearly  vertical  heat  and  mass  fluxes, 
and  other  meteorological  fluctuations  with  height  [Ref.  3S:  p.  236].  Fluctuations  in  wind 
variations,  especially  at  sea  may  reach  fifty  percent  of  the  mean  wind  speed 
[Ref.  31:  p.  557].  The  surface  layer  is  also  characterized  by  intense  small  scale  turbu¬ 
lence  which  is  mechanically  produced  by  surface  roughness  or  friction  from  waves  on  the 
ocean  surface. 6  The  lowest  portion  of  the  surface  layer  is  termed  the  roughness  layer  and 
is  where  the  surface  has  the  greatest  effect  on  generating  turbulence  and  influencing  fluid 
motion.  Mean  flow  in  the  roughness  layer  in  largely  three  dimensional  and  in¬ 
homogeneous  in  nature.  Small  scale  turbulence  dominates  flow  near  the  surface.  As  al¬ 
titude  increases,  small  scale  turbulence  decreases  rapidly  as  large  scale  turbulence  formed 
by  convective  conditions  in  the  atmosphere  prevail.  Several  roughness  heights  above  the 
surface,  the  mean  flowfield  becomes  nearly  horizontally  homogeneous  and  two  dimen¬ 
sional.  Over  sea.  the  mean  surface  layer  and  especially  the  roughness  layer  are  more 
inhomogeneous  and  three  dimensional  than  over  land. 

2.  Turbulence 

The  roughness  length  z0  is  a  measure  of  surface  roughness  as  a  function  of  the 
mean  wind  velocity  and  height  above  the  surface.  Surface  layer  turbulence  length  scale 
and  intensity  can  be  empirically  determined  from  the  aforementioned  parameters.  For 
typical  open  ocean  roughness  lengths  in  the  range  10  3<z,:,  <  10~2  meters,  turbulence 
intensities  arc  on  the  order  of  13  to  17  percent  at  a  10  meter  elevation  and  mean  wind 
speed  of  approximately  25  m  sec  [Ref.  39].  These  high  turbulence  intensities  are  for  high 
crosswind  velocities  and  rough  seas;  extreme  conditions  which  may  be  encountered. 
Tieleman  states  an  average  turbulence  intensity  of  6  to  12  percent  to  exist  below  a  20 
meter  height  above  the  ocean  surface  [Ref.  32], 

From  previous  discussions,  the  turbulence  length  scale  to  body  scale  ratio  is  an 
important  parameter  for  flowfield  behavior.  For  the  conditions  stated  in  the  preceding 
paragraph,  the  longitudinal  turbulence  length  scales  Lu  are  in  the  range  SO  <  Lu  <  90 
meters  (or  262  <  Lg  <  295  feet).  For  a  typical  missile  with  a  1.1  foot  diameter,  the  tur¬ 
bulence  to  mjssile  diameter  length  scale  ratio  is  about  Lu:Ld  s:  280:1  and  would  have 


6  This  small  scale  turbulence  is  larger  than  the  missile  length. 


little  effect  on  boundary  layer  development.  Vet  the  possibility  exists  where  boundary 
layer  and  vortex  scale  turbulence  may  be  present  due  to  the  cascade  effect  from  large 
scale  turbulence  and  from  crosswind  interaction  with  the  ship's  superstructure,  lending 
themselves  to  vortex  and  boundary  layer  interactions. 

3.  Crossttind 

A  typical  VLSAM  launching  with  a  lOg  acceleration  reaches  an  altitude  of  ap¬ 
proximately  56  feet  travelling  at  a  vertical  velocity  of  Vv  3:164  ft  sec.  The  VLSAM  is  still 
well  within  the  surface  layer  environment  (50  m)  and  subject  to  crosswind  and  turbu¬ 
lence  effects.  The  preceding  paragraph  indicates  turbulence  length  scales  are  relatively 
long  in  the  surface  layer.  In  such  a  case,  the  large  scale  flowfield  or  crosswind  may  dic¬ 
tate  aerodynamic  effects  on  the  VLSAM.  As  an  example,  if  a  ship  travels  at  20  knots 
and  a  mean  wind  speed  of  20  m  sec  are  combined  with  the  VLSAM  launch  velocity,  the 
resultant  vector  is  at  191  ft  sec  at  an  angle  of  31°  from  the  missile  direction.  The  result 
is  that  the  VLSAM  flies  at  an  effective  angle  of  attack  of  31°  at  191  ft  sec.  placing  it  in 
Regime  III,  the  asymmetric  vortex  region,  as  it  travels  through  the  surface  layer. 
[Ref.  3] 

A  missile  at  a  lower  altitude,  closer  to  the  ship,  will  have  a  lower  velocity  and 
therefore  it  will  be  flying  at  an  even  greater  effective  angle  of  attack.  Gregoriou  states 
that  a  VLSAM  that  pitches  over  towards  a  target  may  experience  an  effective  angle  of 
attack  up  to  or 3:50°  (Ref.  2].  These  examples  should  illustrate  the  distinct  possibility  of 
asymmetric  vortex  induced  side  forces  on  a  VLSAM  during  launch  phase. 

4.  Ship  Airwake. 

Changing  flowfields  and  turbulence  may  be  dictated  by  the  ship's  hull  and 
superstructure.  This  flowfield  or  ship  airwake  may  increase  the  crosswind  velocity,  cause 
significant  crosswind  gradients,  and  non-uniformities  in  the  flowfield  increasing  turbu¬ 
lence  intensities  while  decreasing  turbulence  length  scales.  These  effects  occur  in  and 
around  the  vicinity  of  the  ship  which  acts  as  a  turbulence  generator.  Results  of  earlier 
discussions  produced  the  following  parameters  in  the  marine  surface  layer; 

Mean  crosswind  velocity  -  2U  <  <30  (m  sec) 

Turbulence  intensity  -  1 3  <  Tu  <17  (percent) 

Turbulence  length  scale  -  SO  <  Lu  <  90  (meters) 
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These  values  may  change  due  to  ship  airwake  effects  in  favor  of  higher  turbulence  in¬ 
tensities  and  smaller  length  scales,  conditions  prone  to  affect  the  missile  boundary  layer 
causing  delayed  separation.  (Ref.  39]  j 

5.  Additional  launch  considerations 

Other  considerations  which  may  influence  the  YLSAM  at  launch  are  dependent 
on  the  launch  platform,  the  mechanics  of  a  canister  launch,  and  the  aerodynamic  effects 
of  the  the  missile  controlling  itself.  These  conditions  are  essential  to  the  overall  aero¬ 
dynamic  effects  on  the  YLSAM  as  it  exits  the  canister. 

The  launch  platform  undergoes  constant  roll,  pitch  and  yaw  motion  and  trans¬ 
mits  this  motion  to  the  YLSAM  as  it  travels  through  the  the  guide  mechanism  in  the 
canister.  Increased  severity  of  the  ship's  motion  adds  complexity  to  the  initial  velocity 
vector  as  the  YLSAM  leaves  the  guide. 

The  launch  feature  unique  to  a  YLSAM  is  that  the  canister  used  lor  storage  is 
also  the  launcher.  Unlike  a  rail  launched  missile  where  motor  exhaust  dissipates  into  the 
surroundings,  a  canister  or  tube  constrains  a  portion  of  the  rocket  motor  gases  resulting 
in  a  muzzle  blast  effect.  As  the  YLSAM  rocket  motor  ignites,  most  of  the  exhaust  gases 
are  typically  vented  via  a  plenum  and  uptake  arrangement.  Still,  there  is  a  significant 
amount  of  backflow  which  Hows  around  the  YLSAM  as  it  leaves  the  canister.  The 
rocket  motor  exhaust  from  the  uptake  may  also  interact  with  the  flowficlds  about  the 
missile.  The  environment  of  a  rocket  motor  ignition  is  severe  in  itself,  with  high  mag¬ 
nitudes  of  heat,  shocks,  and  vibration.  (Ref.  40] 

The  YLSAM  itself  may  i  npose  some  form  of  flight  control,  aerodynamic  or 
thrust  vector,  to  maintain  flight  attitude  during  the  launch  phase.  Should  the  missile 
change  flight  attitude,  the  flowfield  around  it  will  be  altered. 

These  considerations  are  introduced  to  bring  out  the  numerous  factors  which 
may  influence  the  YLSAM  flowfield  and  cause  in-plane  and  out-of-plane  forces  and 
moments  during  the  launch  phase.  They  will  add  complexity  to  a  full}  comprehensive 
flowfield  analysis  of  the  YLSAM  which  is  beyond  the  scope  of  the  present  investigation. 
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II.  EXPERIMENT  AND  PROCEDURES 

A.  APPARATUS 

The  experimental  investigation  was  conducted  at  the  Naval  Postgraduate  School 
wind  tunnel  test  facility.  The  VLSAN1  model  was  mounted  on  a  six-component  strain 
gage  balance  and  sting  support  assembly.  Turbulence  generating  grids  and  the  test  fa¬ 
cility  steady-state  data  acquisition  hardware  software  complete  the  experimental  appa¬ 
ratus  discussed  in  this  section. 

I.  Wind  Tunnel 

The  horizontal-flow  low-speed  tunnel  used  for  this  experiment  is  located  in  the 
Naval  Postgraduate  School  wind  tunnel  test  facility  and  shown  in  Figure  .5.  The  wind 
tunnel  is  powered  by  a  one-hundred  horsepower  electric  motor  driving  a  three  blade 
variable-pitch  propeller  via  a  four-speed  Dodge  truck  transmission.  Immediately  fol¬ 
lowing  the  propeller  blades  are  a  set  of  stator  blades  to  help  straighten  the  flow.  Turning 
vanes  installed  at  each  comer,  and  two  fine  wire  mesh  screens  placed  at  the  entrance  to 
the  settling  chamber  smoothen  the  airflow  and  reduce  turbulence.  The  settling  chamber 
to  test  section  contraction  ratio  is  approximately  10:1.  A  heavy  wire  mesh  screen  is  in¬ 
stalled  at  the  end  of  the  test  section  diffuser  to  prevent  foreign  object  damage  to  the  fan 
section.  [Ref.  41:  pp.3-1  to  3-2] 

The  test  section  operates  at  atmospheric  pressure  and  measures  45  inches  by  32 
inches.  Corner  fillets  house  the  test  section  lighting  and  help  reduce  synergetic  boundary 
laser  effects  at  the  wall  interactions.?  To  counteract  a  reduction  in  frccstreum  pressure 
due  to  boundary  layer  growth,  the  walls  of  the  test  section  are  slightly  diverged.  Im¬ 
mediately  downstream  of  the  test  section  are  breather  slots  around  the  perimeter  of  the 
wind  tunnel.  These  slots  allow  air  to  enter  and  leave  the  tunnel  unhindered,  replenishing 
air  lost  through  leaks  and  ensuring  a  uniform  test  section  pressure.  A  reflection  plane 
located  in  the  center  of  the  test  section  provides  a  horizontal  surface  to  the  flowfield  and 
reduces  the  test  section  available  height  to  2S  inches.  A  flush  turntable  at  the  center  al¬ 
lows  model  the  operator  to  change  model  pitch  or  angle  of  attack  during  wind  tunnel 
operation.! Ref.  41:  pp.  3-3  to  3-4] 


"  The  fillets  reduce  the  area  of  the  test  section  from  10  ft.2  to  9.88  ft.2  Fillets  are  employed 
throughout  the  wind  tunnel  at  wall  interactions. 
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Figure  5.  Naval  Postgraduate  School  Low  Speed  Wind  Tunnel:  [Ref.  41:  p.  3-7] 

Flow  measurement  components  utilized  with  the  wind  tunnel  are  a  dial 
thermometer,  manometer,  and  pitot  static  tube.  Temperature  in  the  wind  tunnel  is 
measured  by  a  dial  thermometer  extending  into  the  settling  chamber.  Dynamic  pressure 
in  the  wind  tunnel  is  measured  by  the  static  pressure  difference  between  the  test  section 
and  the  settling  chamber  using  a  water  filled  manometer.  Static  pressure  in  the  settling 
chamber  is  measured  by  four  static  pressure  taps,  one  on  each  wall.  The  test  section  has 
another  four  static  pressure  taps,  one  on  each  wall,  upstream  of  the  test  section  to  pre¬ 
clude  interference  from  the  model.  The  pressure  taps  at  each  section  arc  connected  via 
a  common  manifold  prior  to  feeding  into  the  manometer.  The  manometer  measures 
pressure  differences  in  centimeters  of  water,  which  is  the  test  section  dynamic  pressure, 
and  is  converted  to  the  actual  wind  tunnel  velocity  fcv  equation  [5]. 
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Where: 


L'm  =  measured  velocity  (-^r) 

2.046  =  conversion  factor 

cm  H20  =  manometer  reading  in  cm  of  HaO 

0.93  =  calibration  factor. 8 
lbm 

p  =  air  density  (— rp) 


A  pressure  transducer  circuit,  connected  to  the  manometer  displays  test  section  dynamic 
pressure  on  a  digital  readout.  Once  calibrated,  the  digital  manometer  has  a  greater  ac¬ 
curacy  than  the  manometer  and  was  used  to  measure  the  dynamic  pressure.  The  pitot 
static  tube  is  read  from  a  aircraft  airspeed  indicator  and  is  used  only  for  test  section  ve¬ 
locity  approximations  when  starting  and  stopping  the  wind  tunnel. 

2.  VS  LAM  Model 

The  VLSAM  model  was  designed  to  be  representative  of  a  cruciform  tail-control 
missile  with  vers-  low  aspect  ratio  wings  or  dorsal  fins.  The  model  was  fabricated  from 
6061  and  2024  aluminum  alloy  by  Naval  Postgraduate  School  personnel.  The  model  ac¬ 
commodates  the  one-inch  precision  balance  and  sleeve  used  in  the  experiment  and  has 
internal  reserve  space  for  fitting  pressure  measurement  instrumentation.  Four  sections 
form  the  major  parts  of  the  model.  The  body  section  is  a  hollow  cylinder  with  locating 
pin  attachment  points  for  the  balance,  sleeve,  wings  and  tails.  The  body  mounts  on  the 
balance,  or  by  an  adapter  for  use  without  the  balance  to  mount  the  model  on  the  sup¬ 
port  or  sting.  Body  roll  angle,  with  or  without  the  balance,  can  be  varied  in  90&  incre¬ 
ments.  The  nose  piece  attaches  to  the  body  forward  of  the  balance.  Nose  roll  angle  may 
be  varied  in  45c  increments  independent  of  the  body.  The  nose  is  a  pointed  ogive  shape 
with  a  short  constant  diameter  section  aft  of  the  ogive  which  provides  an  interference 
free  interface  between  the  ogive  shape  and  the  body  section.  Four  very  low  aspect  ratio 
wings  with  strakes  comprise  the  cruciform  wing  section  and  are  mounted  equilaterally 
in  fixed  position  along  the  model  axis.  Likewise,  in  the  tail  section,  are  four  tail  control 
fins  which  are  fixed  in  postion  at  0°  incidence  relative  to  the  model  axis  mounted  aft  of 
each  wing.  All  parts  are  rigidly  connected  to  the  model  body  by  countersunk  machine 


8  Total  pressure  correction  for  the  settling  chamber. 
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screws.  The  surface  of  the  VLSAM  model  is  free  from  proturberunces  and  is  buffed  to 
a  polished  finish.  Detailed  dimensions  of  the  VLSAM  model  are  shown  in  Figure  6. 
A  short  summary  of  ke\  dimensions  are  listed  below; 
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Total  length 
Base  diameter 
Length  diameter  ratio 
Ogive  nose  length 
Ogive  diameter  ratio 
W  ing  span  root  chord 
Tail  span  root  chord 
Center  of  pressure 


22.85  in. 

1.75  in. 

13.06 
4.0  in. 

2.29 

3.13  in. ,13.55  in. 

5.50  in. T.70  in. 

13.5  inches  aft  of  nose  tip  (approx.) 


NAVAL  POSTGRADUATE  SCHOOL 
SURFACE-TO-AIR 
MISSILE  MODEL 


3.  Balance 


Force  and  moment  measurements  were  taken  by  a  six-component,  one-inch. 
|  strain-gage  precision  balance.  The  balance  was  on  loan  from  the  NASA-Ames  Research 

Center  under  the  Navy- NASA  Joint  Institute  of  Aeronautics.  Balance  calibration  was 
performed  by  NASA-Ames.  Calibration  data  and  data  conversion  values  were  also 
provided  and  are  listed  in  Appendix  A  [Ref.  42).  A  Wheatstone  bridge  circuit  capable 
I  of  measuring  positive  and  negative  strain  comprises  each  channel  and  is  oriented  to 

measure  strain  in  a  chosen  direction.  The  directions  relative  to  the  balance  axis  are 
normal,  side,  and  axial  with  one  gage  measuring  rolling  moment.  Channel  capacities 
and  accuracies  are  listed  in  Table  1.  Individual  channel  accuracies  are  based  on  full 
I  scale  (maximum)  load. 


Table  1.  BALANCE  CHANNELS:  [Ref.  42  [. 


Channel 

Component 

Maximum  Load 

°0  Accuracy 

NT 

Normal  Force 

400  lbs. 

0.056 

N2 

Normal  Force 

400  lbs. 

0.049 

SI 

Side  Force 

200  lbs. 

0.115 

S2 

Side  Force 

200  lbs. 

0.132 

A 

Axial  Force 

100  lbs. 

O.I53 

R 

Roll  Moment 

21  f't-lbs. 

0.204 

A  machined  sleeve  provides  a  close  tolerance  fit  between  the  balance  and  the 
interior  of  the  VLSAM  model.  The  model  is  seated  to  the  balance  by  locating  pins  such 
that  the  balance  center  coincides  with  the  approximate  center  of  pressure  of  the  model. 
Since  the  balance  supports  the  model,  it  provides  the  mechanical  interface  between  the 
model  and  sting  support. 

4.  Model/Balance  Support 

The  VLSAM  model  and  balance  support,  also  called  the  sting  mount,  is  rigidly 
fixed  in  the  wind  tunnel  test  section  by  the  reflection  plane  turntable  at  the  base  and  an 
aluminum  reinforced  clear  plexiglas  section  at  the  top.  The  sting  mount  changes  model 
pitch  in  the  horizontal  direction  with  the  pivot  point  located  to  coincide  with  the  ap- 
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proximate  centers  of  the  VLSAM  model  and  balance.  Heavy  duty  construction  in 
stainless  steel  and  aluminum  keeps  deflections  and  vibrations  to  a  minimum.  The  pre¬ 
sented  cross-sectional  area  of  the  sting  mount  to  the  test  section  How  is  minimal  for  all 
pitch  angles  to  reduce  blockage  effects.  The  reflection  plane  turntable  rests  on  a 
heavy-weight  pedestal  to  prevent  deflections  and  vibration.  The  entire  assembly  rotates 
via  a  chain  gear  drive  powered  by  an  electric  motor.  Pitch  and  pitch  rate  are  set  by  the 
operator  on  the  motor  controller,  or  by  a  relay  actuator  controlled  by  the  data  acquisi¬ 
tion  system.  A  ring  scale  on  the  pedestal  indicates  model  pitch  in  1/10  degree  increments. 


Figure  7.  Photograph  of  VLSAM  model  in  the  test  section. 


Figure  8.  Planview  of  VLSAM  model  in  the  test  section.  •> 

5.  Turbulence  Grids 

The  turbulence  grids  were  designed  for  use  in  the  low  speed  wind  tunnel  to 
generate  turbulence  in  varying  intensities  and  length  scales.  Each  of  the  lour  grius  or  • 

screens  is  mounted  in  a  wooden  frame  and  placed  incites  forward  of  the  mode!  sting 
mount  pivot  point.  The  specifications  for  the  four  grids  arc  Listed  in  Table  2.  Three  of 
the  four  grids  are  square-mesh  square-bar  ami  fabricated  from  wood.  They  employ  bi-  j 

plane  construction  which  generates  a  nearly  isotropic  homogeneous  turbulence.  The 
fourth  grid  is  also  square  mesh,  but  was  fabricated  from  round  wire.  Two  wire  screens  | 

are  attached  together  for  a  fine  mesh.  Turbulence  intensities  and  length  scale  for  each  j 

of  the  grids  arc  shown  in  Figure  9  and  Figure  10.  These  grids  were  constructed  by 
Roane  who  also  measured  the  turbulence  intensities  and  length  scales.  |> 
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Table  2.  GRID  SPECIFICATIONS:  (Ref. 


?  :  P-2S1 


Grid 

Mesh  Width  (in.) 

Bar  Diameter  (in.) 

Mesh/Diameter 

Material 

One 

5.  oo 

1 .00 

X 

Wood 

Tsvo 

3.75 

0.75 

5 

Wood 

Three 

2.50 

0.50 

5 

Wood 

Four 

1.00 

0.0625 

16 

Wire 

DISTANCE  FROM  GRID  (INCHES) 


Figure  9.  Grid  Turbulence  Intensity:  The  dashed  line  indicates  the  model  pivot 
axis.  (Ref.  3:  pp.  45] 
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Figure  10.  Grid  Turbulence  Length  Scales:  The  dashed  line  indicates  the  model 
pivot  axis.  [Ref.  3:  pp.  48] 


Table  3.  GRID  TURBULENCE  PARAMETERS:  At  the  model  phot  axis. 
_ [Ref.  3:  pp.  44-49) _ 


Grid 

Intensity 

(percent) 

Length  Scale 
(in.) 

T  urbulence/ 
Model  Dia. 

Dynamic  Pressure 
(lb/ft-'l 

One 

3.31 

1.84 

1.05 

15.35 

Two 

2.7S 

1.56 

0.89 

14.88 

Three 

1.88 

1.08 

0.62 

I6.3S 

Four 

0.47 

0.27 

0.15 

15.61 

None 

0.23 

- 

- 

15.  S5 

The  dynamic  pressures  shown  in  Table  3  are  the  true  dynamic  pressure  in  the  test  sec¬ 
tion  from  Roane  [Ref.  3:  pp.  44-49].  The  grid  generated  turbulence  length  scales  to 
missile  length  ratios  at  the  test  section  are  smaller  than  1:12  which  suggests  that  the 
turbulence  effects  on  the  flowfield  and  vortices  may  be  due  to  boundary  layer  scale 
and  or  vortex  scale  turbulence.  The  turbulence  length  scales  to  body  diameter  ratios  are 
on  the  order  of  and  smaller  than  one  as  seen  in  Table  3.  Conclusions  on  the  effects  of 
the  grid  generated  turbulence  with  regards  to  changing  length  scales  at  constant  intensity 
or  changing  intensities  with  constant  length  scale  are  not  possible  with  the  present  grid 
turbulence  parameters. 


Figure  1 1.  Turbulence  Grids  2,  3,  and  4.:  Clockwise  from  lower  left  hand  comer; 
Grid  4,  Grid  3,  Grid  2. 
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Figure  12.  Turbulence  Grid  1  installed  in  the  wind  tunnel. 

6.  Data  Acquisition  Hardware 

The  data  acquisition  hardware  consists  of  test  facility  components  to  operate 
the  precision  balance  in  the  wind  tunnel.  Each  strain-gage  in  the  balance  is  supported 
by  a  signal  conditioning  circuit  which  supplies  excitation  voltage  and  allows  zeroing  and 
calibration  of  the  balance  circuit.  Each  balance  channel  from  the  signal  conditioner  unit 
is  fed  into  a  Hewlett-Packard  relay  multiplexer  which  sequences  through  each  channel 
as  required  by  the  data  acquisition  software.  Prior  to  being  measured,  the  balance  signal 
is  fed  through  a  1000  gain  low-noise  amplifier  to  improved  resolution  for  small  magni¬ 
tude  signals.  A  Hewlett-Packard  digital  multimeter  converts  the  analog  voltage  signal 
into  a  digital  signal  for  use  by  the  data  acquisition  program.  An  1BM-AT  microcom¬ 
puter  executes  the  data  acquisition  software  and  stores  the  collected  data.  Other  items 
utilized  in  the  data  acquisition  process  are  the  wind  tunnel  measurement  apparatus  dis¬ 
cussed  earlier.  Data  acquisition  hardware  components  are  shown  in  Figure  13. 
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Figure  13.  Data  Acquisition  Hardware 
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7.  Data  Acquisition  Software 

The  data  acquisition  software  consists  of  a  shell  program  to  control  the 
Hewlett-Packard  instruments  and  a  data  conversion  program  to  convert  the  raw  balance 
readings  into  a  usable  form  and  store  them.  All  readings  are  time  averaged  and  represent 
steady  state  measurements.  The  initial  balance  reading  in  volts  is  manipulated  with  the 
balance  calibration  constants  to  give  force  values  in  pounds  and  moment  values  in 
foot-pounds.  Another  program  translates  the  force  and  moment  values  into  a  coefficient 
form  accounting  for  test  conditions  and  corrections  in  the  wind  tunnel.  The  data  ac¬ 
quisition  programs  are  listed  in  Appendix  C  and  Appendix  D. 

B.  EXPERIMENTAL  CONDITIONS 

The  introduction  discussed  many  variables  which  could  affect  flow  separation  and 
vortex  structure  at  high  angles  of  attack.  To  simplify  the  various  outcomes  of  the  test 
due  to  numerous  perturbations,  the  following  variables  were  kept  as  constant  as  exper¬ 
imental  conditions  would  allow: 

1.  Reynolds  number 

2.  Test  section  dynamic  pressure. 

3.  Test  section  velocitv. 

4.  Nose  geometry  and  roll  angle. 

The  experimental  reference  pressures  used  in  the  turbulence  mapping  by  Roane 
(Ref.  3]  were  duplicated  in  this  experiment  in  order  to  use  the  turbulence  grid  length 
scales  and  intensities  shown  in  Table  3,  Reference  dynamic  pressures  for  the  grid  and 
no  grid  conditions  are  listed  below.  These  reference  pressures  ensure  a  similar  test  section 
velocity  since  the  total  pressure  difference  measured  by  manometer  is  in  error  due  to  the 
disturbances  caused  by  the  screens.  Since  actual  dynamic  pressure  in  the  test  section  is 
determined  by  the  test  section  velocity,  the  velocities  in  the  test  section  from  Roane's 
turbulence  mapping  were  used  to  calculate  the  actual  dynamic  pressures.  These  refer¬ 
ence  pressures  are: 

NO  GRID:  qR  =  7.20  (cm.  H.O) 

GRID  1:  qR  —  10.00  (cm.  H20) 

GRID  2:  qR  =  10.00  (cm.  H20) 

GRID  3:  qR  =  10.00  (cm.  H:0) 

GRID  4:  qR  =  5.40  (cm.  H20) 


The  dynamic  pressure  conditions  for  the  experiment  would  give  a  subcritical  Reynolds 
number  of  Re^l.lXlO5  and  actual  dynamic  pressures  in  the  test  section  shown  in 
Table  3.  Test  conditions  for  individual  runs  are  listed  in  Appendix  B.  Other  variables 
were  altered  to  make  comparisons: 

1.  Angle  of  attack. 

2.  Wing  tail  configuration. 

3.  Afterbody  roll  angle. 

Three  different  VLSAM  model  afterbody  configurations  were  chosen  to  compare  the 
effects  of  turbulence  on  the  VLSAM  model  with  and  without  wings  and  at  two  different 
roll  angles. 

1.  BODY  A  -  Wings  and  tails  at  zero  roll  angle  in  a  "  +  "  configuration. 

2.  BODY  B  -  No  wings  or  tails.  Body  roll  angle  set  at  </>K  =  45°  . 

3.  BODY  C  -  Wings  and  tails  at  <f>R  =  45°  roll  angle  in  a  "X"  configuration. 


Figure  14.  Body  configurations  and  reference  system. 
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In  all  cases,  as  the  afterbody  roll  angle  changes,  the  nose  roll  angle  remains  fixed  to  the 
reference  plane.  BODY  A  and  BODY  C  configurations  were  chosen  to  be  certain  that 
measurements  would  indicate  the  presence  of  vortex  asymmetry  by  any  side  force  mag¬ 
nitude.  The  wings  and  tails  in  the  "  + "  and  "X"  roll  angles  are  the  only  wing  attitudes 
where,  in  the  absence  of  asymmetric  vortices,  zero  side  forces,  rolling  moments,  and 
yawing  moments  occur.  JRef.  43J 

The  goal  of  the  experiment  was  to  measure  the  effects  of  the  grid  generated  turbu¬ 
lence  on  the  high  angle  of  attack  aerodynamics  on  the  VLSAM  model,  in  particular  the 
affect  on  side  force  magnitude  and  direction. 

C.  EXPERIMENTAL  PROCEDURE 

1.  Balance  Calibration 

The  precision  balance  was  calibrated  by  the  NASA  Ames  Calibration  Labora¬ 
tory  which  determined  balance  constants  for  the  data  acquisition  program  and  gave  in¬ 
dividual  channel  or  circuit  accuracy.  The  balance  was  fixed  to  the  model/balance  support 
in  the  test  section  without  the  sleeve  and  VLSAM  model.  After  all  electrical  connections 
were  made  and  checked  the  data  acquisition  hardware  was  energized  and  balance 
excitation  voltage  for  each  channel  was  set  to  5.00  ±  0.003  DC  volts  via  the  signal  con¬ 
ditioner.  The  entire  data  acquisition  system  was  allowed  to  warm  up  and  settle  out  for 
one  hour.  Next,  the  amplifier  was  zeroed  at  a  1000  gain  setting  followed  by  zeroing  each 
balance  circuit  to  0  ±  200  nanovolts  across  the  Wheatstone  bridge.  The  balance  sleeve 
was  fitted  along  with  a  cradle  for  static  weights.  The  data  acquisition  software  was  run 
to  obtain  the  initial  tare  readings.  Finally  calibrated  weights  were  added  to  the  cradle  and 
measured  by  the  balance  to  compare  readings  with  balance  calibration  specifications. 
Normal  and  side  force  readings  were  found  to  be  accurate  to  within  0.05%  of  full  scale 
readings.  After  this  initial  test  was  completed,  the  signal  conditioner  was  never  de¬ 
energized  to  maintain  the  excitation  voltages  to  the  balance  circuits. 

2.  Prerun  Calibration  and  Test 

Prior  to  each  run.  the  data  acquisition  system  was  fully  energized  and  allowed 
to  settle  out  for  half  an  hour.  During  the  settling  period,  the  data  acquisition  system 
was  cycled  repeatedly  by  the  software  to  bring  the  relay  multiplexer,  digital  multimeter, 
and  amplifier  to  operating  temperature  conditions.  Next,  the  static  weight  test  was  re¬ 
peated  with  the  mounted  VLSAM  model,  to  ensure  the  balance  circuits  were  remaining 
within  calibration  accuracy.  If  not,  the  model  was  removed,  the  balance  circuits  were 
zeroed,  and  the  static  weight  test  procedure  completed  without  and  with  the  model. 
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3.  Data  Collection 


When  the  prerun  calibration  results  were  satisfactory,  the  wind  tunnel  was  started  and 
test  section  reference  pressure  set  at  1 10%  of  the  data  run  reference  pressure.  The  model 
pitch  or  angle  of  attack  was  slowly  cycled  from  -10°  to  110°  and  back  at  approximately 
2°  per  second.  This  procedure  tested  the  safety  of  the  support  assembly  with  the  balance 
and  the  model  and  let  the  mechanical  interfaces  settle  out.  The  pitch  angles  where 
noticable  model  vibration  occurred  were  noted. 

The  wind  tunnel  was  put  into  a  no  wind  condition  (zero  dynamic  pressure)  with 
the  motor  running  and  the  model  set  at  -5°  pitch.9  The  data  acquisition  software  re¬ 
corded  the  initial  tare  readings.  The  intial  settling  chamber  temperature  and  barometric 
pressure  were  recorded.  The  wind  tunnel  was  brought  up  to  operating  reference  pressure 
and  the  actual  data  acquisition  began.  Missile  angle  of  attack  was  manually  set  and 
entered  into  the  data  acquisition  program  for  each  set  of  readings.  Test  section  reference 
pressure  was  kept  constant  as  the  model  angle  of  attack  increased  adding  blockage.  The 
microcomputer  reduced  the  balance  readings  to  force  and  moment  values  and  stored  the 
data  on  disk.  At  the  end  of  each  run.  reference  pressure  was  set  at  zero  and  several 
readings  recorded  to  measure  any  shift  between  initial  and  final  readings. 

The  force  and  moment  data  are  translated  into  coefficient  form  by  another 
program  which  also  corrects  for  blockage  and  test  conditions.  Data  were  graphed  for 
comparative  analysis  and  graphs  are  found  in  the  results  section. 

4.  Preliminary  Runs 

Previous  experiments  showed  that  variations  in  nose  roll  angle  would  alter  the 
asymmetric  vortex  structure,  thus  changing  side  force  direction  and  magnitude.  Non¬ 
repeatability  of  data  points  is  another  area  often  encountered  in  earlier  investigations. 
Prior  to  conducting  the  actual  experimental  data  runs,  preliminary  runs  were  made  to 
determine  which  nose  roll  angle  produced  the  largest  side  force  magnitude.  Data  runs 
were  also  conducted  to  examine  the  repeatability  of  data  without  changing  test  condi¬ 
tions.  The  preliminary  runs  would  also  indicate  data  trends  and  would  give  some  degree 
of  confidence  in  the  data  collected  in  the  subsequent  runs. 

Preliminary  test  runs  were  short  and  durations  ran  approximately  30  minutes. 
Temperature  increases  were  small,  from  2°F  to  5°F,  giving  a  smaller  temperature  gradi¬ 
ent  and  allowing  test  section  air  density  to  remain  in  a  narrow  range.  When  the  exper¬ 
iment  proceeded  to  the  actual  test  runs,  test  time  durations  typically  lasted  4.5  hours. 

9  This  conditon  will  be  referred  to  as  the  model's  zero  condition. 


In  order  to  keep  test  section  air  density  as  uniform  as  possible,  the  test  procedure  was 
interrupted  at  least  once  during  the  run  cycle  to  cool  the  air  in  the  wind  tunnel.  Settling 
chamber  temperatures  rose  about  10°F  to  20°F  with  larger  increases  for  coarser  turbu¬ 
lence  generating  grids.  When  the  the  settling  chamber  temperature  increased  to  20°  F 
above  the  start  temperature,  the  wind  tunnel  was  cooled  down  by  air  exchange.  Once 
the  tunnel  had  cooled  down  to  within  5°F  of  the  start  temperature,  the  run  continued 
at  2°  below  the  interrupt  angle  of  attack  to  overlap  readings. 

Other  test  conditions  that  were  beyond  the  control  of  the  experiment  include  the 
turbulence  and  vibrations  produced  by  the  wind  tunnel  itself  during  operation  and  vi¬ 
brations  and  deflections  of  the  model  balance  support. 

D.  EXPERIMENTAL  CORRECTIONS 

Model  blockage  corrections  for  changing  angle  of  attack  were  made  by  equations 
[61,  [7],  and  [S],  from  [Ref.  44[ 


q  =  q.\t(i  +  2r) 

{6} 

U  =  L'm(1  +  c) 

{7} 

where: 

q  =  dynamic  pressure  (lb  ft2) 
qM  =  measured  reference  pressure  (lb  ft:) 

U  =  horizontal  velocity  (ft  see) 

UM  =  measured  horizontal  velocity  (ft  sec) 
r  =  blockage  factor 

and: 

1  model  frontal  area  , ri 

£  =  — - : -  S) 

4  test  section  area 

Model  frontal  area  varies  for  each  body  configuration  and  angle  of  attack  with  minimum 
blockage  at  0°  and  maximum  blockage  at  90°  as  shown  in  Figure  15. 

Blockage  equations  for  each  body  configuration  were  developed  as  a  function  of 
model  angle  of  attack.  These  equations  are  shown  in  Figure  15  and  Table  4,  and  are 
implemented  in  the  data  conversion  program  listed  in  Appendix  E. 
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BLOCKAGE  CORRECTION  FACTOR 


Table  4.  VLSAM  MODEL  BLOCKAGE 


Configuration 

Minimum  Blockage 
(  percent) 

Maximum  Blockage 
(percent) 

BODY  A 

3.10 

9.76 

BODY  C 

3.10 

6.74 

BODY  B 

3.10 

6.37 

0.020 


0.014 

0.012 


0.006 


0.004  • ,  ,  , 
-10 


o  -  + 

BODY: 

t  - 

1.26 

E— 04 

*  ALPHA 

+ 

0.007759 

A-  X 

BODY: 

c  - 

1.01 

E— 04 

*  ALPHA 

+ 

0.007759 

o-  o 

BODY: 

c  *= 

9.08 

E— 05 

*  ALPHA 

+ 

0.007759 
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Figure  15.  Blockage  factors 


III.  RESULTS 


Several  important  factors  must  be  kept  in  mind  as  the  results  are  presented.  These 
factors  are  conclusions  from  previous  high  angle  of  attack  studies  on  slender  bodies  in 
low  turbulence  environments  and  may  apply  to  these  results.  The  nose  fineness  ratio  of 
the  VLSAM  model  is  2.29.  Results  by  Keener  [Ref.  11]  dictate  that  the  vortex  cores  from 
the  nose  at  fineness  ratio  3.5  should  not  be  as  crowded  together  as  in  higher  fineness 
ratio  noses  at  5.0.  Since  the  nose  is  pointed,  no  blunt  nose  effects  are  expected.  From 
equation  [2],  vortices  should  become  asymmetric  when  the  angle  of  attack  is  aAV5;240  for 
the  given  apex  angle  of  the  VLSAM  model.  Flow  separation  throughout  the  length  of 
the  nose  should  remain  laminar  at  all  angles  of  attack  for  Red  <  3.0x10'  [Ref.  11].  Con¬ 
ditions  for  all  data  runs  are  stated  in  the  previous  chapter  and  listed  in  Appendix  B  are 
stated  in  the  previous  chapter.  All  data  represent  steady-state  conditions  and  are  ex¬ 
pressed  in  coefficient  form. 

A.  VARIATIONS  IN  NOSE  ROLL  ANGLE 

Previous  research  showed  that  variations  in  nose  geometry  and  roll  angle  altered  the 
vortex  structure  shedding  from  the  nose  and  afterbody  at  high  angles  of  attack.  The 
objective  of  the  first  set  of  preliminary  runs  was  to  determine  the  nose  geometry  which 
produced  the  highest  side  force  magnitude.  Since  the  model  used  in  the  experiment  used 
only  one  nose,  variations  in  nose  roll  angle  became  the  single  variable  affecting  nose 
geometry.  The  nose  with  the  highest  side  force  magnitude  would  be  fixed  in  the  same 
position  regardless  of  afterbody  configuration  for  all  subsequent  runs.  Both  sets  of  runs 
were  conducted  under  identical  conditions,  Re„  =  1. 33x10s  and  test  section  dynamic 
pressure  at  q  =  22.2  lbs  ft:.  Angle  of  attack  was  varied  from  0°  to  90°  in  10°  and  5° 
increments.  .Model  configuration  was  BODY  A  and  no  grids  were  used.  Nose  roll  angle 
was  changed  in  45°  increments  for  each  run.  Results  for  each  set  of  runs  are  plotted  in 
Figure  16  and  Figure  17. 

In  both  cases,  changes  in  nose  roll  angle  alone  produced  variations  in  side  force 
magnitude  and  direction.  The  plots  also  show  a  some  repeatability  of  the  data  points. 
Nose  eight,  in  both  run  sets,  produced  the  largest  side  force  magnitude  in  the  negative 
direction.  The  angle  of  attack  when  side  force  magnitude  became  noticeable  at  a^30° 
was  the  same  for  all  nose  roll  angles  within  each  set  of  runs,  which  is  also  true  for  the 
angle  of  attack  a^SO0  where  side  force  magnitude  decreases  to  negligible  values. 
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Figure  17.  Side  force  variations  with  nose  roll  angle:  Runs  WUA102  to  WOAS02. 


38 


1 


i 


These  side  force  measurements  suggest  that  vortices  become  asymmetric  at  2.AV^?n;  and 
end  at  acvs:S0o.  which  is  true  for  all  roll  angles.  The  angle  of  attack  increments  are  large. 
Smaller  increments  would  give  better  observations  of  data  trends  and  all  subsequent  data 
collection  runs  are  in  one  degree  increments. 

B.  DATA  REPEATABILITY 

A  recurring  problem  dealing  with  side  force  measurements  from  asymmetric  vortices 
is  the  non-repeatability  of  data  except  in  Keener  [Ref.  231  "here  repeatability  of  results 
was  good.  In  this  experiment,  all  readings  are  an  average  of  ten  separate  measurements 
taken  within  a  40  second  interval,  so  any  signs  of  repeatability  could  indicate  that  the 
data  is  statistically  accurate.  As  shown  in  Figure  16  and  Figure  17,  the  nose  roll  angle 
which  gave  the  maximum  side  force  magnitude  for  all  roll  angles  is  nose  eight,  which  is 
true  for  both  run  sets.  The  objective  of  the  repeatability  runs  was  to  determine  if  any 
repeatability  existed,  and  if  not,  what  degree  of  fluctuation  existed.  After  prerun  cali¬ 
brations,  the  model  was  pitched  to  20°  angle  of  attack  from  O'  at  test  dynamic  pressure. 
Readings  were  taken  at  20°  then  the  dynamic  pressure  reduced  to  zero  and  readings 
taken  again.  The  wind  was  turned  on  and  off  for  readings  as  the  model  remained  at  a 
fixed  angle  of  attack. 
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Figure  18.  Data  Repeatability  Runs 
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After  several  readings,  the  model  was  slewed  to  the  the  zero  condition  for  a  final  cali¬ 
bration  check.  This  procedure  was  repeated  for  45*.  70".  and  90'  angles  of  attack  and 
all  results  are  plotted  in  Figure  IS.  Variation  in  side  force  magnitude  is  obvious  at 
20°,  45°.  and  70°.  Measurements  at  90°  show  negligible  lluctuation.  This  agrees  with 
accepted  side  force  fluctuation  phenomenon.  Statistically,  the  mean  value  of  any  data, 
will  be  improved  if  an  increasing  population  of  readings  are  taken.  Time  and  record 
keeping  limitations  allowed  only  ten  readings  per  channel  to  be  recorded.  Some  degree 
of  experimental  accuracy  may  be  lost  with  the  data  fluctuation,  but  the  data  point  trends 
are  useful  in  making  comparisons  between  different  run  conditions.  Comparing 
Figure  16  and  Figure  17,  the  side  force  follows  a  general  trend  which  will  be  useful  for 
comparing  the  data  runs  discussed  in  the  following  sections. 

C.  BODY  IN  ISOLATION 

The  baseline  test  run  was  made  for  BODY  B  with  no  grids  shown  in  Figure  19. 
The  asymmetric  vortices  begin  at  aAV~2G°  which  concurs  with  the  value  jcav^24°  from 
equation  [2J.  The  unsteady  vortices  begin  to  form  at  auv^S0°  .  Between  these  two  an- 
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Figure  19.  R0B802:  Side  Force  Coefficient. 
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gles  of  attack  the  side  force  changes  magnitude  but  not  direction.  The  data  points 
scatter  suggests  a  strong  possibility  of  fluctuations  in  the  steady  state  side  force  readings. 
The  data  trend  shows  the  maximum  side  force  occurring  at  a^56°  and  Cy  =  4.1  then  ta¬ 
pers  off  to  nearly  Cy  =  0  after  which  it  increases  back  to  Cy  =  -2.0  then  tapers  off  into 
the  unsteady  vortex  region.  Since  there  are  no  wings  or  tails,  the  side  force  is  due  en¬ 
tirely  from  the  nose  and  afterbody  vortices.  The  smaller  angle  of  attack  interval  in  Fig¬ 
ure  19  does  not  show  any  data  smoothness  when  compared  to  Figure  16  and 
Figure  17,  which  were  run  at  larger  increments.  The  data  trends  are  similar  which  illus¬ 
trates  the  dominance  of  the  nose  generated  vortices.  The  side  force  magnitudes  in 
Figure  16  and  Figure  17,  are  just  over  30%  of  those  in  Figure  19  which  may  be  due  to 
wing  vortex  effects  discussed  later  in  the  results  section. 10 

The  BODY  B  configuration  with  Grid  1  generated  turbulence  is  shown  in 
Figure  20.  The  onset  of  asymmetric  vortices  is  delayed  to  aAVs;42°  and  unsteady  asym¬ 
metric  vortices  aLV^S2°  .  In  general  the  side  force  direction  remains  the  same  but  there 
is  a  considerable  reduction  in  magnitude  by  about  77%.  One  possible  cause  of  this  de¬ 
crease  in  side  force  is  the  change  in  positioning  of  the  asymmetric  vortices  caused  by  the 
turbulence.  The  turbulence  length  scale  to  missile  diameter  ratio  is  1.08  and  can  be 
classified  as  vortex  scale  turbulence  which  may  suppress  the  asymmetry  of  the  nose 
vortex  pair.  Perhaps  this  is  the  reason  behind  the  delayed  angle  of  attack  for  the  onset 
of  asymmetric  vortices.  Another  reason  may  be  due  to  turbulence  intensity  since 
GRID  1  turbulence  intensity  is  the  highest. 

Comparing  the  results  from  the  remaining  grids  in  Figure  21  thru  Figure  23,  the 
trends  are  similar  to  GRID  1  and  the  vortex  asymmetry  seems  to  be  depressed  by  longer 
turbulence  scales.  For  GRID  4  where  the  turbulence  length  scale  is  the  smallest  and 
side  force  magnitude  is  the  highest,  aAV  occurs  at  an  earlier  angle  of  attack.  The  data 
trend  in  Figure  23  has  a  closer  similarity  to  Figure  19,  although  side  force  magnitude 
is  reduced  by  55%.  Turbulence  intensity  does  not  seem  to  factor  as  much  as  the  length 
scale  to  control  asymmetric  vortices  since  turbulence  of  any  intensity  might  be  expected 
to  eliminate  asymmetric  vortex  effects  and  especially  in  the  case  of  GRID  I  with  the 
highest  turbulence  intensities,  yet  GRID  1  results  show  the  side  force  is  not  eliminated. 


10  The  axisymmetric  wing  and  tail  arrangement  of  BODY'  A  and  BODY'  B  should  not  con¬ 
tribute  any  side  force  lift  from  the  lifting  surfaces. 
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Figure  20.  RIBS01:  Side  Force  Coefficient. 


ANGLE  OF  ATTACK  (DEGREES) 


Figure  21.  R2B801:  Side  Force  Coefficient. 
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Figure  23.  R4B801:  Side  Force  Coefficient. 
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A  correlation  between  turbulence  length  scale  and  turbulence  intensities  to  the  vortices 
is  not  possible  since  the  turbulence  intensities  for  the  grids  are  not  equal. 

In  previous  research.  WarJlaw  produced  2%  turbulence  intensity  by  a  wire  mesh 
grid  [Ref.  12). H  His  results  showed  that  side  force  distribution  was  not  effected  by  the 
turbulence  and  side  force  magnitude  increased  by  10%.  If  the  turbulence  length  scales 
were  small  enough  to  affect  the  boundary  layer,  the  side  force  magnitude  should  de¬ 
crease.  On  the  other  hand,  larger  turbulence  length  scales  should  also  suppress  side 
forces  as  indicated  by  the  results  in  this  experiment.  Since  the  small  scale  low  intensity 
turbulence  from  GRID  4  does  not  appear  to  suppress  side  force  magnitude,  results  from 
increasing  GRID  4  turbulence  intensity  could  settle  whether  turbulence  intensity  con¬ 
tributes  anything  significant  to  the  vortex  structures. 

In  another  investigation  by  Clark,  he  showed  that  a  single  rod  place  upstream  from 
the  stagnation  point  on  an  ogive  nose  may  have  the  same  effect  as  grid  generated  tur¬ 
bulence.  He  injected  a  thin  stream  of  bubbles  into  the  flowfield  upstream  of  the  model 
which  forced  flow  separation  into  steadiness  [Ref.  45], 12  His  results  indicated  that  when  j 

bubbles  were  injected  at  the  point  of  the  nose,  the  vortices  were  forced  to  remain  in  a 
fixed  separation  condition.  When  Clark  injected  the  bubble  stream  on  the  afterbody, 
there  was  no  effect  on  the  vortex  structure  or  side  force  magnitude.  He  concluded  that 
the  wake  from  the  bubble  probe,  rather  than  the  bubbles  themselves,  forced  the  How  ; 

i 

separation  into  steadiness. 

Conjecture  at  this  point  would  indicate  that  turbulence  length  scales  of  vortex  sbe 
may  have  a  greater  influence  on  the  vortex  structure  than  turbulence  intensity.  The 
smaller  scale  low  intensity  turbulence  docs  not  seem  to  attenuate  the  asymmetric  vortex  I 

structure  since  maximum  side  force  magnitudes  are  nearly  equal,  but  docs  seem  to  cause 
a  steadier  flow  evidenced  by  the  data  scatter. 

D.  WINGS  AND  STRAKES  WITHOUT  TURBULENCE. 

The  addition  of  wings  and  strakes  in  the  BODY  A  and  BODY  C  configurations  ! 

show  changing  trends  in  side  force  direction  and  magnitude.  13  BODY  A  in  Figure  24 
shows  larger  side  force,  a  lower  angle  of  attack  for  the  onset  of  asymmetric  vortices. 


-  I 

11  Red  =  1.5x10s.  No  turbulence  length  scale  given. 

12  Red  =  0.4x10s.  No  turbulence  length  scale  or  intensity  given. 

13  The  following  discussion  centers  on  the  wings  and  strakes  since  the  tails  are  not  expected 
to  contnbute  any  significant  effects. 
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Figure  24.  R0A803:  Side  Force  Coefficient. 
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Figure  25.  R0C802:  Side  Force  Coefficient. 
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Figure  26.  R1A804:  Side  Force  Coefficient. 
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Figure  27.  R2A801:  Side  Force  Coefficient. 
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Figure  28.  R3A801:  Side  Force  Coefficient. 
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Figure  29.  R4A80I:  Side  Force  Coefficient. 
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Figure  31.  R2C801:  Side  Force  Coefficient. 
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Figure  32.  R3C801:  Side  Force  Coefficient. 


Figure  33.  R4C801:  Side  Force  Coefficient. 
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larger  side  force  magnitudes  by  12%  when  compared  to  BODY  B  in  Figure  19,  and 
greater  flow  unsteadiness  seen  in  the  data  scatter.  BODY  C  in  Figure  25,  while  showing 
an  equal  degree  of  unsteadiness,  has  even  lower  side  force  magnitudes  than  BODY  B 
and  pronounced  changes  in  the  direction  of  the  side  forces  seen  from  the  positive  side 
force  values  between  a^l8°  and  as:50°.  These  results  suggest  that  the  addition  of  wings 
do  not  significantly  alter  side  force  magnitude.  The  nose  vortices  continue  to  dominate, 
but  the  wings  and  strakes  do  influence  side  force  direction  which  could  be  caused  by- 
other  flow  mechanisms. 

The  afterbody  vortices  of  the  BODY  B  configuration  are  expected  to  contribute  side 
force  to  the  overall  magnitude  which  raises  questions  as  to  what,  if  any,  are  the  position 
of  afterbody  vortices  in  the  BODY  A  and  BODY  C  configurations.  The  symmetry  of 
BODY  A  and  BODY  C  configurations  are  not  expected  to  contribute  any  side  force  lift 
from  the  wings  and  tail  in  the  absence  of  asymmetric  vortices.  Since  side  force  magni¬ 
tude  and  direction  trends  are  generally  similar  but  not  alike  for  all  cases,  there  must  be 
other  mechanisms  involved  from  the  vortices  generated  by  the  leading  edges  of  the 
strakes  and  wings. 

The  first  possible  effect  is  the  strength  and  positioning  of  the  wing  and  strake 
vortices.  The  nose  vortices  arc  always  expected  to  dominate  any  downstream  vortices 
and  may  drive  the  wing  and  strake  vortices  to  asymmetry  contributing  to  the  overall  side 
force  magnitude.  These  wing  and  strake  vortices  could  be  stronger  at  inducing  side 
forces  than  the  afterbody  vortices  and  may  be  the  reason  why  BODY  A  has  larger  side 
force  magnitudes  than  BODY  B.  The  second  effect  may  be  due  to  the  strake  and  wing 
vortices  pushing  the  nose  vortices  further  away  from  the  body  causing  a  weaker  asym¬ 
metric  vortex  effect  and  reducing  the  overall  side  force  magnitude  which  may  explain 
why  BODY  C  has  smaller  side  force  magnitudes  than  BODY  B  without  turbulence.  It 
may  also  be  a  combination  of  the  two  and  could  be  dependent  on  the  angle  of  attack. 
The  third  effect  may  be  due  to  the  nose  and  afterbody  vortices  causing  side  force  lift  on 
the  leeward  wing.  As  the  nose  vortex  pair  becomes  asymmetric,  the  larger  vortex  could 
roll  up  and  over  the  leeward  wing  causing  a  side  force  from  vortex  lift,  a  phenomenon 
enhanced  by  the  strake  and  low  aspect  ratio  wing.  Side  force  vortex  lift  could  also 
contribute  to  the  rolling  moment.  Finally,  the  wings  and  strakes,  depending  on  roll  an¬ 
gle.  may  act  as  splitter  plates  which  form  a  barrier  between  the  vortex  pairs  preventing 
or  attenuating  an  interaction. 


High  angle  of  attack  experiments  by  Deane  [Ref.  27]  and  Jorgensen  [Ref.  46]  with 
monoplane  configurations  show  no  increase  in  side  force  magnitudes  when  compared  to 
the  body  in  isolation.  Further  studies  with  monoplane  wing  configurations  at  o'  and 
90°  roll  angles  would  be  invaluable  to  study  these  possible  effects.  1  he  disposition  and 
strength  of  the  asymmetric  vortices  would  have  to  be  quantified  to  determine  how  the 
actual  vortex  interactions  occur. 

E.  WINGS  AND  STRAKES  WITH  TURBULENCE. 

The  effects  of  adding  turbulence  to  BODY  A  and  BODY  C  were  similar.  In  both 
cases  Figure  26  through  Figure  33,  the  longer  turbulence  scales  suppressed  the  side 
force  magnitudes,  but.  for  BODY  C  the  data  show  more  unsteadiness  with  the  smaller 
scale  turbulence,  as  opposed  to  the  BODY  B  results.  The  larger  scale  turbulence  delays 
the  angle  of  attack  for  the  onset  of  asymmetric  vortices  for  BODY  C  to  aAVs:23c  The 
results  for  BODY  A  are  again  close  to  BODY  B  with  the  onset  of  asymmetric  vortices 
occurring  earlier  at  aAV^22°  for  the  no  turbulence  condition.  The  effects  of  the  larger 
scale  turbulence  in  delaying  aAV  are  apparent  when  Gif  ID  1  results  in  Figure  26  and 
Figure  30  are  compared  to  GRID  4  results  in  Figure  29  and  Figure  33.  Changes  in  side 
force  magnitude  are  mixed  although  the  general  trend  shows  the  lower  side  force  mag¬ 
nitude  at  longer  turbulence  length  scales  except  for  BODY  A  and  GRID  3  in 
Figure  28  which  has  a  higher  side  force  magnitude.  The  side  force  directional  tendencies 
are  also  suppressed  suggesting  longer  turbulence  length  scales  brings  more  steadiness  to 
the  asymmetric  vortices.  Turbulence  intensity  effects  from  the  graphs  show  that  the 
higher  intensities  suppress  side  force  magnitude  but  do  not  seem  to  change  flow  sepa¬ 
ration  unsteadiness  because  of  the  data  point  fluctuation.  The  general  trend  of  the  data 
points  arc  similar  for  all  cases  examined,  especially  when  comparing  each  separate  body 
configuration  to  each  other.  These  trends  indicate  the  dominating  influence  of  the  nose 
vortices.  The  side  force  differences  between  the  body  configurations  could  be  attested  to 
the  presence  of  the  wings  and  their  roll  angle. 

Vortex  scale  turbulence  may  be  responsible  for  changing  the  position  of  the  vortices 
which  explains  the  changing  direction  and  magnitude  of  the  side  forces  which  is  true  for 
the  case  of  BODY  B  with  GRID  1.  Here  the  vortex  scale  turbulence  may  contribute  to 
the  general  unsteadiness  of  the  vortex  structure  evidenced  by  the  scattered  data  points. 
This  contribution  to  unsteadiness  can  be  compared  to  a  missile  experiencing  unstead¬ 
iness  with  comparable  turbulence  length  scales.  Boundary  layer  scale  turbulence  may 
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exist  due  to  the  cascade  effect  but  does  not  seem  to  factor  into  the  experimental  results 
since  side  force  magnitudes  were  always  larger  for  smaller  scale  turbulence. 

F.  NORMAL  FORCE  RESULTS 

The  normal  force  results  are  plotted  in  Figures  33  through  47  and  show  the  charac¬ 
teristic  trend  of  a  steady  increase  to  a  maximum  where  the  values  level  off.  The  normal 
force  is  largest  for  the  BODY  A  configuration  and  the  smallest  for  BODY  B.  The 
maximum  normal  force  value  of  the  BODY  A  configuration  is  C„  =  32.3  which  is  22% 
higher  than  BODY  C  at  Cn  =  26.4  and  101%  higher  than  BODY  B  at  C„  =  16.1  The 
normal  force  curve  for  BODY  A  is  much  steeper  and  levels  off  at  the  higher  value  and 
lower  angle  of  attack  than  the  other  configurations.  This  may  be  due  to  the  wings 
producing  more  lift  at  lower  angles  of  attack  and  more  drag  at  higher  angles  of  attack 
than  for  BODY  C.  The  addition  of  grid  generated  turbulence  causes  a  decrease  in  slope 
of  the  normal  force  curves.  The  curves  without  turbulence  show  a  sharp  turn  at  a  max¬ 
imum  value  where  the  curve  levels  off.  The  curves  become  rounded  with  the  turbulence. 
There  does  not  appear  to  be  any  significant  effect  on  the  overall  magnitude  of  the 
normal  force  in  the  presence  of  turbulence.  The  variation  in  turbulence  length  scale  and 
intensity  do  not  appear  to  affect  the  normal  force  values  as  in  the  results  for  side  forces. 
There  is  no  noticeable  fluctuations  in  the  trend  of  the  data  points  since  the  normal  force 
does  not  fluctuate  with  or  without  the  presence  of  turbulence. 

G.  SIDE  FORCE  TO  NORMAL  FORCE  RATIO 

The  side  force  to  normal  force  ratio  is  important  to  slender  bodies  with  monoplane 
wings  such  as  cruise  missiles  and  aircraft.  As  this  ratio  increases,  the  slender  body  must 
be  able  to  cope  with  the  increasing  side  forces,  thus  requiring  large  amounts  of  yaw 
control.  Since  yaw  control  does  not  equal  lift  control  due  to  the  available  control  surface 
area  in  monoplane  configurations,  side  forces  can  be  a  problem.  Cruciform  missiles 
should  be  able  to  cope  with  high  side  force  to  normal  force  ratios  since  their  wing  and 
tail  arrangement  gives  equal  lift  capability  in  the  pitch  and  yaw  directions. 

Since  the  normal  force  components  are  not  as  effected  by  the  addition  of  grid  gen¬ 
erated  turbulence,  the  side  force  to  normal  force  ratios  should  follow  the  side  force 
trends.  The  lowest  normal  force  values  of  all  the  tested  body  configurations  was 
BODY  B  which  also  had  a  substantial  side  force  value,  and  therefore  should  be  expected 
to  have  the  highest  side  force  to  normal  force  ratio  which  is  the  case  as  shown  in  the 
plots.  BODY  B  has  a  maximum  ratio  of  Cy Cn  =  -0.369  which  is  112%  higher  than 
BODY  A  at  Cy  Cn  =  -0.174  and  367%  higher  than  BODY  C  at  Cy  C„  =  -0.079  . 
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The  side  forte  to  normal  force  ratios  for  the  turbulence  conditions  in  Figures  5u 
through  63.  should  follow  the  side  force  (with  turbulence)  trends  since  the  normal  force 
components  are  not  affected  to  the  same  degree  as  the  side  force  components.  The 
C\  Cn  plots  show  some  resemblance  to  the  trends  in  the  side  force  plots.  The  wide 
degree  of  scatter,  especially  for  the  turbulence  from  GRID  4  is  representative  of  the 
side  force  trends.  The  cruciform  missile  should  be  able  to  cope  with  higher  side  force  to 
normal  force  ratios  due  to  it's  wing  and  tail  configuration.  From  these  results,  the 
winged  configurations.  BODY  A  and  BODY  C  have  smaller  side  force  ratios  than  the 
body  in  isolation  of  BODY  B.  Although  the  ratio  is  small,  this  may  be  misleading  since 
the  BODY  A  and  BODY  C  configurations  must  be  able  to  cope  with  higher  normal  and 
side  forces  than  BODY  B  which  can  be  easily  accomplished  by  their  cruciform  wing 
configurations. 
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NORMAL  FORCE  COEFFICIENT 


Figure  36.  R2BS01:  Normal  Force  CoelTicient. 
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NORMAL  FORCE  COEFFICIENT 
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Figure  38.  R4BS01:  Normal  Force  CoclTicient. 
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NORMAL  FORCE  COEFFICIENT 
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Figure  40.  R0C802:  Normal  Force  Coefficient. 
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sure  41.  R1A804:  Normal  Force  Coefficient. 
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NORMAL  FORCE  COEFFICIENT 
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Figure  43.  R3A801:  Normal  Force  Coefficient. 
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Figure  44.  R4A801:  Normal  Force  Coefficient. 
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NORMAL  FORCE  COEFFICIENT 
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Figure  45.  R1C802:  Normal  Force  Coefficient. 
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Figure  46.  R2C80I:  Normal  Force  Coefficient. 
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NORMAL  FORCE  COEFFICIENT 


Figure  47.  R3C801:  Normal  Force  Coefficient. 
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Figure  48.  R4C801:  Normal  Force  Coefficient. 
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Figure  53.  R4D801:  Side  Force  to  Normal  Force  Ratio. 
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57.  R3A801:  Side  Force  to  Normal  Force  Ratio. 
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Figure  58.  R4A801:  Side  Force  to  Normal  Force  Ratio. 
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Fissure  59.  RUC802:  Side  Force  to  Normal  Force  Ratio. 
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Figure  61.  R2C801:  Side  Force  to  Normal  Force  Ratio. 
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Figure  62.  R3C801:  Side  Force  to  Normal  Force  Ratio. 
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Figure  63.  R4C801:  Side  Force  to  Norma!  f  orce  Ratio. 


IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  VLSAM  in  the  marine  environment  may  encounter  turbulence  during  the 
launch  phase  at  high  angles  of  attack.  The  turbulence  intensities  in  the  marine  envi¬ 
ronment  are  higher  than  produced  in  the  experiment,  as  are  the  turbulence  length  scales 
which  are  on  on  the  order  of  Lu:Ld  ss  280: 1  for  a  typical  VLSAM.  Cascading  and  or  ship 
airwake  may  produce  smaller  low  intensity  turbulence  on  the  order  of  missile  diameter 
or  vortex  size.  These  turbulence  scales  relative  to  the  model  diameter  may  be  produced 
in  the  wind  tunnel. 

Subsonic  wind  tunnel  investigations  were  performed  on  a  cruciform  wing  missile 
model  with  different  grid  generated  turbulence  length  scales  and  intensities.  The  model 
was  tested  with  three  different  body  configurations  with  and  without  wings  at  0C  to  90° 
angle  of  attack.  The  results  from  the  experimental  runs  with  no  turbulence  yield  the 
following  conclusions'. 

1.  Side  force  magnitude  and  direction  change  with  variations  in  nose  roll  angle. 

2.  There  is  a  general  trend  in  the  repeatability  of  data  points  and  to  some  degree  for 
all  data  points  (for  the  steady  state  condition)  which  suggests  the  dominating  influence 
of  the  nose  generated  vortices. 

3.  Strake  and  wing  vortices  may  play  an  important  role  to  influence  the  disposition 
of  nose  generated  vortices.  While  the  low  pressure  on  the  top  of  the  wing  may  pull  the 
nose  vortices  closer  to  the  body,  leading  edge  vortices  from  strakes  and  wings  may  push 
them  away. 

4.  The  maximum  side  force-to-normal  force  ratio  was  observed  for  the  body  in  iso¬ 
lation  configuration,  although  normal  and  side  force  values  were  actually  lower  than  the 
winged  configurations. 

The  results  of  the  turbulence  experiments  indicate  that  turbulence  effects  for  side 
forces  may  vary  with  length  scale  and  intensity.  Two  turbulence  length  scales,  boundary 
layer  scale  turbulence  and  vortex  scale  turbulence,  are  particularly  important  as  the 
former  may  have  large  effect  on  flow  separation  mechanisms  which  generate  the  asym¬ 
metric  vortices  and  the  latter  on  the  disposition  of  the  vortices  in  the  flowfield.  Turbu¬ 
lence  results  conclude: 

1.  Longer  turbulence  length  scales  decrease  side  force  magnitude  and  possibly  any 
fiow  unsteadiness  in  the  asymmetric  vortex  structure. 
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2.  The  angle  of  attack  for  the  onset  of  asymmetric  vortices  is  postponed  by  longer 
turbulence  length  scales.  Since  these  length  scales  are  quite  large  when  compared  to  the 
boundary  layer,  the  asymmetric  vortex  structure  may  be  suppressed  by  the  interaction 
between  the  inviscid  vortex  fiowfield  and  the  vortex  scale  free  stream  turbulence. 

3.  Higher  turbulence  intensities  cause  less  flow  unsteadiness  for  missiles  with  wings. 

4.  Turbulence  had  little  effect  on  the  normal  force  values.  Maximum  lift  values  were 
nearly  identical  and  there  was  no  fluctuation  in  normal  force  values  for  all  turbulence 
conditions. 

Further  recommendations  for  a  continued  study  of  wing  and  turbulence  interactions 
include: 

1.  Perform  flow  visualization,  hotwire,  surface  flow,  laser  doppler,  and  water  tunnel 
tests  to  quantify  the  various  flowfields  about  the  missile  model  at  high  angles  of  attack. 

2.  Test  monoplane  wing  configurations  at  0°  and  90c  roll. 

3.  Compare  the  effects  of  turbulence  length  scales  at  the  same  intensities  and  with 
longer  and  shorter  turbulence  length  scales  than  those  used  in  this  experiment. 

4.  Compare  the  effects  of  turbulence  intensities  at  the  same  length  scale  and  with 
higher  intensities  than  those  used  in  this  experiment. 

5.  Model  the  rocket  motor  plume  to  study  the  interaction  with  flowfields  about  the 
missile. 
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APPENDIX  A.  BALANCE  CALIBRATION  CONSTANTS 
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Figure  64.  Balance  Calibration  Constants:  [Ref.  42] 
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=  0.00006*00 
=-1.20726-05 
=  2.73256-06 
=-6. 22176-03 
=  0,00  30  6*00 


RM/N2*92*  =-1.15126-07 
RM/4  *4  ♦  =  0.00306*00 
7M/S1*S1*  =  0.00006*00 
RH/$2*S2*  *  5.15605-08 


=  5.263  0  F  —  0 2 
«  3.71765-03 
»  5.26195-03 
»- 7 . 29156-03 
=- 5. 35636-03 

*  5.25815-0? 
*-3.5945^-04 

*  o.oooo?*oc 

*  o.ooon:*oo 

*  0.00036*00 
0.  0000  6*03 

r. 94995-07 
0.  30  0-5*0? 
1.  0  670  c  -  36 


Sl/PM- 
S2/91- 
S2/N2- 
S2/4  - 
S2/S1- 
S2/9H- 
P9/91- 
P  9/ 92- 
P9/4  - 
P8/S1- 
P9/S2-  = 

91/92*92- 
91/4  *4  - 
91/SI*S1- 
9 1/S  2*3’- 
91/R  «**M- 

92/>!l*Jl- 
92/5  *4  - 
92/Sl*Sl- 
92/52*52- 
92/K*t#9M- 
4  /  9 1  *  9 1  - 
4  /N  2*92  - 
4  /31*il- 
4  /3  2* S 1  - 
4  /P8*R“- 
S 1/71 1*91- 
Sl/9  2*92  - 
Sl/4  *4  - 
51/S2*S2- 
Sl/P9*3“- 
52/91*91- 
32/92*9'- 
S2/4  *4  - 
S2/S1*S1- 
S2/R  8*2»- 
R*/Nl*91 - 
RM/N2*9? - 
R»'/4  *4  - 
RM/S1*S1- 

RP./SZ4S2- 


=  3.00005*00 

=  3.  '■3206-04 
=-=.2:975-06 
r-1.  3*  576-05 
=  4.94936-07 
=  0.  00  00  6*00 
=  1.17735-03 
=  4.2547--06 
=- 4 . 59465-06 
=  0.  03  0C5*C0 
=  D. 0000500 
=  7.  5"  016  —  04 
=  1. 2723" -05 
=  0. 00  O'  6  ♦  00 
=  4.33456-35 
=  0.0)302*00 
«  5.39695-34 
=  3.2113  6-37 
=  5. 6  2  3  5  =  -  06 
=  -3.  I'  545-05 
=  -9.  98  30  2-06 
=-S. 0  075-34 
=-1.5 •975-37 
=  3.00006*00 

*  5.00006*00 
=  0.OC025*0O 

*  0.03006*00 
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Figure  65.  Balance  Calibration  Constants:  [Ref.  42] 
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APPENDIX  B.  RUN  MATRIX 

The  run  matrix  lists  all  runs  conducted  for  this  thesis  research.  All  runs  were  con¬ 
ducted  by  the  author  in  the  Naval  Postgraduate  School  Wind  Tunnel  Test  Facility.  Not 
all  of  the  runs  were  discussed  in  the  results  section  and  were  essential  trial  runs  to  fa¬ 
miliarize  the  author  with  the  data  acquisition  system  and  observe  the  general  trend  of 
the  data. 

These  runs  are  listed  by  a  sequence  code  as  in: 

R0B802 

where: 

1.  The  first  letter  is  a  run  code. 

"W"  -  Nose  roll  angle  runs. 

"E"  -  Data  repeatability  run. 

"R"-  Data  collection  run. 

2.  The  next  digit  is  the  grid  number 

0  -  NO  GRID 

1  -  GRID  I 

2  -  GRID  2 

3  -  GRID  3 

4  -  GRID  4 

3.  The  next  digit  is  the  nose  number. 

Nose  1  is  at  0°  roll  angle. 

Each  nose  represents  a  45°  roll  angle  increment. 

4.  The  last  two  digits  are  the  run  number. 
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RUN 

BODY 

SCREEN 

X 

n 

Is. 

DYNAMIC  PRESSURE 
(Ills /IV) 

VELOCITY 
(ft /see) 

(Jo) 

JHfLLVHTdfVTL 

ANCLE  OF  ATTACK 
(degrees) 

INCREMENT 

(degrees) 

WO  A 1 0 1 

A 

0 

1.33 

22.2 

130.2 

68.0 

0—90 

10.5 

W'O  A  201 

D 

0 

1.33 

WISM 

68. 0 

0— »9o 

10.5 

WO  A  301 

IB 

aaa 

130.2 

68.0 

0— *90 

10.5 

W0A401 

A 

A 

1.33 

22.2 

130.2 

68.0 

(>-♦90 

10.5 

WOA50 1 

1.33 

22.2 

130.2 

68.0 

0-*90 

10.5 

W0A6OI 

A 

A 

AA 

130.2 

6S.0 

0-*90 

10.5 

W0A701 

22.2 

130.2 

68.0 

0— *90 

10.5 

W'OASOl 

A 

0 

1.33 

22.2 

130.2 

68.0 

0—90 

lo.5 

W0A101 

D 

0 

1.33 

ESS 

130.2 

68.0 

0— *9u 

lo.5 

WO  A  20 1 

A 

D 

1.33 

B 

130.2 

68.0 

0—90 

■HUB 

WUA301 

na 

B 

130.2 

68.0 

0—90 

10.5 

Wo  A  401 

1.33 

B 

130.2 

68.0 

0—90 

B 

WO  A  501 

D 

1.33 

B 

130.2 

68.0 

0—90 

B 

W0A601 

A 

0 

. 

1.33 

22.2 

130.2 

68.0 

0-90 

10.5 

WO  A  70 1 

A 

0 

1.33 

->•>  ■> 

130.2 

68.0 

0—90 

10.5 

WO  A  SOI 

D 

0 

1.33 

1 

130.2 

.  _ 

68.0 

0—90 

10.5 

EOAS251 

B 

. 

0 

1.13 

1 5.  S' 

. 

109.9 

62.5 

25 

1:0  A  8451 

A 

1.13 

15.8 

110.9 

63.5 

45 

A 


0 


1.13 


LnASTol 


1 5.  S 


109.9 


.5 
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EOA8901 


R0AS01 


R0AS02 


R1A801 


R1A802 


R1A803 


R1AS04 


R2A801 


R3AS01 


R4AS01 


ROC SO I 


ROCS 


R1CS01 


R1CS02 


R2C801 


R3CS01 


R4CS01 


ROB  SO  I 


ROB So 2 


R1BSU1 


R2BS01 


R3B801 


R4BS01 


A  0 


A  0 


110.1 

64.0 

110.6 

66.3 

1 10.6 

66.5 

1 10.5 

68.0 

109.9 

79.0 

109.2 

72.5 

109.5 

75.0 

110.2 

72.0 

107.7 

74.0 

112.7 

71.0 

109.4 

66.ii 

110.S 

68.0 

110.1 

64.0 

109.1 

71.0 

109.0 

_  .  .... 

70.5 

107.5 

72.0 

113.0 

74.0 

109.6 

67.0 

1 10.0 

63.0 

1 10.2 

<0.5 

108.9 

68.5 

107.4 

71.0 

112.7 

71.0 

109.2 

64.0 

90 


0~*92 


•  1  o— ►  100 


-5-95 


0-92 


92— -4 


-5-95 


-95 


-5—95 


0—90 


-5-95 


-5-95 


APPENDIX  C,  NUMERICAL  DATA  LISTING 


The  following  data  was  collected  during  the  data  collection  runs.  A  complete  list  of 
all  the  runs  and  respective  experimental  conditions  is  found  in  the  RUN  MATRIX  of 
Appendix  B. 
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RUNFILE>  W0A101 


AOA 

Cn 

Cy 

Cy/Cn 

Cax 

Cpt 

Crl 

Cyw 

***** 

****** 

****** 

****** 

****** 

****** 

**-. 

,-*** 

0.  0 

-0. 

091 

-0. 

050 

0.  554 

-0. 

009 

0.  001 

0. 

000 

-0. 

002 

1.  0 

-0. 

137 

-0. 

060 

0.  440 

-0. 

010 

0.  000 

0. 

000 

-0. 

002 

2.  0 

-0. 

144 

-0. 

062 

0.  430 

-0. 

021 

-0.  001 

0. 

000 

-0. 

002 

10.  0 

2. 

600 

-0. 

398 

-0.  153 

0. 

238 

-0. 023 

0. 

,  031 

-0. 

010 

20.  0 

8. 

254 

-0. 

024 

-0. 003 

0. 

282 

-0.  061 

-0. 

,  006 

-0. 

013 

30.  0 

15. 

948 

0. 

195 

0.  012 

0. 

139 

0.  064 

-0. 

,  051 

-0. 

020 

35.  0 

20. 

241 

0. 

271 

0.  013 

-0. 

021 

0.  214 

-0. 

.  066 

-0. 

019 

40.  0 

23. 

307 

-1. 

338 

-0. 057 

-0. 

132 

0.  420 

0. 

,  057 

0. 

038 

45.  0 

26. 

572 

-2. 

147 

-0. 081 

-0. 

162 

0.  497 

0. 

,  103 

0. 

039 

50.  0 

26. 

997 

-2. 

990 

-0. Ill 

-0. 

168 

0.  523 

-0. 

.  141 

-0. 

013 

55.  0 

27. 

879 

-0. 

432 

-0. 016 

-0. 

333 

0.  530 

-0. 

.  062 

0. 

054 

60.  0 

28. 

805 

-1. 

807 

-0. 063 

-0. 

489 

0.  661 

-0, 

.  502 

0. 

051 

65.  0 

27. 

747 

-2. 

249 

-0. 081 

-0. 

528 

0.  651 

-0 

.  553 

0. 

053 

70.  0 

26. 

748 

-2. 

026 

-0. 076 

-0. 

620 

0.  459 

-0 

.  710 

0. 

046 

80.  0 

26. 

091 

-0. 

634 

-0. 024 

-0. 

766 

0.  138 

0 

.  006 

0. 

035 

90.  0 

25. 

256 

-0. 

351 

-0. 014 

-0. 

.  227 

-0. 781 

0 

.  030 

0. 

012 

RUNFILE> 

AOA 

W0A201 

Cn 

Cy 

Cy/Cn 

****** 

Cax 

Cpt 

Crl 

Cyw 

****** 

•Mririrte sfr 

****** 

VoWcVc'-Zr* 

****** 

****** 

0.  0 

-0. 065 

-0. 044 

0.  686 

-0. 009 

0.  000 

0.  001 

-0. 003 

0.  0 

-0. 306 

-0. 545 

1.  784 

0.  255 

-0. 002 

0.  052 

-0. 010 

10.  0 

2.  458 

-0. 374 

-0.  152 

0.  216 

-0. 023 

0.  029 

-0. 016 

20.  0 

8.  685 

-0. 709 

-0. 082 

0.  213 

-0. 077 

0.  003 

-0. 029 

30.  0 

16. 323 

-0. 741 

-0. 045 

0.  095 

0.  046 

-0. 055 

-0. 041 

35.  0 

20. 455 

-0. 845 

-0. 041 

-0. 063 

0.  219 

-0. 010 

-0. 058 

40.  0 

24. 255 

-1. 379 

-0. 057 

-0. 157 

0.  394 

0.  007 

-0. 100 

45.  0 

26. 546 

-0. 255 

-0. 010 

-0. 146 

0.  503 

0.  040 

-0. 107 

50.  0 

27. 550 

1.  616 

0.  059 

-0. 181 

0.  527 

0.  093 

-0. 071 

55.  0 

28. 214 

-1.  099 

-0. 039 

-0. 341 

0.  544 

0.  058 

-0.  068 

60.  0 

28.  591 

0.  539 

0.  019 

-0.  419 

0.  638 

0.  655 

-0. 059 

65.  0 

28. 311 

-1.  821 

-0. 064 

-0. 561 

0.  666 

0.  093 

-0. 075 

70.  0 

26. 655 

-1.  068 

-0. 040 

-0. 617 

0.  436 

-0. 106 

-0. 061 

80.  0 

26. 367 

-0. 946 

-0. 036 

-0.  758 

0.  084 

0.  026 

-0.  088 

90.  0 

25.  102 

-1.  091 

-0.  043 

-0.  150 

-0. 885 

0.  064 

-0. 102 
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RUNFILE>  W0A301 


AOA 

Cn 

iclrfrtrk 

*** 

*** 

0. 

0 

-0. 

043 

0. 

0 

-0. 

215 

10. 

0 

2. 

601 

20. 

0 

8. 

595 

30. 

0 

16. 

874 

35. 

0 

21. 

081 

40. 

0 

25. 

104 

45. 

0 

27. 

283 

50. 

0 

28. 

326 

55. 

0 

29. 

232 

60. 

0 

29. 

462 

65. 

0 

28. 

826 

70. 

0 

27. 

417 

80. 

0 

26. 

820 

90. 

0 

25. 

566 

RUNFILE>  W0A401 
AOA  Cn 

I'nHrinV  ****** 


17. 113 
21.  443 
24. 944 
27.  542 
28. 297 


G  •] 


28. 311 
27.  768 
26. 989 
25. 423 


cy 

Cy/Cn 

****** 

****** 

-0. 

.  004 

0. 

083 

-0. 

.  501 

2. 

332 

-0. 

.  287 

-0. 

110 

-0 

.  175 

-0. 

020 

-1, 

.  122 

-0. 

067 

-0, 

.  215 

-0. 

010 

-2. 

.  107 

-0. 

084 

-0 

.  669 

-0. 

025 

-0 

.  274 

-0. 

010 

-3. 

.  079 

-0. 

105 

-3. 

.  712 

-0. 

126 

-1. 

.  684 

-0. 

058 

-4. 

.  210 

-0. 

154 

-1. 

.  355 

-0. 

051 

-1. 

.  463 

-0. 

057 

Cy  Cy/Cn 

iWWWWWp 


figs  *j 


0.  059 


-0.580  -0.021 

-0.770  -0.030 


Cax 

Cpt 

****** 

**i 

WrVnV 

0. 

002 

-0. 

001 

0. 

270 

-0. 

003 

0. 

242 

-0. 

024 

0. 

246 

-0. 

066 

0. 

125 

0. 

106 

0. 

015 

0. 

262 

-0. 

149 

0. 

485 

-0. 

168 

0. 

574 

-0. 

204 

0. 

593 

-0. 

370 

0. 

624 

-0. 

511 

0. 

707 

-0. 

601 

0. 

721 

-0. 

655 

0. 

493 

-0. 

786 

0. 

180 

-0. 

277 

-0. 

730 

Cax  Cpt 

******  ****** 

-0.004  -0.002 

0.282  -0.002 


g 


c 


Crl  Cyw 


iWWrsteWr 

**-. 

0. 

000 

0. 

000 

0. 

062 

-0. 

010 

0. 

041 

-0. 

005 

-0. 

01C 

0. 

001 

0. 

012 

0. 

132 

-0. 

043 

0. 

176 

0. 

013 

0. 

225 

-0. 

010 

0. 

178 

-0. 

014 

0. 

191 

-0. 

016 

0. 

201 

-0. 

084 

0. 

182 

0. 

095 

0. 

150 

-0. 

565 

0. 

133 

0. 

016 

0. 

093 

0. 

097 

0. 

059 

Cr 1  Cyw 

******  ***** 


C.  001 


0.  090 


Sllfaralllww 


-0.  774 
-0. 138 


0.  087 
-0. 933 


RUNFILE>  W0A501 


AOA 

Cn 

Cy 

Cy/Cn 

Cax 

Cpt 

Crl 

jrM&frk 

****** 

•teMrSrSrf r 

ibHbhWf 

****** 

****** 

0.  0 

0.  000 

0.  027 

-96. 508 

0.  002 

-0. 002 

-0. 002 

0.  0 

-0. 030 

-0. 222 

7.  336 

0.  300 

-0.  003 

0.  049 

10.  0 

2.  778 

-0. 138 

-0.  050 

0.  271 

-0.  027 

0.  037 

20.  0 

9.  184 

-0. 235 

-0.  026 

0.  284 

-0. 102 

0.  012 

30.  0 

16. 839 

0.  779 

0.  046 

0.  186 

0.  043 

-0. 048 

35.  0 

21.  164 

-0.  141 

-0. 007 

0.  048 

0.  217 

-0. 099 

40.  0 

25. 155 

2.  221 

0.  088 

-0. 073 

0.  400 

-C. 179 

45.  0 

27. 802 

1.  814 

0.  065 

-0. 081 

0.  496 

-0. 065 

50.  0 

28. 500 

2.  378 

0.  083 

-0.  171 

0.  507 

0.  117 

55.  0 

29. 545 

-1.  806 

-0. 061 

-0. 340 

0.  586 

-0. 118 

60.  0 

29. 907 

0.  590 

0.  020 

-0. 465 

0.  691 

0.  598 

65.  0 

29.  149 

1.  626 

0.  056 

-0. 569 

0.  672 

0.  535 

70.  0 

27.  617 

1.  542 

0.  055 

-0. 661 

0.  433 

0.  760 

80.  0 

27.  235 

-0.  518 

-0.  019 

-0.  804 

0.  133 

0.  032 

90.  0 

25. 562 

-0. 764 

-0.  030 

-0. 189 

-0. 848 

0.  070 

RUNFILE>  W0A601 


AOA 

Cn 

Cy 

Cy/Cn 

Cax 

Cpt 

****** 

Crl 

***** 

****** 

****** 

****** 

****** 

0.  0 

0.  029 

0.  004 

0.  124 

0.  007 

0.  001 

0.  001 

0.  0 

-0. 081 

-0. 479 

5.  884 

0.  287 

-0. 007 

0.  054 

10.  0 

3.  009 

-0. 225 

-0. 075 

0.  293 

-0.  021 

0.  026 

20.  0 

8.  930 

-0.  255 

-0. 029 

0.  319 

-0. 077 

0.  016 

30.  0 

16.  993 

0.  406 

0.  024 

0.  177 

0.  062 

-0. 067 

35.  0 

21. 059 

-0.  914 

-0. 043 

0.  006 

0.  239 

-0. 099 

40.  0 

24. 996 

1.  930 

0.  077 

-0. 070 

0.  376 

-0. 125 

45.  0 

26. 694 

4.  708 

0.  176 

-0. 076 

0.  481 

-0. 130 

50.  0 

27. 084 

3.  092 

0.  114 

-0. 165 

0.  457 

0.  211 

55.  0 

28. 965 

1.  533 

0.  053 

-0. 342 

0.  429 

0.  865 

60.  0 

28. 538 

1.517 

0.  053 

-0. 419 

0.  678 

0.  805 

65.  0 

2  7.  466 

0.  793 

0.  029 

-0. 507 

0.  617 

0.  785 

70.  0 

27.  061 

1.  276 

0.  047 

-0. 628 

0.  482 

0.  824 

80.  0 

26. 666 

-0. 638 

-0. 024 

-0.  779 

0.  114 

0.  042 

90.  0 

25. 381 

-0.  832 

-0. 033 

-0. 194 

-0. 8S1 

0.  075 

7M 


Cyw 

Wr 

0.  002 
0.  001 
-0. 009 
-0. 022 
-0. 005 
-0.  041 
-0.  091 
-0. 123 
-0. 065 
-0. 063 
-0. 040 
-0. 039 
-0. 099 
-0. 090 
-0. 114 


Cyw 

****** 
0.  002 
-0. 004 
0.  006 
-0. 001 
-0. 031 
0.  006 
-0. 039 
-0. 030 
0.  029 
0.  053 
0.  055 
0.  0^3 
0.  036 
0.  023 
-0. 006 


RUNFILE>  W0A701 


AOA 

Cn 

Cy 

Cy/Cn 

Cax 

Cpt 

Crl 

Cyw 

****** 

**■ 

jh'rtWr 

****** 

WnhWdf 

****** 

** 

sww 

0.  0 

-0. 

064 

-0. 

.  002 

0. 

031 

-0. 

002 

-0. 002 

0. 

001 

0. 

001 

0.  0 

0. 

063 

-0. 

.  294 

-4. 

671 

0. 

281 

-0.  006 

0. 

047 

-0. 

004 

10.  0 

3. 

082 

-0, 

.  165 

-0. 

053 

0. 

272 

-0. 019 

0. 

032 

0. 

003 

20.  0 

8. 

978 

-0 

.  295 

-0. 

033 

0. 

286 

-0. 066 

-0. 

002 

0. 

015 

30.  0 

16. 

764 

0, 

.  359 

0. 

021 

0. 

172 

0.  054 

-0. 

046 

0. 

038 

35.  0 

20. 

937 

0 

.  199 

0. 

00  9 

-0. 

009 

0.  210 

-0. 

043 

0. 

012 

40.  0 

24. 

792 

1 

.  287 

0. 

052 

-0. 

093 

0.  418 

-0. 

141 

0. 

068 

45.  0 

26. 

860 

-3 

.  952 

-0. 

147 

-0. 

103 

0.  499 

0. 

119 

0. 

111 

50.  0 

27. 

828 

-3 

.431 

-0. 

123 

-0. 

147 

0.  502 

-0. 

106 

0. 

083 

55.  0 

28. 

785 

-1 

.  112 

-0. 

039 

-0. 

325 

0.  528 

-0. 

266 

0. 

061 

60.  0 

29. 

201 

-2. 

.  495 

-0. 

085 

-0. 

451 

0.  658 

-0. 

641 

0. 

064 

65.  0 

28. 

198 

-2 

.  276 

-0. 

081 

-0. 

535 

0.  662 

-0. 

458 

0. 

061 

70.  0 

27. 

327 

-2 

.  445 

-0. 

089 

-0. 

612 

0.  451 

-0. 

728 

0. 

052 

80.  0 

26. 

521 

-0. 

.  468 

-0. 

018 

-0. 

755 

0.  114 

0. 

018 

0. 

045 

90.  0 

25. 

352 

-0 

.  624 

-0. 

025 

-0. 

157 

-0. 874 

0. 

051 

0. 

022 

RUNFILE> 

AOA 

W0A801 

Cn 

Cy 

Cy/Cn 

Cax 

Cpt 

****** 

Crl 

Cyw 

***** 

***** 

****** 

fei triririrk 

l/HWrsWWf 

****** 

0.  0 

-0. 009 

-0. 005 

0.  592 

0.  002 

-0. 001 

0.  000 

0.  000 

0.  0 

-0. 073 

-0. 441 

6.  017 

0.  269 

-0. 010 

0.  049 

-0. 007 

10.  0 

2.  959 

-0. 365 

-0. 123 

0.  263 

-0. 015 

0.  033 

-0. 003 

20.  0 

9.  295 

-0.  595 

-0. 064 

0.  295 

-0.  074 

0.  027 

0.  037 

30.  0 

17. 193 

-0.  725 

-0. 042 

0.  173 

0.  051 

-0. 018 

0.  049 

35.  0 

21.435 

-1.  701 

-0. 079 

-0. 003 

0.  260 

0.  047 

0.  069 

40.  0 

25. 132 

-2. 326 

-0. 093 

-0. 097 

0.  422 

0.  116 

0.  114 

45.  0 

27.  504 

-4.  127 

-0. 150 

-0. 0S3 

0.  511 

0.  125 

0.  086 

50.  0 

28. 173 

-5. 973 

-0. 212 

-0. 182 

0.  481 

'0. 127 

0.  108 

55.  0 

29. 492 

-5. 202 

-0. 176 

-0. 352 

0.  493 

-0. 761 

0.  084 

60.  0 

29.473 

-2.  778 

-0. 094 

-0. 449 

0.  724 

-0. 838 

0.  098 

65.  0 

27.  741 

-3.  041 

-0.  110 

-0.  523 

0.  55  3 

-0. 768 

0.  081 

70.  0 

27. 842 

-3. 034 

-0. 109 

-0. 656 

0.  481 

>0. 790 

0.  090 

80.  0 

26.  911 

-0. 946 

-0. 035 

-0.  798 

0.  105 

0.  028 

0.  080 

90.  0 

25.  780 

-0.  775 

-0. 030 

-0. 163 

-0.  912 

0.  071 

0.  058 

SO 


RUNFILE>  W0A102 


AOA 

Cn 

cy 

Cy/Cn 

Cax 

Cpt 

Crl 

Cyw 

frffWnV 

****** 

****** 

****** 

****** 

0.  0 

0. 

019 

-0. 

Oil 

-0. 

568 

0. 

003 

0.  001 

0. 

001 

0. 

001 

0.  0 

-0. 

118 

-0. 

413 

3. 

512 

0. 

270 

0.  004 

0. 

054 

0. 

002 

10.  0 

3. 

004 

-0. 

092 

-0. 

031 

0. 

245 

-0. 028 

0. 

,  022 

-0. 

004 

20.  0 

8. 

933 

'0. 

711 

-0. 

080 

c. 

229 

-0. 057 

0. 

.  006 

-0. 

016 

30.  0 

17. 

164 

-0. 

699 

-0. 

041 

0. 

090 

0.  053 

-0. 

.  072 

-0. 

014 

35.  0 

21. 

827 

'0. 

7  84 

-0. 

036 

-0. 

019 

0.  222 

-0. 

.  023 

-0. 

013 

40.  0 

25. 

762 

-3. 

610 

-0. 

140 

-0. 

126 

0.  435 

0. 

185 

0. 

018 

45.  0 

28. 

625 

-4. 

455 

-0. 

156 

-0. 

092 

0.  511 

0, 

.  031 

0. 

022 

50.  0 

30. 

416 

-4. 

048 

-0. 

133 

-0. 

188 

0.  560 

-0. 

.  147 

-0. 

027 

55.  0 

32. 

201 

-2. 

69  3 

-0. 

084 

-0. 

380 

0.  572 

-0. 

.  614 

0. 

057 

60.  0 

32. 

570 

-2. 

894 

-0. 

089 

-0. 

509 

0.  724 

-0. 

.  676 

0. 

074 

65.  0 

31. 

178 

-2. 

576 

-0. 

083 

-0. 

553 

0.  659 

-0. 

.  778 

0. 

060 

70.  0 

31. 

233 

-2. 

337 

-c. 

075 

-0. 

674 

0.  483 

-0. 

.  801 

0. 

062 

80.  0 

31. 

019 

-0. 

583 

-0. 

019 

-0. 

808 

0.  065 

-0, 

.  006 

-0. 

047 

90.  0 

29. 

868 

-1. 

006 

-0. 

034 

-0. 

054 

-1.  114 

0 

.  064 

-0. 

055 

0.  0 

0. 

861 

-0. 

8h4 

-0. 

981 

0. 

131 

-0. 009 

0 

.  043 

0. 

001 

RUNFILE> 

AOA 

W0A202 

Cn 

Cy 

Cy/Cn 

'jfWnWnW? 

Cax 

Cpt 

Crl 

Cyw 

***** 

****** 

****** 

****** 

VrvViWrVnV 

0.  0 

-0.  061 

-0. 252 

4.  145 

0.  273 

-0. 003 

0.  047 

-0. 003 

10.  0 

2.  879 

-0. 369 

-0.  128 

0.  253 

-0. 034 

0.  021 

-0. 009 

20.  0 

8.  783 

-0. 237 

-0.  027 

0.  256 

-0. 086 

-0. 006 

0.  009 

30.  0 

17.  026 

*1.  090 

-0. 064 

0.  130 

0.  035 

-0. 002 

0.  025 

35.  0 

21.  632 

-0. 548 

-0. 025 

-0. 013 

0.  198 

-0. 048 

0.  046 

40.  0 

25.  550 

-0. 402 

-0. 016 

-0. 106 

0.  402 

-0. 030 

0.  013 

45.  0 

28.  532 

1.  469 

0.  051 

-0. 125 

0.  518 

-0. 022 

0.  045 

50.  0 

30. 303 

1.  282 

0.  042 

-0. 237 

0.  545 

0.  110 

0.  06b 

55.  0 

31.  328 

0.  011 

C.  000 

-0. 437 

0.  624 

0.  269 

0.  074 

60.  0 

31.  976 

0.  674 

0.  021 

-0. 572 

0.  737 

0.  498 

0.  080 

65.  0 

32.  018 

-0. 370 

-0. 012 

-0. 707 

0.  722 

0.  098 

0.  081 

70.  0 

30. 740 

-3. 143 

-0. 102 

-0.  7t>2 

0.  482 

-0. 535 

0.  062 

80.  0 

30. 829 

-1.  011 

-0. 033 

-0. 882 

0.  039 

0.  007 

0.  057 

90.  0 

28. 345 

-0. 837 

-0. 029 

-0. 113 

-1.  120 

0.074 

0.  026 

0.  0 

0.  333 

-0. 053 

-0. 158 

0.  032 

-0. 028 

0.  012 

-0. 005 

81 


RUNFILE>  W0A302 


AOA 

Cn 

cy 

Cy/Cn 

Cax 

Opt 

Crl 

Cyw 

*** 

** 

****** 

****** 

**** 

****** 

****** 

****  .* 

**: 

V*** 

0. 

0 

0.  030 

0, 

.  017 

0. 

562 

0. 

010 

0.  002 

0, 

.  000 

-0. 

002 

0. 

0 

-0.  018 

-0 

.  381 

21. 

059 

0. 

278 

0.  004 

0, 

.  056 

0. 

000 

10. 

0 

2.  835 

-0. 

.  396 

-0. 

140 

0. 

240 

-0. 036 

0. 

.  023 

0. 

004 

20. 

0 

8.  893 

-0 

.  226 

-0. 

025 

0. 

257 

-0. 083 

0. 

.  010 

0. 

010 

30. 

0 

17. 045 

0, 

.452 

0. 

027 

0. 

100 

0.  040 

-0, 

,  114 

0. 

026 

35. 

0 

21.  627 

0. 

.  510 

0. 

024 

-0. 

034 

0.  208 

-0. 

.  091 

0. 

042 

40. 

0 

25. 781 

-1. 

.  440 

-0. 

056 

-0. 

137 

0.  422 

-0, 

,  024 

0. 

061 

45. 

0 

28. 590 

0. 

.  878 

0. 

031 

-0. 

155 

0.  512 

-0. 

.  063 

0. 

053 

50. 

0 

29. 771 

2, 

.  273 

0. 

076 

-0. 

271 

0.  556 

0. 

.  094 

0. 

063 

55. 

0 

31.  914 

-1. 

.  837 

-0. 

058 

-0. 

475 

0.  627 

-0. 

.  061 

-0. 

031 

60. 

0 

32. 395 

0 

.  449 

0. 

014 

-0. 

644 

0.  774 

0. 

.  331 

0. 

085 

65. 

0 

31.  817 

-0. 

.  809 

-0. 

025 

-0. 

742 

0.  718 

0. 

.  066 

0. 

085 

70. 

0 

30.  741 

-3. 

.  424 

-0. 

111 

-0. 

796 

0.  492 

-0. 

.  273 

0. 

076 

80. 

0 

30. 122 

-0. 

.  621 

-0. 

021 

-0. 

939 

0.  061 

0 

.  023 

0. 

062 

90. 

0 

28. 647 

-0, 

.  878 

-0. 

031 

-0. 

169 

-1.  092 

0 

.  055 

0. 

035 

0. 

0 

0.  233 

-0. 

.  877 

-3. 

765 

-0. 

014 

-0. 003 

0. 

.  059 

0. 

006 

RUNFILE> 

AOA 

W0A402 

Cn 

Cy 

Cy/Cn 

Cax 

Cpt 

Crl 

Cyw 

***** 

****** 

i WnVvWnV 

****** 

****** 

****** 

0.  0 

-0. 012 

0.  004 

-0. 307 

-0. 003 

0.  000 

0.  001 

0.  001 

0.  0 

-0. 216 

-0. 331 

1.  531 

0.  228 

-0. 006 

0.  041 

-0. 006 

10.  0 

2.  666 

-0. 505 

-0. 190 

0.  219 

-0. 025 

0.  025 

-0. 008 

20.  0 

8.  329 

-0. 453 

-0. 054 

0.  198 

-0. 079 

-0. 017 

-0. 002 

30.  0 

16. 756 

-0. 653 

-0. 039 

0.  071 

0.  032 

-0. 079 

0.  029 

35.  0 

21.  215 

-0. 322 

-0. 015 

-0. 089 

0.  210 

-0. 012 

0.  002 

40.  0 

26. 011 

-2. 282 

-0. 088 

-0.  216 

0.  421 

0.  115 

-0.  011 

45.  0 

28.  622 

-1.  983 

-0. 069 

-0. 245 

0.  521 

0.  016 

-0. 007 

50.  0 

29. 554 

-3. 638 

-0.  123 

-0. 265 

0.  556 

-0. 133 

-0. 020 

55.  0 

31.  913 

-2. 071 

-0. 065 

-0. 554 

0.  523 

-0. 652 

0.  0t?4 

60.  0 

31.  676 

-2.  313 

-0. 073 

-0. 690 

0.  744 

-0.  773 

0.  072 

65.  0 

31.  362 

-3.  265 

-0. 104 

-0. 811 

0.  702 

-0.  517 

0.  070 

70.  0 

30.  758 

-1. 723 

-0. 056 

-0. 894 

0.  484 

-0.  709 

0.  062 

80.  0 

29.  772 

-1.  142 

-0. 038 

-1. 037 

0.  049 

-0. 017 

0.  057 

90.  0 

28.  623 

-0. 892 

-0. 031 

-0. 252 

-1. 081 

0.  056 

-0. 054 

0.  0 

-0.  016 

-0.  988 

60. 145 

-0. 104 

0.  007 

0.  049 

0.  005 

RUNFILE>  W0A502 


AOA 

Cn 

cy 

Cy/Cn 

Cax 

Cpt 

Crl 

Cyw 

***** 

•Mcirttek 

iHrtWnWr 

****** 

VrsWnWr* 

•irtrtriricis 

****** 

**- 

0.  0 

-0. 

077 

-0. 048 

0. 

625 

-0. 

009 

0.  001 

-0. 

002 

-0. 

001 

0.  0 

-0. 

063 

-0. 352 

5. 

622 

0. 

251 

0.  009 

0. 

052 

0. 

006 

10.  0 

2. 

839 

-0. 423 

-0. 

149 

0. 

213 

-0. 027 

0. 

020 

-0. 

006 

20.  0 

8. 

690 

-0. 523 

-0. 

060 

0. 

214 

-0. 082 

-0. 

002 

-0. 

019 

30.  0 

16. 

959 

0.  641 

0. 

038 

0. 

079 

0.  027 

-0. 

111 

-0. 

036 

35.  0 

21. 

235 

-0. 390 

-0. 

018 

-0. 

132 

0.  238 

-0. 

117 

0. 

023 

40.  0 

25. 

089 

1.  487 

0. 

059 

-0. 

269 

0.  405 

-0. 

132 

-0. 

024 

45.  0 

28. 

142 

3.  148 

0. 

112 

-0. 

222 

0.  498 

-0. 

09  7 

0. 

016 

50.  0 

29. 

573 

2.  443 

0. 

083 

-0. 

392 

0.  527 

0. 

118 

-0. 

040 

55.  0 

30. 

916 

-0. 356 

-0. 

012 

-0. 

557 

0.  625 

0. 

104 

-0. 

026 

60.  0 

31. 

881 

3.  536 

0. 

111 

-0. 

762 

0.  749 

0. 

484 

-0. 

024 

65.  0 

31. 

249 

2.  376 

0. 

076 

-0. 

856 

0.  707 

0. 

566 

-0. 

021 

70.  0 

30. 

491 

1.  249 

0. 

041 

-0. 

946 

0.  485 

0. 

782 

-0. 

020 

80.  0 

29. 

884 

-0. 970 

-0. 

032 

-1. 

065 

0.  084 

0. 

001 

-0. 

038 

90.  0 

28. 

417 

-1.  067 

-0. 

038 

-0. 

274 

-1.  093 

0. 

050 

-0. 

050 

0.  0 

0. 

092 

-0. 845 

-9. 

152 

-0. 

140 

0.  012 

0. 

013 

0. 

006 

RUNFILE> 

AOA 

W0A602 

Cn 

cy 

Cy/Cn 

Cax 

Cpt 

Crl 

Cyw 

jric'irfrie 

****** 

****** 

****** 

****** 

****** 

****** 

****** 

0.  0 

0.  052 

0.  005 

0.  093 

-0. 001 

-0. 002 

0.  000 

O.  002 

0.  0 

-0. 105 

-0. 344 

3.  280 

0.  249 

0.  007 

0.  052 

0.  004 

10.  0 

2.  953 

-0. 327 

-0. Ill 

0.  234 

-0. 035 

0.  024 

-0. 001 

20.  0 

8.  965 

-0. 323 

-0. 036 

0.  215 

-0. 079 

-0. 023 

-0. 022 

30.  0 

17. 203 

-0. 820 

-0. 048 

0.  068 

0.  027 

*0. 018 

0.  013 

35.  0 

21.  505 

0.  154 

0.  007 

-0. 078 

0.  216 

-0. 114 

-0. 043 

40.  0 

25. 362 

2.  513 

0.  099 

-0. 200 

0.427 

-0. 210 

0.  023 

45.  0 

28. 067 

3.  503 

0.  125 

-0. 192 

0.  462 

-0. 095 

0.  007 

50.  0 

29. 227 

4.  441 

0.  152 

-0. 359 

0.  490 

0.  252 

-0. 031 

55.  0 

31.  312 

1.  949 

0.  062 

-0. 537 

0.  490 

1.  006 

0.  001 

60.  0 

31. 776 

1.  238 

0.  039 

-0. 636 

0.  775 

0.  911 

-0.  003 

65.  0 

30. 524 

1.  946 

0.  064 

-0.  736 

0.  640 

0.  869 

0.  000 

70.  0 

30. 770 

1.  245 

0.  040 

-0. 898 

0.  532 

0.  857 

-0.  014 

80.  0 

30. 375 

-0. 853 

-0. 028 

-1.  017 

0.  063 

0.  049 

-0. 036 

90.  0 

28.  691 

-0. 915 

-0. 032 

-0. 269 

-1. 065 

0.  075 

-0.  051 

0.  0 

0.  335 

-0. 937 

-2.  793 

-0. 103 

0.  015 

0.  042 

0.  009 

RUNFILE>  W0A702 


AOA 

Cn 

Cy 

Cy/Cn 

Cax 

Cpt 

1 

3rl 

Cyw 

•frfcfrtrk 

*** 

*** 

**- 

ih'rfn't 

VwV’ 

*** 

:*** 

****** 

****** 

**-• 

.-*** 

0.  0 

-0. 

020 

0. 

009 

-0. 

447 

-0. 

002 

0.  002 

-0. 

.  001 

0. 

001 

O.  0 

-0. 

056 

-0, 

.  379 

6. 

753 

0. 

249 

0.  007 

0. 

052 

0. 

000 

10.  0 

2. 

897 

-0. 

,  384 

-0. 

133 

0. 

227 

-0. 031 

0. 

.  017 

-0. 

007 

20.  0 

8. 

870 

-0. 

.  290 

-0. 

033 

0. 

217 

-0. 088 

0. 

005 

-0. 

029 

30.  0 

16. 

961 

-0. 

,  037 

-0. 

002 

0. 

094 

0.  053 

-0. 

,  058 

-0. 

033 

35.  0 

21. 

542 

-0, 

.  499 

-0. 

023 

-0. 

094 

0.  251 

-0. 

.  025 

-0. 

013 

40.  0 

25. 

526 

-1, 

.  963 

-0. 

077 

-0. 

261 

0.  409 

0. 

,  040 

-0. 

013 

45.  0 

28. 

431 

-3 

.  326 

-0. 

117 

-0. 

257 

0.516 

0. 

.  075 

0. 

021 

50.  0 

29. 

908 

-3, 

.  194 

-0. 

107 

-0. 

331 

0.  564 

-0. 

,  050 

-0. 

Oil 

55.  0 

31. 

696 

'  0. 

.  845 

0. 

027 

-0. 

514 

0.  624 

-0. 

.  108 

-0. 

034 

60.  0 

32. 

462 

-1 

.  204 

-0. 

037 

-0. 

697 

0.  762 

0, 

,  010 

-0. 

028 

65.  0 

31. 

970 

-1. 

.  763 

-0. 

055 

-0. 

789 

0.  738 

-0, 

.  187 

-0. 

036 

70.  0 

30. 

668 

-3 

.  114 

-0. 

102 

-0. 

878 

0.  490 

-0. 

.  771 

-0. 

049 

80.  0 

30. 

570 

-0, 

.  787 

-0. 

026 

-V. 

996 

0.  078 

0 

.  002 

-0. 

037 

90.  0 

29. 

139 

-0. 

.  612 

-0. 

021 

-0. 

225 

-1. 093 

0, 

.  035 

-0. 

046 

0.  0 

0. 

099 

-0 

.  972 

-9. 

769 

-0. 

090 

0.  014 

0. 

.  040 

0. 

005 

RUNFILE> 

AOA 

W0A802 

Cn 

Cy 

Cy/Cn 

Cax 

Cpt 

Crl 

Cyw 

***** 

****** 

****** 

i Ir***** 

****** 

****** 

****** 

****** 

0.  0 

-0. 026 

-0. 014 

0.  526 

-0. 007 

-0. 001 

0.  000 

0.  000 

0.  0 

-0. 186 

-0.  441 

2.  366 

0.  240 

0.  005 

0.  055 

0.  000 

10.  0 

2.  767 

-0. 461 

-0. 167 

0.  224 

-0. 033 

0.  023 

-0. 008 

20.  0 

8.  803 

0.  018 

0.  002 

0.  216 

-0. 092 

-0. 031 

-0. 025 

30.  0 

17. 013 

0.  474 

0.  028 

0.  091 

0.  022 

-0. 123 

-0.  041 

35.  0 

21. 567 

-1.  157 

-0. 054 

-0. 094 

0.  218 

-0. 030 

-0. 030 

40.  0 

25. 496 

-2.  157 

-0. 085 

-0. 195 

0.  413 

0.  076 

0.  051 

45.  0 

28. 330 

-3.  053 

-0. 108 

-0. 221 

0.  501 

0.  025 

0.  030 

50.  0 

29. 706 

-5. 590 

-0. 188 

-0. 315 

0.  511 

-0. 245 

-0.  019 

55.  0 

31.  620 

-3. 679 

-0. 116 

-0. 546 

0.  548 

-0.7  64 

-0.  064 

60.  0 

31.  772 

-3. 244 

-0.  102 

-0. 632 

0.  771 

-0. 956 

-0. 069 

65.  0 

30. 272 

-3.  491 

-0. 115 

-0. 723 

0.  610 

-0.  841 

-0. 068 

70.  0 

30. 949 

-2.  403 

-0. 078 

-0. 888 

0.  514 

-0.  791 

-0. 058 

80.  0 

30. 145 

-0. 597 

-0. 020 

-1.  019 

0.  075 

-0. 009 

-0.  041 

90.  0 

28. 951 

-0. 788 

-0. 027 

-0. 221 

-1.  091 

0.  041 

-0. 055 

0.  0 

0.  106 

-0. 818 

-7.  683 

-0. 075 

0.  006 

0.  042 

0.  006 

84 


RUNFILE>  E0A825 


AOA 

Cn 

Cy 

Cy/Cn 

Cax 

Cpt 

Crl 

Cyw 

iWrtWrff 

****** 

****** 

•iririririrk 

****** 

****** 

****** 

0.  0 

-0. 036 

0.  017 

-0. 474 

-0.  007 

0.  000 

0.  003 

-0. 002 

0.  0 

0.  055 

-0.  006 

-1.  728 

0.  269 

-0. 002 

0.  010 

0.  000 

25.  0 

11.  701 

0.  170 

0.  015 

0.  241 

-0. 065 

-0. 045 

-0. 007 

25.  0 

0.  081 

-0. 014 

-0. 171 

0.  006 

0.  001 

0.  001 

0.  001 

25.  0 

11.  756 

0.  013 

0.  001 

0.  236 

-0. 048 

-0. 008 

0.  011 

25.  0 

0.  105 

-0. 030 

-0. 286 

-0. 005 

0.  0u2 

0.  002 

-0. 001 

25.  0 

11.  866 

-0. 352 

-0. 030 

0.  234 

-0. 052 

0.  003 

-0. 007 

25.  0 

0.  213 

0.  024 

0.  114 

0.  005 

0.  001 

0.  002 

0.  001 

25.  0 

12. 130 

0.  607 

0.  050 

0.  245 

-C. 059 

-0. 029 

"  vy  .  -L«+ 

0.  0 

0.  227 

-0.  112 

-0. 491 

0.  263 

-0. 004 

0.  005 

-0. 003 

0.  0 

0.  233 

0.  001 

-0. 003 

-0.  012 

-0. 002 

-0. 004 

-0.  001 

RUNFILE>  E0A845 


AOA 

Cn 

cy 

Cy/Cn 

Cax 

Cpt 

Crl 

Cyw 

***** 

****** 

****** 

****** 

****** 

****** 

****** 

****** 

0.  0 

-0. 070 

-0. 038 

0.  545 

0.  001 

0.  000 

0.  002 

-0. 002 

0.  0 

-0. 031 

-0. 064 

2.  039 

0.  289 

-0. 003 

0.  003 

-0.  005 

45.  0 

23. 703 

-5. 013 

-0. 212 

0.  062 

0.  417 

0.  055 

0.  031 

45.  0 

0.  010 

-0. 029 

-3.  049 

0.  074 

-0. 001 

-0.  006 

0.  008 

45.  0 

23.  787 

-3.  719 

-0.  156 

0.  096 

0.  399 

-0.  026 

0.  048 

45.  0 

0.  077 

-0. 020 

-0.  261 

0.  091 

-0. 005 

'0. 008 

0.  019 

45.  0 

23. 583 

-3. 904 

-0. 166 

-0. 068 

0.  365 

0.  001 

0.  055 

45.  0 

0.  074 

-0. 083 

-1. 113 

0.  091 

-0. 007 

-0. 011 

0.  016 

45.  0 

23. 971 

-4. 599 

-0.  192 

-0. 040 

0.  389 

-0. 019 

0.  060 

0.  0 

0.  179 

-0. 056 

-0. 314 

0.  355 

-0. 010 

|»Hj  wH 

0.  013 

0.  0 

0.  225 

0.  092 

0.  411 

0.  057 

-0. 011 

0.  011 

RUNFILE>  E0A870 


AOA 

Cn 

cy 

Cy/Cn 

Cax 

Cpt 

Crl 

Cyw 

***** 

****** 

****** 

ihWnhWf 

****** 

****** 

****** 

UV..W. 

0.  0 

-0. 017 

-0. 015 

0.  911 

-0.  009 

-0.  004 

0.  000 

0.  001 

0.  0 

-0. 068 

-0. 252 

3.  711 

0.  255 

-0. 002 

0.  014 

-0. 005 

70.  0 

30. 017 

-2. 283 

-0. 076 

-0.  760 

0.  522 

-0. 613 

-0.  060 

70.  0 

0.  179 

0.  001 

0.  008 

0.  029 

0.  012 

0.  004 

-0. 004 

70.  0 

29. 783 

-1.  754 

-0. 059 

-0. 772 

0.  485 

-0. 669 

-0. 060 

70.  0 

0.  314 

0.  010 

0.  031 

0.  034 

0.  013 

0.  004 

-0. 007 

70.  0 

30. 018 

-0. 497 

-0.  017 

-0. 749 

0.  491 

-0. 652 

-0. 064 

70.  0 

0.  271 

-0. 091 

-0. 335 

0.  024 

0.  012 

0.  000 

-0.  012 

70.  0 

30. 248 

-2.  926 

-0. 097 

-0. 742 

0.  493 

-0. 681 

-0. 062 

0.  0 

0.414 

-0.  096 

-0. 232 

0.  259 

0.  001 

0.  012 

-0. 014 

0.  0 

0.  390 

-0. 085 

-0. 217 

-0. 014 

0.  004 

0.  003 

-0. 013 

RUNFILE> 

AOA 

E0A890 

Cn 

cy 

Cy/Cn 

Cax 

Cpt 

Crl 

Cyw 

***** 

****** 

****** 

****** 

****** 

****** 

'trfcJn’n’c’k 

****** 

0.  0 

0.  035 

0.  022 

0.  620 

0.  009 

-0. 003 

0.  001 

0.  000 

0.  0 

0.  062 

-0. 210 

-3. 373 

0.  280 

-0. 003 

0.  018 

-0. 002 

90.  0 

27.  764 

-0. 039 

-0. 001 

-0. 065 

-1.  066 

0.  000 

-0. 053 

90.  0 

0.  167 

0.  038 

0.  229 

0.  062 

-0. 001 

0.  003 

-0. 001 

90.  0 

27. 962 

0.  120 

0.  004 

-0. 069 

-1. 074 

0.  006 

-0. 051 

90.  0 

0.  222 

0.  055 

0.  246 

0.  065 

0.  000 

0.  003 

-0. 001 

90.  0 

28. 331 

0.  139 

0.  005 

-0. 028 

-1.  090 

0.  015 

-0. 054 

90.  0 

0.  351 

0.  051 

0.  145 

0.  083 

-0. 002 

0.  006 

0.  000 

90.  0 

28. 076 

0.  045 

0.  002 

-0. 011 

<1.  096 

0.  007 

-0. 055 

0.  0 

0.  378 

-0. 050 

-0. 132 

0.  326 

-0. 011 

0.  012 

-0. 011 

0.  0 

0.  369 

0.  053 

0.  144 

0.  049 

-0.  006 

0.  004 

-0.  011 

S6 


RUNFILE>  R0A803 


AOA 

Cn 

cy 

Cy/Cn 

Cax 

Cpt 

Cr  1 

Cyw 

***** 

****** 

****** 

****** 

****** 

id rMcirk 

****** 

-5.  0 

0.  019 

0.  009 

0.  500 

0.  002 

0.  001 

0.  000 

0.  000 

-5.  0 

-1.  571 

-0. 022 

0.  014 

0.  262 

0.  074 

0.  005 

-0. 002 

-4.  0 

-1.  133 

-0.  010 

0.  008 

0.  297 

0.  083 

0.  016 

0.  002 

-3.  0 

-0. 863 

-0. 081 

0.  091 

0.  297 

0.  078 

0.  014 

0.  000 

-2.  0 

-0. 712 

-0. 156 

0.  219 

0.  289 

0.  052 

0.  012 

-0. 002 

-1.  0 

-0. 485 

-0. 137 

0.  282 

0.  289 

0.  007 

0.  013 

-0. 001 

0.  0 

-0. 097 

-0. 095 

0.  965 

0.  274 

0.  001 

0.  012 

-0. 002 

1.  0 

0.  062 

-0. 145 

-2. 344 

0.  260 

-0. 012 

0.  009 

0.  000 

2.  0 

0.  362 

-0. 176 

-0. 485 

0.  260 

-0. 034 

0.  013 

0.  001 

3.  0 

0.  685 

-0. 145 

-0. 212 

0.  256 

-0. 090 

0.  014 

-0. 002 

4.  0 

1.  050 

-0. 062 

-0. 059 

0.  250 

-0. 094 

0.  009 

-0. 002 

5.  0 

1.  391 

-0. 011 

-0. 008 

0.  253 

-0. 101 

0.  005 

-0. 005 

6.  0 

1.  734 

-0. uu9 

-0.  005 

0.  258 

-0. 099 

-0. 005 

0.  001 

7.  0 

2.  168 

0.  097 

0.  045 

0.  252 

-0.  077 

0.  004 

-0. 004 

8.  0 

2.  318 

-0. 011 

-0. 005 

0.  242 

-0. 067 

-0. 006 

-0. 004 

9.  0 

2.  629 

-0. 103 

-0. 039 

0.  252 

-0. 057 

0.  007 

-0. 027 

10.  0 

3.  110 

-0. 143 

-0. 046 

0.  257 

-0. 048 

0.  007 

-0. 003 

11.  0 

3.  457 

0.  008 

0.  002 

0.  263 

-0. 040 

-0. 007 

0.  004 

12.  0 

3.  824 

-0. 232 

-0.  Oftl 

0.  277 

-0. 037 

0.  002 

-0. 006 

13.  0 

4.  556 

-0. 174 

-0.  038 

0.  286 

-0. 035 

0.  003 

0.  008 

14.  0 

5.  034 

-0. 102 

-0. 020 

0.  290 

-0. 045 

0.  008 

-0. 081 

15.  0 

5.  781 

0.  026 

0.  005 

0.  309 

-0. 062 

-0. 011 

0.  012 

16.  0 

6.  377 

0.  017 

0.  003 

0.  309 

-0. 077 

-0. 004 

0.  013 

17.  0 

6.  965 

-0. 064 

-0. 009 

0.  373 

-0. 096 

0.  004 

-0. 009 

18.  0 

7.  645 

0.  000 

0.  000 

0.  362 

-0.  Ill 

-0.  001 

0.  010 

19.  0 

8.  399 

0.  249 

0.  030 

0.  332 

-0. 110 

-0. 005 

0.  002 

20.  0 

9.  357 

-0. 095 

-0. 010 

0.  344 

-0. 101 

0.  007 

0.  010 

21.  0 

10. 089 

0.  116 

0.  011 

0.  339 

-0. 103 

0.  014 

0.  017 

22.  0 

10. 830 

-0. 007 

0.  001 

0.  330 

-0. 096 

-0.  013 

-0. 016 

23.  0 

11.  729 

0.  233 

0.  020 

0.  327 

-0. 093 

-0. 017 

-0. 010 

24.  0 

12. 480 

-0.  421 

-0. 034 

0.  324 

-0. 087 

-0. 021 

0.  021 

25.  0 

13.  198 

-0.  562 

-0. 043 

0.  326 

-0. 079 

0.  000 

0.  019 

26.  0 

13.  879 

0.  Ill 

0.  008 

0.  331 

-0. 084 

-0. 020 

0.  019 

27.  0 

14.  984 

0.  690 

0.  046 

0.  310 

-0. 050 

-0. 031 

0.  020 

28.  0 

15.  652 

0.  171 

0.  011 

0.  305 

-0. 031 

-0. 048 

0.  031 

29.  0 

16.  661 

0.  540 

0.  032 

0.  269 

-0. 005 

-0. 011 

0.  02 0 

30.  0 

17.  315 

0.  623 

0.  036 

0.  268 

-0.  003 

-0.  010 

0.  016 

31.  0 

18.  474 

0.  428 

0.  023 

0.  249 

0.  046 

-0. 052 

0.  003 

32.  0 

19.  390 

1.  139 

0.  059 

0.  240 

0.  061 

-0. 153 

0.  042 

33.  0 

20. 121 

-0. 332 

-0. 016 

0.  203 

0.  124 

-0. 036 

-0. 022 

34.  0 

20.  609 

-1. 704 

-0. 083 

0.  163 

0.  133 

0.  076 

0.  003 

35.  0 

21.  787 

-2. 381 

-0. 109 

0.  155 

0.  211 

0.  142 

0.  020 

36.  0 

22.  638 

-1.  079 

-0. 048 

0.  130 

0.  264 

0.  055 

0.  021 

37.  0 

23. 468 

-2. 219 

-0.  095 

0.  163 

0.  299 

0.  117 

0.  026 

38.  0 

24. 326 

-2. 246 

-0. 092 

0.  130 

0.  317 

0.  162 

0.  021 

39.  0 

24. 990 

-2. 310 

-0. 092 

0.  136 

0.  351 

0.  112 

0.  003 

40.  0 

25. 772 

-3. 255 

-0. 126 

0.  101 

0.  378 

0.  099 

0.  039 

41.  0 

26. 707 

-2. 418 

-0.  091 

0.  131 

0.  405 

0.  099 

0.  073 

42.  0 

26. 933 

-1.  499 

-0. 056 

0.  066 

0.  413 

0.  109 

0.  030 

43.  0 

27.  541 

-3.  006 

-0. 109 

0.  Ill 

0.  446 

0.  065 

0.  034 

44.  0 

27.  894 

-3.  189 

-0. 114 

0.  087 

0.  466 

0.  103 

0.  026 

45.  0 

28. 377 

-3.  959 

-0. 140 

0.  114 

0.  519 

0.  045 

0.  033 

S9 

46. 

0 

28. 

521 

-4. 

553 

-0. 

160 

0. 

059 

0. 

495 

-0. 

046 

0. 

010 

47. 

0 

29. 

113 

-4. 

111 

0. 

141 

0. 

069 

0. 

469 

0. 

003 

0. 

017 

48. 

0 

29. 

123 

-4. 

188 

-0. 

144 

0. 

058 

0. 

493 

-0. 

137 

-0. 

015 

49. 

0 

29. 

985 

-2. 

544 

-0. 

085 

0. 

003 

0. 

477 

-0. 

171 

-0. 

024 

50. 

0 

30. 

283 

-5. 

268 

-0. 

174 

-0. 

044 

0. 

464 

-0. 

371 

0. 

072 

51. 

0 

30. 

627 

343 

-0. 

077 

-0. 

064 

0. 

446 

-0. 

391 

0. 

065 

52. 

0 

31. 

124 

-2. 

938 

-0. 

094 

-0. 

123 

0. 

442 

-0. 

551 

0. 

06 1 

53. 

0 

31. 

663 

-2. 

681 

-0. 

085 

-0. 

233 

0. 

385 

-0. 

793 

0. 

055 

54. 

0 

31. 

377 

-3. 

442 

-0. 

110 

-0. 

253 

0. 

492 

-0. 

821 

0. 

049 

55. 

0 

31. 

470 

-4. 

262 

-0. 

135 

-0. 

322 

0. 

528 

-0. 

806 

0. 

054 

56. 

0 

32. 

256 

-3. 

603 

-0. 

112 

-0. 

331 

0. 

574 

-0. 

878 

0. 

058 

57. 

0 

32. 

078 

-2. 

561 

-0. 

080 

-0. 

367 

0. 

622 

-0. 

938 

0. 

064 

58. 

0 

32. 

058 

-2. 

845 

-0. 

089 

-0. 

382 

0. 

689 

-0. 

823 

0. 

065 

59. 

0 

32. 

015 

-2. 

225 

-0. 

069 

-0. 

392 

0. 

694 

-0. 

622 

0. 

070 

60. 

0 

32. 

287 

-1. 

550 

-0. 

048 

-0. 

398 

0. 

735 

-0. 

541 

0. 

072 

61. 

0 

31. 

894 

-4. 

297 

-0. 

135 

-0. 

434 

0. 

735 

-0. 

521 

0. 

075 

62. 

0 

31. 

463 

-1. 

218 

-0. 

039 

-0. 

455 

0. 

725 

-0. 

517 

0. 

074 

63. 

0 

31. 

136 

-2. 

010 

-0. 

065 

-0. 

466 

0. 

730 

-0. 

526 

0. 

065 

64. 

0 

30. 

627 

-3. 

069 

-0. 

100 

-0. 

492 

0. 

660 

-0. 

633 

0. 

062 

65. 

0 

30. 

363 

-3. 

339 

-0. 

110 

-0. 

496 

0. 

602 

-0. 

664 

0. 

061 

66. 

0 

30. 

044 

-3. 

228 

-0. 

107 

-0. 

506 

0. 

509 

-0. 

712 

0. 

059 

67. 

0 

30. 

268 

_2. 

573 

-0. 

085 

-0. 

534 

0. 

495 

-G. 

760 

0. 

068 

68. 

0 

30. 

598 

_  o 

744 

-0. 

090 

-0. 

576 

0. 

496 

-0. 

781 

0. 

069 

69. 

0 

30. 

547 

-1. 

068 

-0. 

035 

-0. 

618 

0. 

463 

-0. 

632 

0. 

061 

70. 

0 

30. 

541 

_2. 

695 

-0. 

085 

-0. 

648 

0. 

463 

-0. 

661 

0. 

059 

71. 

0 

30. 

397 

-3. 

285 

-0. 

108 

-0. 

667 

0. 

447 

-0. 

383 

0. 

063 

72. 

0 

30. 

379 

0. 

461 
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0.  054 

0.  050 

1U2 


46. 

0 

22. 219 

-2. 

178 

47. 

0 

22. 238 

-3. 

791 

48. 

0 

22. 551 

-2. 

973 

49. 

0 

22.  791 

-2. 

327 

50. 

0 

23.  130 

-2. 

326 

51. 

0 

23.  516 

-1. 

489 

52. 

0 

23. 656 

-2. 

645 

53. 

0 

23. 790 

-2. 

102 

54. 

0 

23. 902 

-2. 

489 

55. 

0 

24. 324 

-1. 

600 

56. 

0 

24.  434 

-2. 

683 

57. 

0 

24. 299 

-2. 

028 

58. 

0 

24. 431 

-1. 

118 

59. 

0 

24. 213 

753 

60. 

0 

24.414 

-1. 

536 

61. 

0 

24.  485 

0. 

067 

62. 

0 

24. 131 

“  2, , 

193 

63. 

0 

23.  917 

0. 

499 

64. 

0 

23.  386 

-0. 

929 

65. 

0 

23. 285 

-1. 

575 

66. 

0 

23.  036 

-1. 

731 

67. 

0 

22.  921 

-1. 

506 

68. 

0 

23.  093 

-1. 

532 

69. 

0 

23. 362 

-1. 

432 

70. 

0 

23. 330 

-1. 

220 

71. 

0 

23.  460 

-1. 

441 

72. 

0 

23.  627 

.0 

4.  « 

593 

73. 

0 

23.463 

-0. 

035 

74. 

0 

23. 408 

-0. 

145 

75. 

0 

23.475 

0. 

236 

76. 

0 

23. 133 

0. 

240 

77. 

0 

23.  197 

0. 

181 

78. 

0 

22.  989 

0. 

035 

79. 

0 

23. 200 

0. 

263 

80. 

0 

22. 921 

0. 

131 

81. 

0 

22. 93b 

0. 

113 

82. 

0 

22. 485 

0. 

151 

83. 

0 

22. 283 

0. 

120 

84. 

0 

22. 138 

0. 

235 

85. 

0 

22.  152 

0. 

073 

86. 

0 

22. 078 

0. 

244 

87. 

0 

22.  413 

0. 

182 

88. 

0 

22. 472 

0. 

029 

89. 

0 

21.  989 

0. 

001 

90. 

0 

22. 391 

-0. 

040 

91. 

0 

22. 264 

0. 

159 

92. 

0 

22.  348 

0. 

084 

93. 

0 

22. 285 

0. 

066 

94. 

0 

22. 233 

-0. 

115 

95. 

0 

22. 216 

0. 

145 

95. 

0 

1.  531 

0. 

058 

-5. 

0 

1.466 

0. 

014 

0. 

0 

1.  544 

0. 

059 

-0.098  0.230 

-0.170  0.212 

-0.132  0.211 

-0. 102  0. 194 

-0.  101  0.  154 

-0.063  0.127 

-0. 112  0. 119 

-0.088  0.084 

-0.104  0.033 

-0. 066  0. 025 

-0.110  0.000 

-0.083  -0.004 

-0.046  -0.010 

-0.114  -0.012 

-0.063  -0.113 

0.003  -0.094 

-0.049  -G.  114 

0.021  -0.137 

-0.040  -0.137 

-0. 068  -0. 150 

-0.075  -0.151 

-0.066  -0.171 

-0.066  -0.187 

-0.061  -0.212 

-0.052  -0.223 

-0.061  -0.263 

-0.110  -0.290 

-0.001  -0.308 

-0.006  -0.339 

0.010  -0.347 

0.010  -0.343 

0. 008  -0. 339 

0.002  -0.324 

0.012  -0.324 

0.006  -0.317 

0.005  -0.269 

0.007  -0.238 

0.  005  -0. 194 

0.  011  -0.  149 

0.003  -0.084 

0.011  -0.018 

0.  008  0. 035 

0.001  0.132 

0.  000  0.  202 

-0.002  0.225 

0.007  0.224 

0.004  0.251 

0.003  0.248 

-0.  005  0.  240 

0.007  0.247 

0.038  0.311 

0.  009  0. 244 

0.039  0.257 
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0.  390 

0. 

077 

0. 

0.  348 

-0. 

014 

0. 

0.  330 

0. 

041 

0. 

0.  348 

-0. 

055 

0. 

0.  364 

-0. 

099 

0. 

0.  366 

-0. 

259 

-0. 

0.  350 

-0. 

324 

-0. 

0.  363 

-0. 

402 

-0. 

0.  322 

-0. 

603 

-0. 

0.  334 

-0. 

685 

-0. 

0.  384 

-0. 

653 

-0. 

0.  426 

-0. 

644 

-0. 

0.  478 

-0. 

576 

-0. 

0.  499 

-0. 

556 

-0. 

0.  503 

-0. 

421 

-0. 

0.  500 

-0. 

263 

-0. 

0.  497 

-0. 

329 

-0. 

0.  496 

-0. 

272 

-0. 

0.  471 

-0. 

354 

-0. 

0.  413 

-0. 

460 

-0. 

0.  366 

-0. 

554 

-0. 

0.  338 

-0. 

524 

-0. 

0.  304 

-0. 

503 

-0. 

0.  293 

-0. 

516 

-0. 

0.  286 

-0. 

466 

-0. 

0.  277 

-0. 

266 

-0. 

0.  263 

-0. 

256 

-0. 

0.  257 

-0. 

Oil 

-0. 

0.  227 

0. 

053 

0. 

0.  194 

0. 

043 

-0. 

0.  121 

0. 

045 

-0. 

0.  034 

( . 

036 

-0. 

0.  039 

0. 

Oil 

0. 

0.  011 

0. 

021 

-0. 

-0. 074 

0. 

018 

-0. 

-0. 138 

0. 

027 

-0. 

-0. 213 

0. 

029 

-0. 

-0. 257 

0. 

030 

-0. 

-0.  348 

0. 

030 

-0. 

-0. 452 

0. 

044 

-0. 

-0.  513 

0. 

029 

-0. 

-0. 623 

0. 

028 

-0. 

-0. 754 

0. 

049 

-0. 

-0.  809 

0. 

019 

-0. 

-0. 866 

0. 

036 

-0. 

-0. 954 

0. 

032 

-0. 

-1.  006 

0. 

037 

-0. 

-1.  072 

0. 

039 

-0. 

-1. 118 

0. 

038 

-0. 

-1.  178 

0. 

048 

-0. 

-0. 049 

0. 

037 

0. 

-0. 049 

0. 

037 

0. 

-0. 049 

0. 

038 

0. 

036 

025 

019 

Oil 

004  ■ 

005 

007 

026 

027 
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020  | 

020  ■ 

020 

017 

008 
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004 

CIO  Q 

017 
019 
016 
013 
013 
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009  • 

007  : 

003  i 
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000  j 
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RUNFILE>  R0C801 


AOA 

Cn 

cy 

Cy/Cn 

Cax 

Cpt 

Crl 

Cyw 

***** 

****** 

****** 

JcSfiricVcte 

****** 

****** 

0.  0 

-0. 008 

*0.  013 

1.  695 

-0. 003 

0.  001 

0.  000 

-0. 002 

0.  0 

0.  126 

-0. 342 

-2.  715 

0.  349 

-0. 004 

0.  063 

0.  000 

10.  0 

3.  616 

0.  052 

0.  014 

0.  271 

-0. 130 

-0. 023 

-0. 002 

20.  0 

7.  671 

-0.  204 

-0.  027 

0.  375 

0.  124 

-0. 002 

0.  006 

30.  0 

12. 631 

0.  200 

0.  016 

0.  347 

0.  394 

-0. 014 

0.  001 

35.  0 

15. 264 

-0. 486 

-0. 032 

0.  271 

0.  493 

-0. 045 

0.  006 

40.  0 

18. 592 

*0. 498 

-0. 027 

0.  161 

0.  416 

-0. 229 

-0. 031 

45.  0 

21.  397 

1.  037 

0.  048 

0.  078 

0.  451 

-0.  449 

-0. 012 

50.  0 

24. 332 

-0. 133 

-0. 005 

-0. 081 

0.  545 

-0. 385 

0.  015 

55.  0 

26. 783 

-2. 625 

-0. 098 

-0. 230 

0.  670 

-0.  297 

-0.  007 

60.  0 

29. 142 

0.  815 

0.  028 

-0. 372 

0.  867 

0.  044 

-0. 043 

65.  0 

30. 947 

-0. 818 

-0. 026 

-0. 571 

1.  058 

-0. 725 

0.  045 

70.  0 

31.  683 

-2.475 

-0. 078 

-0. 701 

1.  155 

-0. 947 

-0. 079 

80.  0 

31.  932 

0.  073 

0.  002 

-0. 780 

0.  800 

-0. 054 

0.  096 

90.  0 

30. 062 

0.  292 

0.  010 

-0. 241 

-1.  041 

-0. 028 

0.  122 
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RUNFIL£>  R0C802 


AOA 

Cn 

Cy 

Cy/qn 

****** 

Cax 

Cpt 

Crl 

Cyw 

***** 

****** 

****** 

****** 

Hn’nWri Wr 

-sVsWwWoV 

-5.  0 

0.  053 

0.  013 

0.  245 

-0. 006 

0.  001 

0.  000 

0.  001 

-5.  0 

-1.  357 

-0. 056 

0.  041 

0.  250 

0.  109 

0.  002 

-0.  004 

-4.  0 

-0.  939 

-0. 063 

0.  067 

0.  250 

0.  087 

0.  004 

0.  001 

-3.  0 

-0. 637 

-0.  131 

0.  205 

0.  247 

0.  073 

0.  012 

0.  001 

-2.  0 

-0.  442 

-0.  161 

0.  364 

0.  220 

0.  046 

0.  008 

-0. 004 

-1.  0 

-0. 021 

-0. 080 

3.  773 

0.  220 

0.  023 

0.  C07 

-0. 019 

0.  0 

0.  110 

-0. 064 

-0. 583 

0.  212 

0.  008 

0.  008 

-0. 006 

1.  0 

0.  498 

-0. 088 

-0. 177 

0.  209 

-0. 014 

0.  010 

0.  000 

2.  0 

0.  820 

0.  047 

0.  057 

0.  220 

-0. 032 

0.  003 

-0. 003 

3.  0 

1.  082 

0.  025 

0.  023 

0.  215 

-0. 059 

0.  004 

-0. 005 

0.  0 

1.  279 

*0. 149 

-0. 117 

0.  213 

-0. 077 

0.  001 

-0. 006 

4.  0 

1.  292 

-0. 156 

-0. 121 

0.  213 

-0. 082 

0.  005 

-0.  007 

5.  0 

1.  610 

0.  063 

0.  039 

0.  227 

-0.  101 

-0. 009 

-0. 010 

6.  0 

1.  991 

-0.  051 

-0.  026 

0.  234 

-0. 114 

-0. 006 

-0. 005 

7.  0 

2.  328 

0.  087 

0.  037 

0.  235 

-0.  121 

-0.  001 

-0. 009 

8.  0 

2.  622 

-0. 212 

-0. 081 

0.  225 

-0. 125 

0.  001 

-0. 007 

9.  0 

2.  973 

-0. 036 

-0. 012 

0.  246 

-0. 124 

-0. 001 

-0. 004 

10.  0 

3.  183 

-0. 382 

-0. 120 

0.  224 

-0. 113 

0.  006 

-0. 013 

11.  0 

3.  439 

-0. 035 

-0. 010 

0.  241 

-0. 105 

0.  007 

-0. 015 

12.  0 

3.  984 

-0. 056 

-0. 014 

0.  273 

-0. 097 

>0. 004 

-0. 005 

13.  0 

4.  334 

-0. 052 

-0. 012 

0.  269 

-0. 085 

-0.  001 

-0. 006 

14.  0 

4.  731 

-0. 060 

-0. 013 

0.  304 

-0. 064 

0.  016 

-0. 006 

15.  0 

4.  880 

0.  034 

0.  007 

0.  304 

-0.  044 

0.  002 

-0. 011 

16.  0 

5.  266 

-0.  037 

-0. 007 

0.  321 

-0. 026 

-0. 003 

-0. 010 

17.  0 

5.  592 

0.  001 

0.  000 

0.  346 

0.  005 

0.  002 

-0. 006 

18.  0 

5.  918 

0.  076 

0.  013 

0.  355 

0.  023 

0.  001 

-0. 010 

19.  0 

6.  395 

0.  227 

0.  035 

0.  367 

0.  052 

0.  004 

-0. 003 

20.  0 

6.  317 

0.  259 

0.  041 

0.  364 

0.  063 

-0. 009 

-0. 015 

21.  0 

6.  735 

0.  038 

0.  006 

0.  389 

0.  095 

-0. 025 

-0. 013 

22.  0 

7.  303 

0.  194 

0.  027 

0.  417 

0.  119 

-0. 013 

-0. 003 

23.  0 

7.  617 

-0. 049 

-0. 006 

0.  435 

0.  135 

-0. 034 

-0.  016 

24.  0 

7.  948 

-0. 517 

-0. 065 

0.  444 

0.  144 

-0. 005 

-0. 002 

25.  0 

8.  35S 

0.  295 

0.  035 

0.  460 

0.  154 

-0. 057 

-0, 010 

26.  0 

8.  409 

0.  236 

0.  028 

0.  435 

0.  181 

-0.  024 

-0. 003 

27.  0 

9.  041 

0.  081 

0.  009 

0.  462 

0.  203 

-0. 047 

0.  000 

28.  0 

9.  473 

0.  512 

0.  054 

0.  465 

0.  229 

-0. 058 

-0. 008 

29.  0 

9.  701 

0.  125 

0.  013 

0.  452 

0.  245 

-0. 021 

-0. 016 

30.  0 

10.  470 

0.  873 

0.  083 

0.  454 

0.  276 

-0. 073 

-0.  014 

31.  0 

10.  410 

0.  012 

0.  001 

0.  415 

0.  307 

-0. 016 

-0. 017 

32.  0 

10. 906 

0.  334 

0.  031 

0.  399 

0.  336 

-0. 018 

-0. 025 

33.  0 

11.  086 

0.  598 

0.  054 

0.  384 

0.  352 

-0. 066 

-0. 024 

34.  0 

11.  180 

0.  640 

0.  057 

0.  456 

0.  311 

0.  006 

-0. 083 

35.  0 

12. 185 

0.  603 

0.  049 

0.  370 

0.  373 

-0. 047 

-0. 019 

36.  0 

12.  780 

0.  302 

0.  024 

0.  366 

0.  374 

-0. 126 

-0. 024 

37.  0 

12.  671 

0.  687 

0.  054 

0.  433 

0.  315 

-0. 054 

-0. 099 

38.  0 

13. 622 

1.  026 

0.  075 

0.  334 

0.  357 

-0. 123 

-0. 047 

39.  0 

14. 113 

2.  940 

0.  208 

0.  312 

0.  356 

-0. 164 

-0. 061 

40.  0 

14. 465 

-0. 212 

-0.  015 

0.  302 

0.  357 

-0. 088 

-0. 028 

41.0 

15.  264 

-0.  182 

-0. 012 

0.  312 

0.  336 

-0. 191 

-0. 030 

42.  0 

15. 868 

-0.  130 

C.  008 

C.  301 

0.  330 

-0. 081 

-0. 047 

43.  0 

15.  947 

0.  222 

0.  014 

0.  263 

0.  347 

-0. 235 

-0. 070 

44.  0 

16. 521 

1.  160 

0.  070 

0.  261 

0.  342 

-0. 249 

-0. 0h3 

105 
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RUNFILE>  R1C801 


AOA 

Cn 

Cy 

Cy/Cn 

Cax 

Cpt 

Crl 

Cyw 

***** 

y ,  j.  s. 

****** 

****** 

****** 

iSr?hwwV» 

****** 

-5.  0 

-0. 056 

-0. 052 

0.  938 

-0. 009 

0.  002 

0.  000 

-0. 001 

-5.  0 

-1.  182 

-0. 127 

0.  108 

0.  287 

0.  050 

-0. 005 

0.  002 

-4.  0 

-1.  042 

-0. 197 

0.  189 

0.  271 

0.  040 

-0.  004 

-0. 001 

-3.  0 

-0. 643 

-0. 184 

0.  286 

0.  261 

0.  032 

-0. 002 

-0. 002 

-2.  0 

-0.  276 

-0. 159 

0.  576 

0.  259 

0.  016 

-0. 003 

0.  003 

-1.  0 

0.  056 

-0. 132 

-2. 354 

0.  258 

0.  007 

-0. 002 

0.  001 

0.  0 

0.  359 

-0. 186 

-0.  519 

0.  252 

-0. 003 

-0. 003 

0.  001 

1.  0 

0.  945 

-0. 011 

-0. 011 

0.  293 

-0. 015 

-0. 002 

0.  009 

2.  0 

1.  186 

-0. 015 

-0. 013 

0.  300 

-0. 025 

0.  002 

0.  004 

3.  0 

1.  558 

■0. 032 

0.  020 

0.  333 

-0. 036 

0.  001 

0.  009 

4.  0 

1.  832 

0.  058 

0.  032 

0.  340 

-0. 052 

0.  003 

0.  009 

5.  0 

2.  132 

0.  127 

0.  060 

0.  363 

-0. 062 

0.  001 

0.  011 

6.  0 

2.  542 

0.  107 

0.  042 

0.  365 

-0. 074 

0.  002 

0.  011 

7.  0 

2.  844 

0.  129 

0.  045 

0.  390 

-0. 083 

-0. 002 

0.  013 

8.  0 

3.  266 

0.  093 

0.  028 

0.  400 

-0. 090 

0.  004 

0.  012 

9.  0 

3.  552 

0.  169 

0.  048 

0.  422 

-0. 089 

-0. 003 

0.  014 

10.  0 

3.  751 

0.  120 

0.  032 

0.  422 

-0. 092 

-0. 002 

0.  013 

11.  0 

4.  065 

0.  144 

0.  035 

0.  422 

-0. 086 

-0. 002 

0.  012 

12.  0 

4.  503 

0.  240 

0.  053 

0.  455 

-0. 082 

-0. 006 

0.  014 

13.  0 

4.  768 

0.  225 

0.  047 

0.  454 

-0. 071 

-0. 001 

0.  012 

14.  0 

5.  046 

0.  283 

0.  056 

0.  468 

-0. 062 

-0. 001 

0.  012 

15.  0 

5.  535 

0.  364 

0.  066 

0.  515 

-0. 058 

0.  000 

0.  015 

16.  0 

5.  869 

0.  269 

0.  046 

0.  507 

-0. 046 

-0. 001 

0.  014 

17.  0 

6.  311 

0.  448 

0.  071 

0.  542 

-0. 043 

0.  008 

0.  017 

18.  0 

6.  416 

0.  448 

0.  070 

0.  524 

-0. 033 

0.  001 

0.  012 

19.  0 

6.  811 

0.  288 

0.  042 

0.  536 

-0. 028 

0.  011 

0.  011 

20.  0 

7.  202 

0.  547 

0.  076 

0.  575 

-0. 022 

0.  009 

0.  018 

21.  0 

7.  552 

0.  608 

0.  080 

0.  567 

-0. 012 

-0. 009 

0.  015 

22.  0 

8.  076 

0.  638 

0.  079 

0.  608 

-0. 011 

-0. 001 

0.  018 

23.  0 

8.  460 

0.  601 

0.  071 

0.  631 

0.  000 
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0.  534 

0.  020 

-0.  001 

32.  0 

3.  757 

-0. 232 

-0. 062 

0.  306 

0.  551 

0.  011 

-0. 004 

33.  0 

3.  942 

-0. 156 

-0. 040 

0.  315 

0.  576 

0.  016 

-0. 002 

34.  0 

4.  145 

-0. 362 

-0. 087 

0.  319 

0.  602 

-0. 029 

-0. 003 

35.  0 

4.  335 

-0. 644 

-0. 148 

0.  317 

0.  619 

-0. 025 

-0. 001 

36.  0 

4.  358 

-0. 648 

-0. 149 

0.  317 

0.  637 

-0. 050 

-0. 004 

37.  0 

4.  670 

-0. 589 

-0. 126 

0.  329 

0.  671 

-0. 016 

-0. 002 

38.  0 

4.  863 

-0. 875 

-0. 180 

0.  323 

0.  690 

-0. 117 

-0. 003 

39.  0 

5.  066 

-0. 923 

-0. 182 

0.  334 

0.  709 

-0. 092 

0.  001 

40.  0 

5.  156 

-0. 327 

-0. 063 

0.  314 

0.  730 

-0. 126 

-0. 004 

41.  0 

5.  496 

-0. 383 

-0. 070 

0.  316 

0.  762 

-0. 141 

-0. 002 

42.  0 

5.  629 

-1.  066 

-0. 189 

0.  307 

0.  785 

-0.  141 

-0. 002 

43.  0 

5.  705 

0.  051 

0.  009 

0.  290 

0.  796 

-0.  070 

-0. 004 

44.  0 

5.  891 

0.  165 

0.  028 

0.  286 

0.  825 

-0.  125 

-0. 005 

45.  0 

6.  014 

-0. 933 

-0. 155 

0.  268 

0.  834 

-0. 162 

-0. 004 

APPENDIX  D.  DATA  ACQUISITION  PKC^RaM 


The  following  computer  program  was  originally  written  by  Sestak  |Ref.  47]  in 
BASIC  and  modified  by  the  author  for  use  with  this  experiment.  This  program  follows 
a  data  aquisition  shell  program  which  allows  the  data  acquisition  program  to  control  the 
Hewlett-Packard  data  acquisition  hardware.  The  original  program  implemented  auto¬ 
matic  data  input  for  angle  of  attack  and  tare  values.  Angle  of  attack  is  manually  input 
to  this  program  by  the  operator. 

Preliminary  test  runs  indicated  a  larger  than  desired  variation  in  the  balance 
readings.  The  program  was  modified  to  average  ten  separate  readings  per  channel  and 
calculate  the  mean  and  standard  deviation.  Any  readings  which  were  greater  than  one 
standard  deviation  from  the  mean  were  discarded.  The  program  was  also  modified  to 
allow  the  operator  to  recycle  the  balance  reading  routine  without  having  to  leave  the 
program.  This  was  necessary  to  enable  the  operator  to  reject  all  balance  readings  during 
a  current  read  cycle  after  viewing  the  data. 

Since  the  missile  model  pitched  in  the  horizontal  direction,  only  one  tare  reading 
was  necessary,  eliminating  the  original  routine  for  recording  tare  values  at  each  angle 
of  attack.  Additional  information  for  the  program  procedure  is  contained  in  the  the 
program  listing. 
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I 


1000  'Program  to  scan  with  the  DMM  and  RELAY.  MUX.  01 
1010  'This  program  was  writen  by  T.  SESTAK  and  modified  by 
1020  ’P.  RABAN’G  for  use  with  the  VLSAM  model. 
lu3C  ’ 

1040  'This  section  after  the  SHELL  program  directs  reading 
1050  'the  voltages  from  the  balance,  computes  forces  measured 
1060  'by  the  strain  guages ,  then  stores  the  values  in  two  arrays, 

1070  'one  for  the  TARE  one  for  FORCE.  This  data  file  can  then 
1060  'be  used  for  graphs  or  other  displays.  Each  test  run 
1090  'will  generate  a  windtun.dat  file  which  should  be  copied 
1100  'under  another  name  before  the  next  test  run  so  that  it 
1110  'will  not  be  overwritten. 

1120  ' 

1130  'dimension  arrays 

1140  DIM  READING! 7] ,F0RCE[ 140,8] ,TARE[ 8] ,TREAD[ 7 , 10] 

1150  COLOR  14,1,1 
1160  CLS 
1170  ' 

1180  AOA=0 

1190  VALUE=5 

1200  LOCATE  12,28 

1210  prinT'setting  UP  DATA  FILES" 

1220  ' 

1230  'The  program  will  write  the  data  to  several  files. 

1240  STATEFILE$  =  "C:  PCIB  WIND. HPC”  'stored  in  PCIB  subdirectory 

1250  DATAFILE$  =  "C:  WINDTUN.DAT"  'stored  on  disc  C 

1260  DISKFILE?  =  "A:  WINDTUN.DAT"  'stored  on  disc  A 

1270  BALANFILE?  =  "C:  MISSILE  BALANCE. DAT  'stored  on  disc  C 
1280  ' 

1290  RELAY. SETTLING.  TIME  =  .8  ' 800  ms 

1300  LOCATE  16,35:  PRINT'D  0  N  E" 

1310  CALL  DELAY (VALUE) 

1320  ' 

1330  CLS:  LOCATE  12,28:  PRINT"lNITIALIZING  INSTRUMENTS" 

1340  CALL  INITIALIZE.  SYSTEMt  STATEFILE$) 

1350  IF  PCIB. ERR  <>  0  THEN  ERROR  PCIB. BASERR 
1360  CALL  ENABLE.  SYSTEM 

1370  IF  PCIB. ERR  <>  0  THEN  ERROR  PCIB. BASERR 
1380  LOCATE  16 , 35: PRINT'D  0  N  E" 

1390  CALL  DELAY( VALUE) 

1400  ' 

1410  'This  part  of  the  program  is  to  preserve  the  data  if 
1420  'if  the  program  is  aborted  in  mid  run.  Parity  errors 
1430  'in  the  Hewlett  Packard  PC  Instruments  setup  caused  by 
1440  'electrical  noise  and  undervoltage  at  NPS  requires 
1450  'this.  A  voltage  regulated,  uninteruptible  power  source 
1460  'would  ameliorate  this  problem.  Just  in  case-  this  little 
1470  'sequence  allows  reentry  into  the  program  and  the  data 
1480  'arrays  with  minimal  inconvenience. 

1490  ’ 

1500  CLS:  LOCATE  12,20:  INPUT"WERE  YOU  INTERRUPTED  (Y  OR  N)";A$ 

1510  IF  A$="Y"  GOTO  1670 
1520  ’ 

1530  'The  next  three  variables  are  counters  in  the  arrays 
1540  'FORCE  and  TARE 
1550  ’ 
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1560  TRIAL  =  0 
1570  TRY  =  0 
1580  ' 

_ >3C  u^cit  the  datafile  sc  ~ cell  caan  c"r.  be  recorded 

1600  ' 

1610  OPEN  DATAFILE?  FOR  OUTPUT  AS  #1 
1620  CLOSE  if\ 

1630  OPEN  BALANFILF.?  FOR  OUTPUT  AS  #3 
1640  CLOSE  it 3 
1650  GOTO  1800 
1660  ' 

1670  OPEN  DATAFILE?  FOR  INPUT  AS  it  1 

1680  INPUT  #1,  TARE(  1)  ,TARE(  2)  ,TARE(  3)  ,TARE(4)  ,TARE(5)  ,TARF.(  6) , 
TARE(7) ,TARE(8) 

1690  FOR  X  =  1  TO  140 

1700  INPUT  it  1,  FORCE(X, 1) , FORCE (X, 2) ,F0RCE(X,3) ,F0RCE(X,4) , 
FORCE(X,5) ,F0RCE(X,6) ,FORCE(X, 7) ,FORCE(X,8) 

1710  NEXT  X 
1720  CT.OSE  HI 
1730  ’ 

1740  GOTO  1800 

1750  'A?  is  used  as  a  marker  for  interrupted  run  sequences 
1760  'in  the  program,  it  is  set  to  <>"Y"  so  the 
1770  'uninterrupted  sequences  are  used  unless  otherwise  directed 
1780  ’ 

1790  A?="N" 

1800  KEY  OFF 
1810  ' 

1820  'prompt  to  begin  each  scan  or  quit  program  if  desired 
1830  ' 

1840  CLS:  LOCATE  12,10 

1850  INPUT  "TO  START  SCAN  ENTER  ANY  KEY  EXCEPT  Q,  Q  TO  QUIT"; ANSWER? 
I860  IF  ANSWER?  =  "Q"  THEN  GOTO  4730 
1870  ’ 

1880  'THIS  ENTERS  THE  AOA  FOR  EACH  TRIAL  AND  DISPLAYS  IT  ON  THE  CRT 
1890  ' 

1900  CLS:  LOCATE  12,10 

1910  FR1NT  "THE  CURRENT  ANGLE  OF  ATTACK  IS  AOA 
1920  ’ 

1930  ' 

1940  LOCATE  15,10:  PRINT  "INPUT  THE  ANGLE  OF  ATTACK  (AOA)  FOR  THE 
NEXT  TRIAL" 

1950  INPUT  AOA  'READING(l) 

1960  GOTO  1970 
1970  ' 

1980  READING( 1)=A0A 
1990  ' 

2000  'This  variable  is  a  marker  in  the  iteration  loop 
2010  ' interaction  equations  for  convergence. 

2020  ' 

2030  CYCLE  =  0 
2040  ' 

2050  'This  loop  scans  the  pitch  angle  and  6  balance  channels 
2060  'and  stores  the  values  in  the  array  READING 
2070  'Each  channel  is  read  ten  times  and  averaged. 

2080  'The  user  may  reject  the  current  readings  and  input  a  new  set. 


I 


2090 

2100 

2110 

2120 

2130 

2140 

2150 

2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

2290 

2300 

2310 


2320 


2330 

2340 

2350 

2360 

2370 

2380 


CLS 


PRINT"****’" 

>***** *********  j) j RECT 

BALANCE 

READINGS  **** 

PRINT" 

CHECK  OF 

SYSTEM 

OPERATION 

PRINT 
"  IN  VOLTS 

N1 

N2 

SI 

S2  A 

R  " 

PRINT 
"  **■****•** 

****** 

■TWrjWoWf 

‘/rA”jV,;V')Wr 

VnV^WrVrVr’1 

’This  file  is  for  storing  the  direct  voltage  readings  and  averages. 
'The  data  file  is  continually  appended. 

'The  data  is  for  further  analysis  of  the  direct  voltage  readings. 
OPEN  BALANF1LE$  FOR  APPEND  AS  #3 


FOR  CNT  =  1  TO  10 
FOR  CHANNEL  =  2  TO  7 

CALL  OUTPUT! RELAY.  MUX.  01,  CHANNEL) 

IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB. BASERR 
CALL  DELAY! RELAY.  SETTLING.  TIME) 

IF  PCIB. ERR  <>  0  THEN  ERROR  PCIB. BASERR 
CALL  MEASURE! DMM.  01,  READING! CHANNEL] ) 

IF  PCIB.ERR  <>  0  THEN  ERROR  PCIB. BASERR 
TREAD! CHANNEL, CnT)  =  READING! CHANNEL) 


NEXT  CHANNEL 
PRINT  USING 

"  +.######  +.mih m 

READING!  2), READING!  3)  ,READING(4)  , READINGS)  .READING!  6) , READING!  7 ) 
PRINT  #3 ,  USING 

"  +m.  #  +.mnm  +.  mm  +.  wm*  +.######  +.mmm 

READING! 1) .READING! 2) .READING! 3) .READING! 4) .READING! 5) .READING! 6) , 
READING! 7) 

NEXT  CNT 

t 


CALL  SUBROUTINE  TO  AVERAGE  READINGS 
GOSUB  4790 


PRINT" . 

!  I 


2390 


2400 


2410 

2420 

2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 


PRINT  USING 


"MEAN  VALUE  +.#«/#  +.  +.######  +.  MUm  +-M- 

READING! 2 ) , READING! 3 ) , READING! 4) .READING! 5 ) , READING! 6 ) , READING! 7 ) 

PRINT  # 3,  USING 

"  +###.#  +.######  +.######  +.######  +.######"•, 
READING! 1) .READING! 2) .READING! 3) .READING! 4) .READING! 5) .READING! 6) , 
READING! 7) 

CLOSE  #3 
PRINT"  BEEP 

PRINT"<CR>  TO  CONTINUE,  "l"  TO  GET  NEW  READINGS" 

INPUT  XYZ 

IF  XYZ=1  GOTO  2040 

f 


'These  equations  take  voltage  readings  from  the  balance, 
'converts  them  to  counts,  then  applys  the  primary  force 
'equations  to  the  results.  These  values  are  applied  to 
'the  balance  interaction  equations.  Each  channel  has 
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'separate  equations  for  positive  and  negative  readings  and 
'may  have  a  "+"  or  reading  on  any  test  run  so  the 
'rather  involved  logic  path  below  is  my  solution  to  the 
'problem.  For  more  information  consult  Calibration  laboratory 
'guidelines  at  NASA  Ames  Research  Facility  for  TASK  balances 

' **************  CONVERT  SIGNAL  TO  FORCES  *************** 

I 

’Direct  balance  readings  are  multiplied  by  a  scale  factor 
'5000000  then  divided  by  the  balance  exitation  voltage  to 
'get  a  readiing  in  COUNTS.  The  program  will  send  each  reading 
'to  the  appropriate  equation  and  convert  to  force  or  moment 
'then  return  to  send  the  next  reading  for  calculation 
'The  data  acquisition  system  for  using  this  program  used  an 
'amplifier  with  1000  gain.  The  scale  factor  is  divided  by  1000. 


VEX=5 

N1=READING(2)*5000!  /VEX 
N2=READING!3)*5000!  /VEX 
S 1=READING( 4)*5000!  /VEX 
S2=READING( 53*5000!  /VEX 
A=READING( 6 )*5000!  /VEX 
R=READING(  73*416.  67,>VVEX 


'Excitation  Voltage 


send  each  reading  to  the  appropriate  equation 


IF  READING(23>0 
IF  READING(3)>0 
IF  REaDING! 4)>C 
IF  READING! 5 ) >0 
IF  READING! 6)>0 
IF  READING! 7 ) >0 


2870  ELSE 
2S3G  ELSE 
2930  ELSE 
2950  ELSE 
2910  ELSE 
2970  ELSE 


GOTO  3020 
GOTO  3040 
GOTO  3080 
GOTO  3100 
GOTO  3060 
GOTO  3120 


' **************  POSITIVE  FORMULAS  ************* 

I 

EN1  =  . 050861*N1  -  5. 4826E-09*(N1*N1) 

GOTO  2790 

EN2  =  .  047211--N2  -  1.  7015E-08*! N2*N23 
GOTO  2800 

EA  =  . 014309*A  -  7. 1962E-10*! A*A) 

GOTO  2830 

ESI  =  . 031309*S1  -  3. 8153E-08*! 31*S1) 

GOTO  2810 

ES2  =  . 030366*S2  -  3. 8607E-08*! S2*S23 
GOTO  2820 

ER  =  . 0030885*R  +  2. 5672E-09*! R*R3 
GOTO  3130 

t 

***************  NEGATIVE  FORMULAS  ************* 


EN1  =  . 051591*N1  +  1. 7157E-08*(N1*N1) 
GOTO  2790 

EN2=.  047763*N2+8.  915299E-09*(N2*N2) 
GOTO  2800 

EA  =  . 01429*A  -  1. 3322E-09*(A*A3 


3070  GOTO  2830 

3080  ESI  =  . 032073*S1  -  8.  931601E-09*(S1*S1) 

30°0  GOTO  2310 

3100  ES2  =  . 031167*S2  -  7. 2517E-09*! S2*S2) 

3110  GOTO  2820 

3120  ER  =  .  0030908--R  -  2. 4769E-09*!R*R) 

3130  ' 

3140  ’ 

3150  'a  heading  for  the  iteration  values 
3160  ' 

3170  PRINT"  " 

3180  PRINT" ********************  FORCE  INTERACTION  ITERATIONS 
****************** 

3190  PRINT"  CHECK  FOR  CONVERGENCE 

3200  PRINT 

"  CYCLE  AOA  N1  N2  SI  S2  A  R  " 

3210  PRINT 

"  #  DEG  POUNDS  POUNDS  POUNDS  POUNDS  POUNDS  FT-LBS 

3220  PRINT 

**  VciWfiWrVf  Vr-iV  ,V  .WrsV  VntoWnWc  sWftfftWwV  VoWr  V'*V:V  *W,'V?:W«V 

3230  ’ 

3240  'The  loop  that  controls  the  balance  interaction 
3250  'equations  and  allows  a  visual  convergence  check 
3260  ’ 

3270  FOR  I  =  1  TO  10 

3280  IF  READING(2)>0  THEN  GOTO  3370  ELSE  GOTO  3570 

3290  IF  READING! 3)>0  THEN  GOTO  3400  ELSE  GOTO  3600 

3300  IF  READING! 4 )>0  THEN  GOTO  3460  ELSE  GOTO  3660 

3310  IF  READING! 5 )>0  THEN  GOTO  3490  ELSE  GOTO  3690 

3320  IF  READING! b)>0  THEN  GOTO  3430  ELSE  GOTO  3630 

3330  IF  READING! 7 )>0  THEN  GOTO  3520  ELSE  GOTO  3720 

3340  1 

3350  ' ************** POSITIVE  FORMULAS************** 

3360  ' 

3370  XN1=  EN1+. 005e036*N2+. 0041655*S1+. 058079*R-7. 1926E-07*(N2*N2) 

+4. 0352E-06*! S1*S1)~. 0006786*!R*R) 

3380  GOTO  3290 
3390  ' 

3400  XN2=  EN2+.  046218*N1-.  0028393*A-. 0081694*S1+. 0041463*S2+. 077279 

*R-6.  8577E-07*!N1*N1) -1.  7755E-05*(A*A)+2. 1719E-06*! S1*S1)+1. 8582E-06 
* !  S  2*  S  2 )  - .  0  0 1 9  2  9  4*  ( R*R  ) 

3410  GOTO  3300 
3420  ' 

3430  XA=  EA+8. 6893E-04*Nl+6. 0359E-04*Sl+7. 7722E-05*S2-. 11115*R 

+4. 45  3  7E - 0  7* ! N 1*N 1 ) +4 . 7936E-06*!S1*S1) -4.  1033E-06*!S2*S2)+2. 0697E-04 
*!R*R) 

3440  GOTO  3330 
3450  ' 

3460  XS1  =  ES1-6. 3459E-04*N1-. 11148*R+5. 535E-06*! N1*N1)+. 0024592*! R*R) 
3470  GOTO  3310 
3480  ' 

3490  XS2=  ES2-.  0024237*N1+.  0022455*A+. 0066785*S1-. 26377*R+1. 7099E-06 
*!N1*N1)+1. 2072E-05*! A*A) -2. 7825E-06*! S1*S1)+. 0062217*! R*R) 

3500  GOTO  3320 
3510  ' 

3520  XR=  ER- 1.  9928E-04*N2-2.  5893E-04*S2+1. 1512E-07*(N2*N2)-5.  156E-08 


*(S2*S2) 

GOTO  3730 

I 

****************  negative  formulas  ************ 

I 

XN'1=  EN1+.  010257*N2-.  0045396*S1  04494*R+7.  9499E-07*(N2*N2) 

-1.  967E-06,V(S1*S1)  -.  0003232*(R*R) 

GOTO  3290 

* 

XN2=  EN2+. 05 1778*N1-. 0044056*A-9. 038499E-03*S1-. 061125*R 

+5. 2897E-06*(N1*N1)+1.  0467E-05*(A*A) -4.  8493E-07*(S1*S1)-.  0011773*(R*R) 

GOTO  3300 

i 

XA=  EA-. 0021217*Nl+9. 1524E-04*N2-. 097l48*R-4. 2547E-06*(N1*N1) 

+4. 5846E-06*(N2*N2)-7. 5001E-04*(R*R) 

GOTO  3330 

i 

XS1=  ESI-. 0071275*N1-. 0089235*A-. 05268*R-1. 2923E-05*(N1*N1) 

-4. 0345E-05*(A*A)-9. 3969E-04*(R*R) 

GOTO  3310 

t 

XS2=  ES2-. 0037176*N1-. 0052619*N2+. 0072915*A+.  006856*S1-.  062581 

*R - 5 . 2 1 1 E - 0  7  * ( N 1 * N 1 ) - 8 . 6265E-06*(N2*N2)+3. 7054E-05*(A*A)+9. 963001E-06 

*( S  1*S 1 )+8. 0007E-04*(R*R) 

GOTO  3320 

i 

XR=  ER+3. 5945E-04*N1+1.  5497E-07*(N1*N1) 

i 

'Shift  all  the  new  variables  back  to  the  old  name 

'for  the  next  iteration 

N1=X.N1 

N2=XN2 

A=XA 

S1=XS1 

S2=XS2 

R=XR 


'A  marker  for  the  interations 
CYCLE  =  CYCLE  +  1 

'print  the  iterations  to  watch  for  convergence 
» 

PRINT  USING 

"  ##  +##. ft  Htftft.  ftft  Hfftft.  ftft  +ftftft.  ftft  +ftftft.  ftft  Hftft.  ftft  +ftftft.  ftft"\ 

CYCLE , A0A , N 1 , N2 , S 1 , S2 , A , R 

NEXT  I 
» 

INPUT  "IF  CONVERGENCE  IS  ADEQUATE  ENTER  Y,  IF  ANOTHER  RUN 
IS  DESIRED  ENTER  N";ANSWR$ 

IF  ANSWR?  =  ”N"  THEN  GOTO  3160 

t 

NORMAL  =  N1  +  N2 
SIDE  =  SI  +  S2 
AXIAL  =  A 

PITCH  =  (N1-N2)  *  .0854 
ROLL  =  (S1-S2)  *  .  0698 
YAW  =  R*12 


3990  ' 

4000  TRIAL  =  TRIAL  +  1 

4010  INPUT;  "IS  THIS  A  TARE  READING,  Y  OR  N"; AN$ 

4020  IF  AN$  <>  "Y"  GOTO  4290 
4030  COLOR  10,4,1 
4040  TRIAL  =  0 
4050  TRY  =  TRY  +  1 
4060  IF  A$<>"Y"  GOTO  4100 
4070  CLS:  LOCATE  12,10 

4080  INPLT"EN'TER  THE  VALUE  OF  THE  LAST  TRY  BEFORE  INTERRUPT";  ITRY 
4090  TRY  =  ITRY  +  1 
4100  TARE(l)  =  TRY 
4110  TARE( 2)  =  AOA 
4120  TARE( 3)  =  NORMAL 
4130  TARE(4)  =  SIDE 
4140  TARE(5)  =  AXIAL 
4150  TARE(6)  =  PITCH 
4160  TARE( 7 )  =  ROLL 
4170  TARE( 81  =  YAW 
4180  ’ 

4190  PRINT"  " 

4200  PRINT”*  *  *  *  *  * 

■jV  Vr  Vr  *  Vr  *  :Y  -.V  * 

4210  PRINT 

"  TRIAL  AOA 

4220  PRINT 

"  il  DEG 

4230  PRINT 

M  VrVrVfjVVr  Vri'.-sWr-v'c 

4240  'A  loop  to  list  all  values  so  far 
4150  FOR  J  =  1  TO  TRY 

4260  PRINT  USING"  M  +M.M  Mil.  ilil  Mil.  ilil  Mil.  M  Mil.  M 
Mil ■  M  Mil •  M" ;  TARE ( 1 ) , TARE ( 2) , TARE ( 3 ) , TARE ( 4 ) , TARE ( 5 ) , TARE (  6 ) , 

TARE(  7  )  ,TARE(  8) 

4270  NEXT  J 

4280  GOTO  4550 

4290  ' 

4300  IF  A$<>"Y”  GOTO  4340 
4310  CLS:  LOCATE  12,10 

4320  INPUT"ENTER  THE  VALUE  OF  THE  LAST  TRIAL  BEFORE  INTERRUPT";  ITRIAL 

4330  TRIAL  =  ITRIAL  +  1 

4340  FORCE(TRIAL, 1 )  =  TRIAL 

4350  FORCE(TRIAL, 2)  =  AOA 

4360  FORCE(TRIAL, 3)  =  NORMAL  -  TARE( 3) 

4370  F0RCE(TRIAL,4)  =  SIDE  -  TARE(4) 

4380  FORCE(TRIAL,5)  =  AXIAL  -  TARE( 5) 

4390  FORCE(TRIAL, 6)  =  PITCH  -  TARE( 6) 

4400  FORCE(TRIAL, 7 )  =  ROLL  -  TARE( 7 ) 

4410  FORCEfTRIAL, 8)  =  YAW  -  TARE( 8 ) 

4420  ' 

4430  'print  the  values  and  store  in  file 
4440  ' 

4450  PRINT"  " 

4460  PRINT"*  *  *  *  *  *  *  *  *  *  *  *  FORCE  CALCULATIONS  *  *  * 

Vr  Vr  V?  Vr  Vc  *  ■>*  Vr  * 

4470  PRINT 


*  *  *  *  *  *  TARE  CALCULATIONS  *  *  *  * 


NORMAL 

SIDE 

AXIAL 

PITCH 

ROLL 

YAW" 

POUNDS 

POUNDS 

POUNDS 

FT-LBS 

FT-LBS 

FT-LBS" 

Vf-.VVf^VVfVf 

*Vr>VVcVc**V 

iWrVr«V“V*jV 

IV. 


p 


1 

"  TRIAL 

A0A 

NORMAL 

SIDE 

AXIAL 

PITCH 

ROLL 

YAW" 

o 

oo 

PRINT 

■ 

"  JL 

ir 

DEG 

POUNDS 

POUNDS 

POUNDS 

FT-LBS 

FT-LBS 

FT-LBS 

■  4490 

PRINT 

,f  Vr-VbV-VVf 

***** 

iWnVsV^V 

****** 

****** 

4500 

4510 

' a  loop  to 

list 

all  values 

so  far 

4520  FOR  J  =  1  TO  TRIAL 

4530  PRINT  USING"  M  +###.#  +###.##  +###.##  +###.##  +###.## 

_  +### .  ##  +### .  #5?" ;  FORCE ( J  ,  1 )  , FORCE ( J ,  2 )  , FORCE ( J ,  3 ) , FORCE(  J  ,4) , 

■  FORCE( J,5) ,FORCE( J,6) ,FORCE(J,7) ,FORCE( J,8) 

4540  NEXT  J 
4550  ' 

4560  'Write  the  data  to  the  output  data  files 
4570  COLOR  14,1,1 

4580  OPEN  DATAFILES  FOR  OUTPUT  AS  #1 

■4590  WRITE  #1,  TARE( 1) ,TARE(2) ,TARE(3) ,TARE(4) ,TARE(5) ,TARE(6) , 

TARE( 7) ,TARE( 8) 

4600  FOR  X  =  1  TO  140 

4610  WRITE  #1,  FQRCE(X , 1 ) , FORCE(X , 2) .FORCE! X , 3 ) , FORCE(X , 4) , 

FORCE  ( X ,  5 ) ,  FORCE  ( X ,  6 ) ,  FORCE ( X ,  7 ) ,  FORCE  ( X ,  8 ) 

4620  NEXT  X 
4630  CLOSE  #1 
^  4640  ’ 

4650  OPEN  DISKFILE$  FOR  OUTPUT  AS  #2 

4660  WRITE  # 2,  TARE( 1) ,TARE(2) ,TARE(3) ,TARE(4) ,TARE(5) ,TARE(6) , 
TAREC  7 ) ,TARE( 8 ) 

4670  FOR  X  =  1  TO  140 

4680  WRITE  ,v 2,  FORCE(X,  1)  ,F0RCE(X, 2)  ,FORCE(X , 3)  .FORCECX ,4) , 
FORCECX , 5 ) , FORCECX , 6 ) , FORCECX , 7 ) , FORCECX , 8 ) 

■  4690  NEXT  X 

^  4700  CLOSE  n 

4710  ' 

4720  'Prompt  for  next  scan 
4730  ’ 

4740  INPUT  "DO  YOU  WANT  ANOTHER  SCAN  CY  OR  N)";ANSW$ 

B  4750  IF  ANSWS  <>"N"  THEN  GOTO  1790 

■  4760  ' 

4770  END 
4780  ' 

4790  'This  subroutine  averages  the  balance  voltage  readings 
I  4800  'bv  computing  the  mean  and  standard  deviation. 

|  4810  'Any  readings  less  or  greater  than  one  standard  deviation 

4820  'are  thrown  out  and  a  new  mean  is  computed 
4830  ' 

;  4840  FOR  CHANNEL  =  2  TO  7 

4850  N=10 
4860  SREAD=0 
4870  SSDEV=0 

#  4880  ' 

4890  'Mean  of  balance  voltage  readings 

4900  FOR  CNT  =  1  TO  10 

4910  SREAD  =  SREAD  +  TREADC CHANNEL, CNT) 

4920  NEXT  CNT 

4930  MEAN  =  SREAD /N 

4940  READING (CHANNEL)  =  MEAN 


137 


10)  THEN  GOTO  5120 


'Standard  deviation  routine 
FOR  CNT  =  1  TO  10 
DIF  =  TREAD( CHANNEL, CNT)  -  MEAN 
SDEV  =  DIF  *  DIF 
SSDEV  =  SSDEV  +  SDEV 
NEXT  CNT 

DEV  =  SQR( SSDEV/N) 

HI  =  MEAN  +  DEV 

LO  =  MEAN  -  DEV 

FOR  CNT  =  1  TO  10 

ARG  =  TREAD (CHANNEL, CNT) 

IF  (ARG  <  HI)  AND  (ARG  >  LO)  THEN  GOTO  5110 
TREAD( CHANNEL, CNT)  =  0! 

N  =  N  -  1 
NEXT  CNT 
NEXT  CHANNEL 
RETURN 
END 
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APPENDIX  E.  COEFFICIENTS  TRANSLATION  PROGRAM 


This  program  was  written  by  the  author  for  use  with  this  thesis  experiment.  The 
program  reads  a  data  file  containing  force  and  moment  values  and  converts  these  into 
a  coefficient  form  taking  into  account: 

Test  conditions. 

Body  configuration  for  blockage. 

Screen  number. 

The  results  are  written  into  a  data  file  which  are  used  to  produce  graphs  and  plots. 


> 


1000 

1010 

1020 

|  1030 

1040 
1050 
1060 
1070 
1080 

.  1090 

I  1100 

1110 
1120 
1130 
1140 
1150 

I  1160 

1  1170 

1180 
1190 
1200 
1210 
1220 

i  1230 

1240 
1250 
1260 
1270 
1280 

i  1290 

1300 
1310 
1320 
1330 
1340 
1350 
1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 


’  PROGRAM  BY  M. P.  RABANG  TO  READ  FORCE  AND  MOMENT  VALUES  FROM 
'  A  DATA  FILE  CREATED  BY  THE  DATA  ACQUISITION  PROGRAM. 

COLOR  10,4,1 

DIM  TARE[ 8] ,  FORCE[ 140,8] ,  COEF{ 140,9] 

CLS 

LOCATE  5 

INPUT" ENTER  THE  NAME  OF  THE  INPUT  FILE": INFILE$ 

INPUT"ENTER  THE  NAME  OF  THE  OUTPUT  FILE";  OUTFILE$ 

INPUT"ENTER  THE  MINIMUM  TEMPERATURE  (F)";TMIN 
INPUT"ENTER  THE  MAXIMUM  TEMPERATURE  (F)";TMAX 
INPUT"ENTER  THE  BODY  TYPE,  A=+,  B=0,  C=X  BOD$ 

IF  BOD$="A"  THEN  BODY=l:  GOTO  1180 
IF  BOD$="B"  THEN  B0DY=3:  GOTO  1180 
IF  BOD$=”C"  THEN  BODY=2:  GOTO  1180 
GOTO  1130 

INPUT" ENTER  THE  SCREEN  NUMBER.  0-4"; SCR 
IF  SCR=0  THEN  GOTO  1250 
IF  SCR=1  THEN  Q=15. 35:  GOTO  1260 
IF  SCR=2  THEN  Q=14. 882:  GOTO  1260 
IF  SCR=3  THEN  Q=16. 38:  GOTO  1260 
IF  SCR=4  THEN  Q=15. 605:  GOTO  1260 
GOTO  1180 

INPUT"ENTER  THE  WIND  TUNNEL  VELOCITY  (cm  H20)";PH20 
INPUT"DO  YOU  WANT  TO  VIEW  OUTPUT"; PANS? 

INPUT”DO  YOU  WANT  A  HARDCOPY"; ANS$ 

’  OPEN  THE  DATA  FILE  SO  EACH  SCAN  IS  RECORDED 

TRANSFILE$="C:  MISSILE  TRANS.  DAT" 

OPEN  OUTFILES  FOR  OUTPUT  AS  #2 

CLOSE  if! 

\ 

'  READ  THE  FORCE  VALUES  FROM  THE  DATA  FILE 
OPEN  INFILES  FOR  INPUT  AS  if 1 

INPUT  if  1,  TARE ( 1 )  , TARE ( 2 ) , TARE ( 3 ) , TARE ( 4 ) , TARE ( 5 ) , TARE (  6 ) , 
TARE( 7) ,TARE( 8) 

FOR  X  =  1  TO  140 

INPUT  ifl,  F0RCE(X,  1) , FORCE (X, 2)  , FORCE (X, 3) , FORCE (X, 4) , 
FORCE(X,5) ,FORCE(X,6) ,FORCE(X,7) ,FORCE(X,8) 

IF  FORCE(X, 1 ) =0 !  THEN  GOTO  1410 
NEXT  X 
CLOSE  in 

IF  PANS$<>"Y"  THEN  GOTO  1610 
! 

'  DISPLAY  FORCE  VALUES 

COLOR  14,1,1 

CLS 

PRINT"FILE  NAME:  "; INFILES 
PRINT"  " 

PRINT"*  *************  FORCE  READINGS 

■>V  *  *  Vr  *  Vf  Vc  *  *  *  *  Vr  -iV  *  *  ->Y  * 

PRINT”  ” 

PRINT 

"  TRIAL  AOA  NORMAL  SIDE  AXIAL  PITCH 


140 


ROLL  YAW" 

1520  PRINT 

,f  jl  nrr 


"  #  DEG  POUNDS 

FT-LBS  FT-LBS" 

. POUNDS 

POUNDS 

FT-LBS 

1530 

PRINT 

******** 

****** **  ******** 1 1 

1540 

1550 

t 

FOR  J  =  1  TO  140 

1560 

IF  FORCE! J,1)=0  THEN  GOTO 

1590 

1570 

PRINT  USING"  Mil  +###.# 

+##.  MM 

+m.  mm 

+M-  MM 

+##.  MM  +M-  MM  +##.  MM";  FORCE! J,  1)  ,F0RCE( J,2)  .FORCE!  J,3) , 

F0RCE(  J ,4) ,F0RCE( J , 5 ) , F0RCE( J , 6 ) , F0RCE( J , 7 ) ,F0RCE( J , 8 ) 

1580  NEXT  J 
1590  ' 

1600  INPUT"ENTER  <CR>  TO  CONTINUE"; INPT$ 

1610  ' 

1620  ’  BEGIN  COEFFICIENTS  CALCULATION 
1630  A=. 0167 
1640  MU=3. 719E-07 
1650  TAV=(TMIN+TMAX)/2! 

1660  RHO=2116.  22/(  1545*(459. 7+TAV)) 

1670  IF  SCR>0  GOTO  1740 

1680  VE=( 2!  *2.  046*PH20)/(.  93--RHO) 

1690  VEL=SQR( VE) 

1700  RED=( RHO*VEL*( 1.  75/12))/MU 
1710  Q=(RHO*VEL*VEL)/2! 

1720  GOTO  1800 
1730  ' 

1740  VE=2*Q/RH0 
1750  VEL=SQR( VE) 

1760  RED=(RH0*VEL*(1. 75/12) )/MU 
1770  PH20=Q*. 93/2. 046 
1780  ' 

1790  '  WRITE  THE  COEFFICIENTS  TO  THE  OUTPUT  FILE 
1800  OPEN  OUTFILE5  FOR  OUTPUT  AS  1U 
1810  OPEN  TRANSFILES  FOR  APPEND  AS  ji 3 
1820  WRITE  #3,  0UTFILE$ 

1830  FOR  X  =  1  TO  140 

1840  FOR  Y  =  3  TO  8 

1650  IF  FORCE! X,1)=0  THEN  GOTO  2040 

1860  '  ROUTINE  TO  CORRECT  THE  DYNAMIC  PRESSURE  FOR  BLOCKAGE 
1870  '  AND  BODY  TYPE 
1880  COEF(X, l)=FORCE(X, 1) 

1890  COEF( X, 2)=FORCE(X, 2) 

1900  ALPHA  =  FORCE(X,2) 

1910  IF  FORCE(X, 2)  <  0  THEN  ALPHA  =  ABS(FORCE(X,2) ) 

1920  IF  FORCE(X, 2)  >  90  THEN  ALPHA  =  180-F0RCE(X,2) 

1930  IF  BODY=l  THEN  EPS=. 0000126*ALPHA+. 007759 
1940  IF  BODY=2  THEN  EPS=.  0000 101*ALPHA+. 007759 
1950  IF  BODY=3  THEN  EPS=.  0000908*ALPHA+.  007759 
1960  D  =  A*Q*( l+( 2*EPS) ) 

1970  COEF(X,Y)  =  FORCE(X,Y)/D 
1980  C0EF(X,9)=C0EF(X,4)/C0EF(X, 3) 

1990  NEXT  Y 
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I 


I 


2000 

2010 

2020 

2030 

2040 

2050 

2060 

2070 

2080 

2090 

2100 

2110 

2120 

2130 

2140 

2150 

2160 

2170 

2180 

2190 

2200 

2210 

2220 


2230 


2240 

2250 

2260 


2270 

2280 

2290 

2300 

2310 

2320 

2330 

2340 

2350 

2360 

2370 

2380 

2390 

2400 

2410 

2420 

2430 

2440 

2450 


IF  F0RCE(X,1)=0  THEN  GOTO  2040 

WRITE  #2,  C0EF(X, 1) ,COEF(X,2) ,COEF(X,3) ,C0EF(X,4) ,COEF(X,9) , 

CO£F(X,5) ,COEF(X,6) ,C0EF(X,7) ,COEF(X,8) 

WRITE  #3,  COEF(X, 1) ,COEF(X,2) ,COEF(X,3) ,C0EF(X,4) ,COEF(X,9) ,  | 

COEF ( X , 5 ) , COEF ( X ,  6  )  , COEF ( X ,  7  )  , COEF ( X ,  8 ) 

NEXT  X 
CLOSE  in 
CLOSE  #3 

IF  PANS$<>"Y"  THEN  GOTO  2300 

f 

'  DISPLAY  ROUTINE  • 

COLOR  10,4,1 
CLS 

PRINT"FILE  NAME:  ";OUTFILE$ 

PRINT"  " 

PRINT" AVERAGE  TEMPERATURE  (F) 

PRINT'VlND  TUNNNEL  VELOCITY  (FT/SEC) 

PRINT"WIND  TUNNNEL  VELOCITY  (cmH20) 

PRINT"AIR  DENSITY  ( LBM/FT3) 

PRINT"REYNOLDS  NUMBER 
PRINT"DYNAMIC  PRESSURE  ( LB/FT2) 

PRINT"  " 

PRINT"*  *  *  *  *  *  ••  ******  J.-ORCE  COEFFICIENTS 


*  *  *  *  * 

“it  “it  'it  “it 

*  *  *  * 

*  * 

PRINT"  " 
PRINT 
"  TRIAL 
ROLL 

PRINT 

AOA 

YAW" 

NORMAL 

SIDE 

AXIAL 

PITCH 

******** 

****** 

**M 

FOR  X  =  1  TO  140 

IF  COEF(X,3)=0  THEN  GOTO 
PRINT  USING"  Uii  +##?;.  it 

2280 

+##. #### 

+m.  un 

+##.  ##>># 

+##.  MM  +##.  MM  +M.  MM" ;  COEF ( X ,  1 ) , COEF ( X , 2 )  ,COEF(X,3) , 

C0EF(X,4) , COEF ( X , 5 ) , COEF ( X , 6 ) , COEF ( X , 7 ) , COEF ( X , 8 ) 

NEXT  X 

I 

PRINT"  " 

IF  ANS$<>"Y"  THEN  GOTO  2640 
'  HARDCOPY  ROUTINE 
LPRINT" INPUT  FILENAME:  "; INFILE$ 

LPRINT" OUTPUT  FILENAME:  ";  OUTFILE$ 

LPRINT"  " 

LPRINT"INITIAL  TEMPERATURE  (F) 

LPRINT"FINAL  TEMPERATURE  (F) 

LPRINT" AVERAGE  TEMPERATURE  (F) 

LPRINT"WIND  TUNNEL  VELOCITY  (FT/SEC) 

LPRINT" WIND  TUNNEL  VELOCITY  (cmH20) 

LPRINT" AIR  DENSITY  ( LBM/FT3) 

LPRINT"REYNOLDS  NUMBER 
LPRINT"DYNAMIC  PRESSURE  ( LB/FT2) 

LPRINT"  " 

LPRINT"  " 

LPRINT"*  *  *  *  *  *  *  *  *  *  *  *  *  *  FORCE  READINGS 


;  TMIN 
";  TMAX 
";  TAV 
VEL 
PH20 
“ ;  RHO 
RED 
Q 


VEL 

PH20 

RHO 

RED 

Q 
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POUNDS 


PITCH 


FT-LBS 


*  *  it  it  it  it  it  it  it  *  *  *  *  it  it  it  it 

2460  LPRINT"  " 

2470  LPRINT 

"  TRIAL  AOA  NORMAL  SIDE  AXIAL  PITCH 

ROLL  YAW" 

2480  LPRINT 

"  #  DEG  POUNDS  POUNDS  POUNDS  FT-LBS 

FT-LBS  FT-LBS" 

2490  LPRINT 

"  *****  *****  ********  ********  ********  ******** 
********  ********" 

2500  FOR  J  =  1  TO  140 

2510  IF  FORCE( J, 1)=0  THEN  GOTO  2540 

2520  LPRINT  USING"  +###.#  HU.  MM  Hit.  MM  +M.M1W 

4##.  ####  +##.  ####  HU.  Him"-,  FORCE(  J ,  1) ,  FORCE(  J ,  2 ) ,  FORCE  (J ,  3 ) 

FORCE( J,4) ,F0RCE( J,5) ,F0RCE(J,6) ,FORCE( J,7) ,F0RCE(J,8) 

2530  NEXT  J 
2540  LPRINT"  " 

2550  LPRINT"  " 

2560  LPRINT"*  *  *  *  *  *  *  *  *  *  *  *  *  *  FORCE  COEFFICIENTS 

Ve  *  -V  Vr  it  it  it  it  it  it  it  9V  it  it  it 

2570  LPRINT"  " 

2580  LPRINT 

"  TRIAL  AOA  NORMAL  SIDE  AXIAL  PITCH 

ROLL  YAW" 

2590  LPRINT 

"  *****  ****  ********  ********  ********  ******** 
********  ******** 1 1 

2600  FOR  X  =  1  TO  140 
2610  IF  C0EF(X , 3)=0!  THEN  GOTO  2640 

2620  LPRINT  USING"  i!M  HiHt.it  HU.  MM  HU.  MM  Hit.  MM 

Hi?.  MM  HU.  MM  Hit.  MM"',  C0EF(X,  1)  ,C0EF(X,2)  ,COEF(X,  3) , 
C0EF(X,4) , COEF ( X ,  5 ) , COEF ( X , 6 ) ,COEF(X,7) ,C0EF(X,8) 

2630  NEXT  X 
2640  ' 

2650  INPUT"D0  YOU  WANT  TO  EXIT  THE  PROGRAM"; AANS$ 

2660  IF  AANS$<>"Y"  THEN  GOTO  1060 
2670  END 
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